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PREFACE 


Following a decade of research in the application of remote 
sensing technology to earth resources measurement and monitoring 
with earth observing satellites, the National Aeronautics and 
Space Administration found a need developing for the advancement 
of our understanding of the physical processes associated with 
the remote sensing. Therefore, at the beginning of the 1980s 
NASA established a continuing research program in fundamental 
remote sensing science that would enhance our ability to learn 
moT-e about the physical and biological processes of our planet 

and would contribute to a base of 
needed to support the planning of new 
flight programs. 

Characterization 
established within this program to improve 
the fundamental relationships of energy interactions between the 
sensor and the surface target, including the effect of the 
atmosphere . 


from space acquired data 
scientific understanding 
and effective sensors and 
and Atmospheric Effects 


The Scene Radiation 
(SRAEC) Project was 
our understanding of 


This report has been prepared to provide a very brief 
overview of the SRAEC Project history and objectives and to 
report on the scientific findings and project accomplishments. 
An attempt has been made to elucidate the overall project 
progress and accomplishments through principal investigator team 
introspection, and the principal investigators have prepared 
succinct summaries of their individual study results achieved 
since the project’s initiation just over three years ago. 


A meeting of the Pis is conducted annually to bring the 
scientists together to present and discuss their results with one 
another. This annual summary report derives from the most recent 
annual meeting held at the Massachusetts Institute of Technology 
in Cambridge, Massachusetts during January 29-31, 1985. The 
overall program management is under the direction of Dr. Robert 
E. Murphy of NASA Headquarters, Code EE, Washington, D.C. 


Donald W. Deering 
GSFC Project Science Manager, 
Fundamental Remote Sensing Science 
Research Program 

Earth Resources Branch (Code 623) 
Goddard Space Flight Center 
Greenbelt, MD 20771 
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I. SRAEC PROJECT BACKGROUND 
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I. A. SRAEC PROJECT ORIGIN AND OBJECTIVES 


During the "Landsat Era” of the 1970* s, the development of 
rapid, empirical approaches to satellite data interpretation took 
precedence over the development of a full scientific understand- 
ing of the earth-atmosphere scene radiation complex. The 
importance of synoptic remote measurements of the earth’s 
surface and the atmosphere has been clearly demonstrated, but the 
need to more fully utilize the powerful satellite remote sensing 
measurement potential has led to the realization that purely 
empirical techniques are no longer sufficient. 

Significant advancements in both s:eaJote «>ensing components 
and systems and in the remote sensing derived underrtanding of 
physical and biological processes on the eart.i’s surface require 
increased understanding of the characteristics of earth scene 
radiation and atmospheric effects at visible, infrared and 
microwave wavelengths. Thus, NASA established the Fundamental 
Remote Sensing Science Research Program within the Land Processes 
Branch of the Earth Science and Applications Division at NASA 
Headquarters . 

The Scene Radiation and Atmospheric Effects Characterization 
Project is conducted as the primary element of the program with 
the general goal being to improve our understanding of the energy 
emitted or reflected from an earth surface target, through the 
intervening atmosphere as measured by a remote sensing system. 
Emphasis has been placed on theoretical model development with 
some observational studies and empirical characterizations 
relating electromagnetic energy reflected and emitted from earth 
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surfaces to their biophysical characteristics. 

Fundamental research needs for future applications research 
in aerospace remote sensing were identified through a series of 
workshops during fiscal year 1980 that involved leading scien- 
tists from universities, industry and government. The research 
issues that were identified through these workshops and reported 
as a final report to NASA Johnson Space Center (under NAS 
9-16016) in August 1980 are summarized in last year’s annual 
SRAEO report (NASA Technical Memorandum 86078, March 1984) and, 
therefore, are not repeated here. 

Proposals from organizations outside NASA were solicited in 
July 1981, and following a NASA Headquarters conducted peer 
review, sixteen investigations were funded during fiscal year 
1982. Most of these were proposed as three-year studies, such 
that this annual report reflects the end-of-study results for 
many of the principal investigators. Since the initial call for 
proposals, other unsolicited proposals have been received and 
peer reviewed. Therefore, there are a few principal investi- 
gators who have only recently begun their studies under this 
project. Even though there is finite funding available, pro- 
posals for new research are continually encouraged, especially in 
those topical areas that were identified in the early planning 
phase as being important research issues that have yet to be 
adequately addressed. 

The product of the SRAEC project is an improved understand- 
ing of the physical radiant energv interaction (remote sensing) 
processes through measurements and modeling that will ultimately 
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enable the proper specification of earth observing space sensor 
systems and the measurement and modeling of the important 
physical, chemical, and biological processes and their inter- 
actions. The success of the individual scientists in the program 
is measured by the quality professional papers that are produced 
and contributed through refereed scientific journals. 
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I.B. PRINCIPAL INVESTIGATOR WORKSHOPS 


Following the selection and announcement of successful SRAEC 
proposals, a "kick-off" meeting was held on June 15, 1982 at NASA 
Goddard Space Flight Center. The objective of this meeting was 
to enable the principal investigators in the SRAEC Project to 
i8£ a V at their new collegues in the project were planning to do 
in their various studies and begin the process of sharing 
iiit >r.fflation and ideas for synergistic benefit to the remote 
sensing science program. Collaboration among principal investi- 
gators was encouraged. 

The next meeting of the SRAEC principal investigators was 
held January 5-7, 1983 at the University of Colorado in Boulder 
in two separate sessions as a part of the National Radio Science 
Meeting. The objectives for this gathering of the Pis was to 
provide NASA management with a review of early progress, to 
initiate the sharing of results between SRAEC investigators and 
to emphasize the importance of reporting the study results with 
the broader remote sensing scientific community. 

The second meeting to report on the research results was 
conducted at Colorado State University in Fort Collins from 
January 9-11, 1984. This second annual meeting was conducted 
with a workshop format and to serve as a forum for sharing 
terminology and articulate the merging of visible, IR, and 
microwave theoretical and experimental techniques. Tutorial 
presentations were made by representatives of these remote 
sensing disciplines. Although time constraints did not allow a 
full development of the ideas, good progress was made in bridging 
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the communications gap between the groups, which has been 
apparent in the most recent meeting. 

The Third Annual Meeting of the SRAEC principal investi- 
gators was held on the campus of the Massachusetts Institute of 
Technology in Cambridge, Massachusetts during January 29-31, 
1985. Dr. Jin Kong of the Research Laboratory of Electronics at 
MIT hosted this meeting and arranged for local scientists Ken 
Hardy and Dave Staelin to address the Pis and stimulate their 
thinking. 

The letter of invitation to the principal investigators and 
the meeting agenda are included in Appendix A. A list of the 
workshop participants, their affiliations and addresses is 
provided in Appendix B. And a telecoaaunications message to the 
meeting participants from Dr. Robert Murphy, Fundamental Remote 
Sensing Science Research Program Manager at NASA Headquarters, is 
included in Appendix C. Dr. Murphy was personally unable to 
attend the meeting due to other committments. 

The current document is the product of the contributions 
made by the SRAEC Principal Investigators, all of whom were 
present, at this annual workshop (Figure 1). This SRAEC annual 
status report is intended to 1) document responses to some of the 
questions addressed during the workshop concerning the progress, 
accomplishments and future directions of the project (Section 
II. A) and 2) to provide a concise summary of the justification, 
goals and objectives, and results of the individual SRAEC 
investigations (Section II. B), including an updated listing of 
the publications resulting from this research (Section II. C). 
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Figure 1. SIABC principal investigators at the Third Annual Work*" ■ p oa the 
ca«pus of MIT on January 30, 1985. The Pis are (front row, left o right) 
Drs. Badhwar , Saith, Diner, Qerstl, Fung, lahle; (second row) Ulaby , Strahler, 
Pearce, Paris, Coael, Moore, Kong; (third row) Irons, Nor»an, Fraser, Laag, 
¥a«derbilt; (foarth row) Miaes. 


II. SHAEC PROJECT RESEARCH RESULTS 
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II. A. PROGRESS AND ACCOMPLISHMENTS 


By the time of the Third Annual SRAEC Principal Investi- 
gator’s Meeting the SRAEC Project was beginning its fourth year 
and most of the investigators had either completed or were 
nearing completion of their third year of research. It therefore 
seemed an appropriate time to look back on what had been done, 
what progress had been made toward addressing the "research 
issues" that had been identified at the inception of the program, 
and what recommendation could be made now concerning the future 
directions of the program. 

A time for discussing the progress and accomplishments of 
the project was reserved on the agenda for the last day of the 
workshop to allow everyone to be updated on the progress of each 
of the investigations. Drs. Fawwaz Ulaby and Jim Smith accepted 
the responsibility for chairing these sessions and summarizing 
the results. 

Four questions were presented to the participants in the 
workshop; three of which dealt with the success of the project 
and one which considered a topic of considerable interest to most 
of the scientists concerning the potential establishing a 
centralized facility for artificial plant canopy experiments. 
After some general discussion, each of the participants was asked 
to individually prepare a written response to four questions. 
The four questions are given below, and the summarized, selected 
responses are given on the following pages of this section. 

In his cover letter (Feb. 20, 1985) reporting on the 
responses. Dr. Ulaby summarized by stating, "...the SRAEC program 
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is regarded by its members as ‘one of the most successful 


programs NASA has ever run’." The principal investigator 
responses that follow contain numerous suggestions and varied 
opinions regarding future directions. 


QUESTION NO. 1 

Excluding your own contributions, what are the three 
most significant technical advances or developments 
that were accomplished under the SRAEC Project? 

The responses to this question have been 
categorized according to: 

(a) General Contributions 

ft) Specific Contrib. — Optical Region 

fo) Specific Contrib. — Microwave Region 


QUESTION NO. 2 

What accomplishments have been made in the SRAEC 
Project relative to the goals defined in its inception 
document? Suggest changes in direction, approach, or 
administration . 

These responses have been grouped by: 

(a) Accomplishments Relative to Goals 

(b) Suggested Changes in Direction 


QUESTION NO. 3 

Are there "gaps" in the program? Should new areas be 
added? Should there be deletions? 


QUESTION NO. 4 

What are your views regarding the creation of a 
"National Canopy Facility"? 
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QUESTION #1 

Excluding your own contributions. Identify the three raost 
significant technical advances or developments that were accomplished 
under the SRAEC Program. 

The answers may be grouped into two classes: (a) general statemeiits 
that address the overall technical direction of the program, and (b) specific 
statements relating to either the Optical or the Microwave portions of the 
program. 

(a) General Contributions 

"Remote sensing represents, in essence, an inverse problem, llie 
characteristics of radiation fields are used to deduce surface properties. 

The two basic keys to the success of this are fundamental understanding of 
(1) the interactive properties of the target (the surface) and (2) the 
influence (if any) of the propagation medium. 

Efforts towards canopy modeling which have resulted in basic models 
for BRDF's and/or backscatter cross sections are therefore a very important 
component to our progress in attacking this inverse problem." 

"Create a nucleus of experts and expertise," 


"The program has helped to bring together investigators from both 
the optical and microwave areas. As a result, both group's knowledge has 
been strengthened by techniques used in the other area. Examples are: 
the optical people were using transport theory — now the microwave 
investigators are using it also. The microwave researchers were more 
aware of the coherence effects and the non-uniform angular radiation 
properties of leaves and surfaces. The optical researchers are 
measuring the BRDE of leaves and surfaces. In addition, the concentration 
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by the whole group on characterization of individual scatterers and 
the architecture of the canopy has led to an increased understanding 
of the relationships between measurements and biophysical parameters." 

"I and other optical modellers are moving toward better under- 
standing of the microwave problem. We all are working toward wavelength 
independent solutions to the scattering/absorbing problem." 


"Interdisciplinary interaction among various groups: optical and 
microwave, atmospheric and surface researchers; experimentalists and 
theoreticians with investigators who develop inversion algorithms." 


"Developing a better understanding of the interaction of radiant 
energy with surfaces, canopies, and atmosphere. Such knowledge enables 
development of better models and inversion algorithms." 


"The bringing together of multiple disciplines has been useful in 
getting the optical surface people to think about atmospheric problems 
and microwave contributions, and vice versa. This has led to greater 
sophistication in scene characterizations." 

"Smith's phase function work helps build an important bridge 
between optical, microwave and atmospheric work." 


"I do think that a very important result was getting these diverse 
pi's together and exposing them to each other's work. I think that the 
communication links that were formed and the sharing of common 
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problems concealed in the different terminologies of the subdisciplines 
was most helpful and will pay benefits for a long time to come." 


"To an outsider just joining this group, there appears to be too 
heavy an emphasis on a few areas (modeling of broadband electromagnetic 
radiation interaction with a few vegetative canopy types and measuring 
of some) to the exclusion of many other problems. I see little spectral 
work, no soils studies, no geology, no forestry, no mixing-model studies, 
no hydrology.... There are fundamental science problems in these areas, 
involving "SRAEC". Maybe nobody proposed to study them. Maybe announcement 
of opportunity needs to be more broadly written and distributed. 

The work that has been and is being done by this program appears 
to be excellent and important." 

"The development of a variety of computational tools for 
addressing the multidimensional terrain media problem, principally 
canopies . " 


"The beginnings of a solid, unified, attack on combined 
theoretical modeling of microwave/optical interactions." 


"Establishment of the communications channel among scientists 
in optics and in microwave field. A fruitful and stimulating exchange 
is being developed." 


"The potential unification of the theoretical frame work, which 
makes use of common mathematical approaches such as the radiative 
transfer theory." 
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(b) Specific Contrlbutlona — Optical Region 

"I am impressed with progress being made in the visible-IR area in 
inversion to determine actual relevant plant canopy properties from data. 
I found Strahler's success with a simple model particularly striking." 


"Established common theoretical basis for systematic evaluation 
of land remote sensing issues, i.e., radiative transfer, BRDF, phase 
functions, cross-sections, etc. Combined surface/atmosphere descriptions." 


"Identified new and important sensitivities of electromagnetic 
radiation (emitted or reflected from land surfaces) to intrinsic 
surface features, e.g., spatial heterogeneities and non-Lambert ian 
angular signatures." 


"Identified angular reflection features that seem invariant 
to atmospheric effects." 


"Tremendous progress has been made in testing, -evaluating and 
then improving canopy reflectance models in a systematic way. And with 
a clear confidence in these models, key biophysical parameters have been 
extracted. The 'Holy Grail' of remote sensing, the solution of the 
so-called 'signature extension' problem has been shown to exist and the 
role of canopy models in providing an understanding of this solution is 
clearly identified." 


"The complete atmospheric-canopy problem has been studied by 
transport methods. Investigators have made, initial attempts at obtaining 
boundary conditions for this model. Atmospheric correction methods of 
an elementary form has been proposed." 
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"The developments in atmospheric modeling coupled to the canopy 
reflectance models and measurements has ameliorated doubts? regarding 
capabilities to observe canopy BDRF's through the atmosphere. Therefore, 
satellite remote sensing remains a viable tool for botanical, land-surface 
climatology, and habitability types of research. 

"Sound practical approaches to atmospheric correction of digital 
image remote sensing data were discussed by Fraser and Dines (i.e., 
deblurring) . The applied remote sensing community has been in desperate 
need for such an approach. I believe this development is immediately 
applicable . " 


"The fact that the canopy reflectance models are demonstratbly 
invertible is most encouraging. The continuing development of these models 
is probably the largest accomplishment of the SRAEC program. 


"Tiie most significant technical advance is that the effect of 
specularly reflected light is now being included in canopy models." 

"The atmosphere/canopy is being treated as one problem rather than 
two. Phase functions pacific to each are being included." 

"I believe that all of the research being done in the space 

program is addressing the goals of the inception document. The 

approaches of these studies are appropriate. The individual researchers 

have developed a strong theoretical background which provides a sound 

basis for exploring extraction techniques of desired canopy characteristics. 

These 'dc.signs' qf canopy characteristics are varied for many similar 

problems inside and outside of NASA as we have discussed. I believe 
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this is the area which will give the SRA.EC program a great deal of 
publicity and future support." 


"Significant advances appear to have been made in both measurement 
and characterization of BDRF's of a number of complex surface types." 

"The potential merge of the vegetation canopy study and the 
atmospheric science, providing a complete framework for the study of 
the complete earth environment including earth terrain and atmosphere." 


"The recognition at optical wavelengths, that besdies ectral 
signature, there are other observables which distinguish di. .-rent crop 
canopies, i.e., angular reflectance signature, temporal greenness 
profile, or scene textural analysis. 

"The recognition that the information content of previous 
generation remote sensors (e.g., Landsat) is limited due to the 
restriction of nadir viewing. It appears that most interesting or 
useful effects appear at view angles greater than 45 degrees. A 
redirection of thinking to consider high off-nadir angles, taking into 
account the advantages as well as technical problems associated with 
such an approach, is a significant step. 

"Time dependency of canopy reflectances seems to be a reliable 
crop Identifier. 


"Jim Smith improved his canopy phase function." 
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"Angular patterns of reflectance is understood only qualitatively 
(bowl-shape, hot-spot, "invariant valley), but not quantiCetively . An 
artificial canopy with the possibility to set LAI, leaf, sizes, etc, 
may give us the hint to quantitative expressions." 

“Unscrambling pixels affected by the adjacency effect." 

“"Identifying distinctive features for remote sensing (hot spot, 
persistent valley, specular reflections in visible IR, branch structure 
and the like for millimeter waves. 

"Canopy reflectance models for inhomogeneous canopies." 


"More detailed models of canopy BRDF that consider leaf and 
soil 13RD1', partial cover and specular reflection." 


"Calculating phase functions for canopy so general RTE code can 
be used to obtain BRDF predictions. Need to do better work relating 
this phase function to actual canopy architecture." 

"The combined radiative transfer modeling, in a mu o more 
realistic way, of the atmosphere/surface problem with remote sensing 
implications . " 


"The potential merge of the vegetation canopy study and the 
atmospheric science, providing a complete framework for the study of 
the complete earth environment including earth terrain and atmosphere." 
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(c) Specific Contributions — Micro w ave Region 


"Improved theoretical descriptions that are beginning to take into 
account the actual geometries better." 


"Measurements of dielectric constants for materials and frequencies 
not well handled in the past." 


"By replacing stems with lossy dielectric rods and leaves with 
lossy dielectric discs, modeling techniques can be used to approximately 
reproduce observed experiments." 


"Rough surface effects have been incorporated into models 
showing Chat surfaces can exhibit strong backscattering at nadir angles." 

"Much better understanding of the dielectric properties of a 
single leaf through measurements." 

"Much improved understanding of the relative contribution of the 
component parts of a canopy." 

"Scattering models with much wider ranges of applicability have 
been developed." 


"Ulaby's measurements fill a great need." 


"Different contributions of different canopy elements to the returned 
and transmitted signals in the microwave region, assessment of importance of 
those elements for microwave remote sensing." 


18 



"Putting microwave sensing on a very firm mechanistic basis. 


especially in the work of Ulaby and Lang." 


"Trend towards understandable and usable models for backscatterlng 
in the microwave region . " 


"Some basic measurements like dielectric properties of soil, 
leaves, stalks, etc." 


"Microwave modelers considering leaf angle distributions and 
stems and branches." 
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qoEsrroM # 2 

Evaluate the a<i«o»olishjmetits of the SRAEC program relative 
to the goals defined in its inception document. Suggest changes 
in direction, approach, or administration. 

(a) Accomplishments Relative to Goals 

"As far as I know, little has been done under controlled conditions 
to solve the stress-response problem directly — yet this is one of the most 
important from an applications point of view. Ulaby has done a lot of 
work on seasonal variations, per se, as have the French and Dutch. 

"ulaby has made extensive measurements of soil dielectric as part 
of this and other programs. This is i\n accomplishment. 


"Goal 4 ('overriding necessity to relate to biophys parm.') 
is attacked by theory and experiment and may find culmination in work 
perceived for artificial canopy." 

"Perhaps the most important justification for NASA to support 
SRAEC is to develop a strong scientific basis for the remote sensing 
programs being developed for the 1990s. The entire radiation spectrum 
that has high transmission through the atmosphere will be utilized." 


'■"Relative to the accomplishments stated above, SRAEC program 
has made very significant progress towards the fundamental approach 
to the remote sensing science in unifying researchers from different 
disciplines, from different electromagnetic spectrum, from different 
schools of thought, and from both theoretical and experimental expertise." 


20 



"The first year« of the SRAEC program have led to advances 
In developing the tools needed to address the question outlined in 
the inception document. While some of the questions may have not 
yet been answered, we are in a significantly better position to address 
them now than we were three years ago. The flexibility of the program 
is a major plus, allowing researchers free freedom to pursue new ideas 
and not be rigidly held to their proposals." 


"I think the SRAEC research program has pretty well followed the 
goals objectives set out for the program at its inception. . The original 
document called for a very broad research program with many topics 
and needs listed. Although not all of these areas have been attacked, 
the main thrust has clearly been along the lines of the plan. It seems 
to be funding a good balance of theoretical work, measurement, and work 
that is a bit more application-oriented. Except for the usual NASA 
bureaucratic problems (funding delays, etc.), the program has been 
well administered. I don't think it really needs much change." 

"Modelers should not start impacting the applied programs. 

A question like: 'what is the best CR model for a given vegetation?' 

1 believe, should be addressed by the group." 


"Program is proceeding toward stated goals at a rapid pace and 
I would not change the approach or administration and would let the 
direction continue to be defined by the original document." 
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"Much better than I would have thought. While individuals may 
come and go, I believe it is important to keep some collective memory. 
I believe administration by NASA has been exemplary. Overall, a first 
class group doing flrs^": class work under enlightened NASA enviornment 
fostering creativity and comraderie." 


"Relative to the accomplishments stated above, SRAEC program 
has made very significant progress towards the fundamental approach 
to the remote sensing science in unifying researchers from different 
disciplines, from different electromagnetic spectrum, from different 
schools of thought, and from both theoretical and experimental 
experties. " 
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(b) Suggested Changes In Direct.ton 


"Perhaps the more glaring gaps in the program are filled with 
the induction of geologists into the group. We need better coverage 
on the properties of rocks and soils. We need thermal band coverage. 
These are addressed to some extent with the recent program addition. 
On the other hand, we have not, as a group, looked at snow or ice. 

We have not looked at the sea although useful observations might be 
possible, especially for coastal waters. 


"Little has been done to improve cross-polarized measurements. Use 
of circular reflectors of parabolic shape leads to significant cross 
polarization. Low-cross-polarization arrays or horns should be used 
in such measuremcMits . Cross-polarized returns have significant potenti-''s, 
as demonstrated years ago with tlie Westinghouse 35 GHz radar, and more 
recently with ERIM, CCl^S .ind .IPl, railars." 


"Nothing has been done in this program with regard to snow and 
ice dielectrics, although some other work relates to them. The problem 
is complex, and perhaps should be left to the Arctic and Air Force 
programs ." 


'I see no need for a change in direction." 


"I agree with Dr. Smith's statement that the SRAEC program needs 
to extend applications of radiative transfer theory to other terrestrial 
media. I am biased towards the study of soils. Other media obviously 
include snow, ice, exposed rock, and even water bodies (e.g., the 
turbid waters of the Great Lakes)." 



"Steps need to be taken to ensure the transfer of SRAEC accompli.slinu'nts 


to applied remote sensors. In particular, these accomplishments are of 
immediate and great utility to those scientists planning retrospective 
and field observations for ISLSP." 


"The greatest strength of the program is the heavy emphasis on 
dev''lopment of better understanding of the canopy radiation field using 
modeling as an approach." 


"As implemented, this emphasis has elevated determination of 
LAI, leaf angle distribution etc. (e.g., a very limited set of 
parameters) as an end objective. This is both good (we need to know 
them) but also bad (we must not lose sight of the goal to obtain 
informat i on . The variables or parameters mean nothing by themselves. 
To be useful, they mu.st be put in some context and they must be used 
with other information. The modeling efforts tend to neglect these 
aspects . 

Thus, there is a need to better relate the modeling activities 
to basic botanical information." 


"SRAEC investigations should commence on all subjects that are not 
now being done properly. Examples may be snow-ice, hydrology, wet 
lands, soil and shore erosion, etc. 


"There were almost no intercomparisons between different models 
addressed to the same topic." 
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"Validation, verification and sensitivity analysis of existing 
models are still needed. This can be done by using numerica simulation 
in the case of rough soil surface and perform basic experiments on 
artificial canopy of known potentials in the case of vegetation layers." 

"Extension of the existing models to heterogeneous media is 
being done but not yet completed." 


"To relate electromagnetic parameters to the classical 
biophysical descriptions — much progress has been made over recent 
years both experimentally and from mcdelling viewpoint." 

"Stronger links between experimental programs and theoretical 
work are still desirable. Should encourage joint efforts or proposals 
from theoreticians and experimentalists." 


"Basic measurement should now be de-emphasized . Some modeling 
effort be directed to check models against observations already made." 

"SRAEC program needs more research related to incoporating as 
many aspects of the soil-plant-atmosphere system as possible. We have 
emphasized radiation exchange, but numerous other processes are closely 
rela, ;d to radiation exchange such as evapotranspiration and dry matter 
production and there should be some attempt to bridge this gap with as 
rigorous an approach as possible," 
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"Further developments are needed in relating classic biophysical 


descriptors to electromagnetic parameters (e.g.. Biomass n6eds to be 
estimated as opposed to LAI and LAD). These needed developments, 
however, will likely follow from the current research direction." 

"One gap is in not using the theoretical radiation caoability 
assembled in SRAEC program to address problems of leaf BRDF as it 
depends on leaf structure - internal (cells, chloryplasts, air gaps, 
water, etc.) and external (epidermis, cuticle, surface irregularities, 
etc . ) . " 


" Much more emphasis needs to be placed on the inverse problem 
both in terms of rigorous theoretical models and simple algorithms that 
can be used efficiently on millions of pixels. This universe problem 
needs to include the atmosphen in the near future." 

"1 think there were not enough model intercomparisons." 
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QUESTION # 3 

Are there gaps in the program? Should new areas be added? 
Should there be deletions? 


"It seems to me that a gap exists in tying the microwave and vis-IR 
response together. Tlie ulitraate application is going to involve sensor 
combinations with combined use of vis-IR and microwave when vis-IR 
is possible, and use of microwave alone when clouds obscure the ground." 


"If natural vegetation is the goal these days, it seems to me that 
we are not doing justice to some forms of natural vegetation. Specifically, 
I see almost no effort on grassland and sparse desert \ getation. Although 
other programs are addressing snow cover, sea ice, etc., at least some 
link to these would be appropriate. Moreover, urban areas are a part of 
tl)e environment — a very big part in some places — yet nothing seems to be 
underway with regard to them except a tiny effort on, our part to look 
at radar rources in wooden buildings." 


"Another area that seems to be receiving too little emphasis is 
uncovered lithology. It is fine that it is starting in the thermal IR, 
but the visible and microwave regions are not being treated. We have 
done some work in the past that shows that one really can get lithologic 
signatures in the microwave, but the surface has barely been touched." 


"Much work needs to be done for grasses and forest canopies in 
both canopy models and hope to extract key biophysical variables. 

A clear need to attack (or combine) canopy reflectance models with 
atmospheric models needs to be done, though a start has been made 
(f'erst.nl ) ." 
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"The important gap in the program is the lack of a systematic 
plan- to verify individual models with carefully controlled- experiments . 

The experiments should be made on both real and artificial canopies. 

All investigators should be notified of the experiment in advance and 
have the opportunity to have ground truth that might need to be measured." 


"In the optical area, there is a need to have the individual 
scatterers in the canopy carefully modeled." 


"Both in optics and microwaves, the inversion of present models 
should be carried forth. Since this is our ultimate goal, the groups 
should start to address this portion of the work mo-re seriously." 


"There are gaps in two areas. There is a crying need for 
research to provide means to obtain the input parameters needed for use 
in canopy radiation models. That is, we need better instruments to 
measure/cstimate LAI, LAD, etc in a timely manner. At present we are 
quite capable of collecting enormous amounts of spectral reflectance data 
on numerous fields but only limited amounts of supporting data on 
only a few fields. Thus we need a machine to measure the canopy 
geometry. We need support for collection of a limited amount cf scanner 
polarization data, preferably polarized AIS type data, in support of 
my area and other modellers . Without these data we will not be able 
to address these questions. 

(1) How is light, scattered and polarized by a plan canopy, transmitted 
through a disturbing atmosphere over a large area. (2) Is polarization 
data another independent source of information about the scene — 
providing another 'channel' of information for discriminating species 
and assessing the condition of the plant canopy on the ground." 
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"An error analysis needs to be made for the entire remote sensing 
system to determine the amount of research that is required for each 
component of the system. The system consists of the surface, canopy, 
atmosphere, sensor, processing the measurements, and finally inversions. , 
If the accuracy of only one component of the system is only five percent , 
say, then it may be unnecessary to achieve much greater accuracy for 
any other component. Sensors should not be built at great expense 
with greater accuracy than is justified by the accuracy of the remaining 
components in the system." 


"Lack of interaction between SRAEC investigators and meteorologists 
and climatologists." 


"Only uhe thermal region is missing in the existent programs." 


"Ihe gaps in the program are related to tlie issues addressed in 
question 1, which are in fact the future directions that we should be 
concentrating on. It v/ill be of fundamental importance to unify the 
optics and the microwave approaches, to include earth terrain media 
other than vegetation canopy, and extend the atmospheric study to 
include both the lower and the upper atmosphere. With such 
interdsiciplinary approaches, we should be able to establish credential 
in accountability and applicability to NASA goals such as the global 
habitability studies, and at the same time remain in the domain 
of fundamental research." 
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"Gaps: The thermal is an obvious gap that has been mentioned. 

Snow surfaces are not being addressed, and should. New work in soils 
and rock spectra are extending our knowledge of scene radiation to 
include phenomena beyond plant canopies. I think the program should 
continue to broaden in this way, but still keep its emphasis on 
fundamental research unknowns. 


"There is now a critical need for more coherent involvement 
of measurement and modeling. The so-called "National Canopy Facility" 
is all manifestation." 

"Inversion of models to obtain biophysical parameters from 
remotely sensed data." 


"Considering tlic funds available and the number of people involved, 
the program is very balanced. The biggest gap I see is a strong research 
program on modeling thermal IR phenomena." 


"Thermal IR studies of the surface, including multi-angular 
studies as well as investigation of IR sounding to characterize the 
atmosphere . " 


"The addition of someone in the lidar field would be useful, to 
get the optical passive and active communities talking to each other, 
as well as to give the laser community a better sense of the type of 
atmospheric information that is needed in order to correct images for 
atmospheric effects." 
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"The gaps in the program are related to the issues addressed in 
Question 1, which are in fact the future directions that we should be 
concentrated on. It will be of fundamental importance to unify the 
optics and the microwave approaches, to include earth terrain media 
other than vegetation canopy, and extend the atmospheric study to include 
both the lower and the upper atmosphere. With such interdisciplinary 
approaches, we should be able to establish credential in accountability 
and applicability to NASA goals such as the global inhability studies, 
and at the same time remain in the domain of fundamental research." 
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QUESTION H 4 

What are your views regarding the creation of a National Canopy 
Facility? 


"Considering the effort requirevi to establish a chosen and 
various LAI and LAD for a canopy of credible size (number of "leaves", 
etc), it would seem to me that it might be better to devote effort to 
devising rapid, comprehensive, and accurate techniques for measurement 
of the geometrical leaf descriptions than to constructing a hardware 
model of large scale. 


"The use of models in experiments is important in many investi- 
gations, but it is also important that individual experimenters use simple 
models at their own sites to answer simple questions quickly. The whole 
job cannot be done effectively with one or two massive models used by 
everyone. Perhaps these are desirable for comparison purposes and for 
classes of experiments where simple models will not do." 


"Geometric scaling should be possible for answering most modeling 
problems. This means that the models can be much smaller than the 
real world, which makes them easier to use in some ways. For microwaves 
this means that a higher frequency must be used than the frequency that 
would be used with the real thing, but this is not hard to do. I see no 
real problems with doing this scaling in the optical-IR regime, and it 
should allow the models to be much more tractable. Moreover, for a dense 
canopy, there is no need to model the entire thickness, since the 
radiation will not penetrate beyond a certain depth anyhow!" 
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"I recomniend a vary inexpensive, simple canopy at first to see 
if tnis idea is really worth putting the time and money into." 

"Absolutely necessary because (1) angular distribution (besides 
time variation) of canopy reflectance seems to be a valuable crop 
identifier and (2) angular pattern of reflectance is understood only 
qualitatively (bowl-shape, hot-spot, invariant valley), but not 
quantitatively. An artificial canopy with the possibility to set 
LAI, LAD leas sizes, etc may give us the hint to quantitative 
expressions . " 


"I believe there is n strong need to establish a "controlled" 
facility for a variety of representative media (not just canopies) . 

Some tests would be very simple/some complex. Considerable thought 
should go into such a program-wide resource — it is important to do it 
right and assure that sufficient resources be allocated for a sustained 
effort — not a "one-shot" deal. I am reasonably confident that 
funds from other agencies could be identified to matdh NASA funds. 

I am quite concerned that such a facility be established and 
maintained by a group/groups that already have existing measurement 
expertise in microwave/optical regions. I think a combined facility 
could be beneficial. I need to think through some steps/ways to get 
things going. I have some significant requirements for the op.. teal case. 


33 



"I believe there is well-defined value, over all programmatic 
interests, in the NCF: models are not verified in full as being 
invertible for enough distinctive canopy features, we need to 
understand "hidden variables" or hidden assumptions about canopy 
structure (especially pair cotrelations of leaves or other scattering 
structures) . 

"I think that construction of a National Canopy Facility will be 
a total waste of scarce resources. For homogeneous and uniform canopies, 
CR models exist which represent the CR quite accurately (within the 
measurement accuracy). For complex canopies, the models are not as 
good but a complex artificial canopy will be hard to make and will 
itself introiluce unwanted nrtifacts. 

I support the concept of nn individual P.I. proposim; in 
artificial canopy experiment as part of his/her proposal." 

"I favor the idea of an artificial canopy being explon'd. However, 
there must be at least a promise of learning something from this that 
we do not already know. Also, I would hope that it is not so conceived 
and constructed that it becomes an end in itselt without contributing 
to the real vegetation remote sensing problem. 

I think that a significant effort should be funded to explore 
designs for such a facility using expertise from a broad cross section 
of the biological and remote sensing communities. This would involve 
materials, questions, logistic questions, radar vs. optical, and outlines 
of some reasonable projects that it would be used for. 
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After this one year or more of study, the question should be 
posed well enough to be evaluated in a peer review system. Whether 
this artificial canopy has naything to offer will depend strongly on 


the design and we need to know specifics that are generated fuom a broader 
cross section of potential users." 


"I don't think large amounts of money should be spent on making 
an artificial canopy closely resembling plants when one coul ^ spend 
69c on a seed packet and grow the real thing! To me, the utility of 
an artificial canopy would be to convincingly demonstrate that canopy 
properties can be recovered uniquely for relatively simple (but non-trival) 
geometries, because if that can't be done, then I don't see how crop 
situations can be inverted. Thus, there is ho pressing need for the 
artifical canopy to resembly plants other than in some gross way. If 
such a simple canopy can't be constructed for under, say, $10,000, we're 
getting into dollar amounts that might be better spent elsewhere, 
especially in this time of budget limitations." 
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II. B. PRINCIPAL INVESTIGATOR SUMMARY REPORTS 


Each of the principal investigators was requested (see 
Appendix A) to bring to the Third Annual SRAEC Principal Investi- 
gator’s Meeting a succinct summary research report that provided 
information on the 1) justification, goals and objectives; 2) 
experimental approach; and 3) study results, with up to three key 
figures to illustrate their findings. In addition, a separate 
listing of publications and presentations on their research was 
also requested (see Section II. C.) 

The remainder of this section contains the individual 
investigator’s summary research reports in the following order: 


Smith Reflectance Modeling (p.37) 

Norman Bidirectional Reflectance Modeling of 

Non-Homogeneous Plant Canopies 
(p.45) 

Kimes Understanding the Radiant Scattering 

Behavior of Vegetated Scenes (p.53) 

Vanderbilt Measurement and Modeling of the Optical 

Scattering Properties of Crop 
Canopies (p.59) 

Strahler Geometric-Optical Modeling of a Conifer 

Forest Canopy (p.66) 

Badhwar Crop Characteristic Research: Growth 

and Reflectance Analysis (p.71) 

Conel Ground Reflectance Measurements from 

Thematic Mapper Data (p.78) 

Kahle Thermal Infrared Geometric Remote 

Sensing Research (Co-I, 

Bartholomew) (p.85) 

Diner Atmospheric and Scene Complexity Effects 

on Surface Bidirectional 
Reflectance (p.92) 
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Fraser 


Pearce 

Gerst 1 

Fung 

Kong 

Lang 

Ulaby 

Moore 

Paris 


Atmospheric Effect on Remote Sensing of 
the Earth’s Surface (p.99) 

The Characterization of Surface Varie- 
gation Effects on Remote Sensing 

(p. 106) 

Multidimensional Modeling of Atmospheric 
Effects and- Surface Heterogeneities 
on Remote Sensing (p.113) 

Scattering Models and Basic Experiments 
in the Microwave Regime {Co-I, 
Blanchard) (p.ll9) 

Remote Sensing of Earth Terrain (p,125) 

Discrete Random Media Techniques for 
Microwave Modeling of Vegetated 
Terrain ( p ^ 129 ) 

Microwave Dielectric and Propagation 
Properties of Vegetation Canopies 
(p. 135) 

Determination of Backscattering Sources 
in Various Targets (p.l4i) 

Active Micowave Properties of Vegetation 
Canopies (p.l47) 
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gSFLECTftMCE MODELING 
DC. James A. Smith 


N86- 13795 


College of forestry & Natural Resources 
Colorado State University 
Fort Collins, Colorado 80523 


PROJECT OBJECTIVES 

The overall goal of our work has been to develop a set of conputa- 
tional tools and media abstractions for the terrain bidirectional reflec- 
tance problan. The modeling of soil and vegetation surfaces has been 
anphasized with a gradual increase in the conplexity of the media 
geometries treated. Pragmatic problems involved in the combined modeling 
of soil, vegetation, and atmospheric effects have been of interest and one 
of our objectives has been to describe the canopy reflectance problem in a 
classical radiative transfer sense permitting easier inclusion of our work 
by other workers in the radiative transfer field. 

APPROACH 

Th’-i ^plication of radiative transfer theory to soil and vegetation 
media requires an abstraction of the scattering and absorbing properties of 
the medium in terms of biogeophysical attributes, a specification of the 
geonetrical boundary of the medium, and a mathematical solution technique 
to the radiative transfer equations. We have abstracted our media to con- 
sist of a collection of discrete scatterers, e.g. leaf or soil facets, with 
individual diffuse and specular optical properties. Bulk properties of the 
media are usually derived frcm suitable averaging over the statistical dis- 
tribution of the individual elenents within the medium. For sinple plane 
parallel nedia geometries, we are using classical solutions to the radia- 
tive transfer equations, e.g. the discrete ordinates/ finite difference 
method of Chandrasekhar, the Adding method of Van de Hulst, or the two 
stream method of Kaufman. For complex geanetries with irregular rough sur- 
faces, we are employing Monte Carlo techniques. 

RESEARCH FINDINC^ 

We will briefly summarize three areas which have been addressed since 
the last SRAEC annual meeting. These are orthogonal ecpansion of our phase 
functions in terras of total scattering angle and exanple use with standard 
radiative transfer code, e.g. Kaufiran's two stream approximation; our Monte 
Carlo reflectance modeling for complex soil surfaces which now include 
Fresnel facets; and our current work in modeling ccsnplex vegetation sur- 
faces. 


flrthogon^ -ExEansion ^ .UfeiliRatlon of CanoEy Phase Functions 

We have previously described our calculational technique for averaging 
over leaf angle distributions and leaf optical properties to estimate the 
canopy scattering diagram as a function of the global incident and exitant 
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angles. ^The jiiase functions we presented were given in terms of incident 
angles, and exitant angles, As Ross has indicated, it appears 
that canopy jdiase functions, unlike some of the more symmetrical atmos- 
pheric cases, may not necessarily be explicitly formulated in terms of oily 
a total scattering angle, 6. We have i^nt seme time pondering this cir- 
cumstance. Recently, we have adapted an approach outlined by Liou for ice 
crystals v^er^y "average" phase functions can be computed. The procedure 
essentially averages over ranges of ju,eJ and ja' which yield the same 
total scattering angle, e, i.e. 

cos 0 - /l' + (1“A1)^^2 (1-/1')^^^ cos {(Sr- >6') (1) 

We term the resulting scattering diagrams, "stylized" jhase functions, 
we do not “believe the results are strictly correct in the fine details, 
because of both fundamental reasons and numerical problems. Hov/ever, as 
Van de Hulst has indicated, for many calculations the fine details may not 
be all that important. Theoretically, the fhase functions are much more 
trachable and be expanded, in Legendre polyncmials, for example, as fol- 
lows; 


P(cos 9) = Pj^(cos e ) + P 2 (cos 0 ) + ... (2) 

and, thence, into a variety of useful forms. The table below shows an 
example calculation for a spherical leaf angle distribution for visible and 
near infrared wavelengths. The jhase function, and, henc«, expansion coef- 
ficients vary with the vegetation optical properties. 



^ = 0.55 Ai m 

^ = 0.65 ju 1 


1.009 

1.018 


0.057 

-0.238 


0.174 

0.180 


-0.022 

-0.053 


0.021 

0.028 


-0.001 

-0.021 


0.021 

0.022 

NOTE; 

u^ should be 1.0 
calculations is 

; sane measure ( 
indicated by thi 


^ = 0.80 M m 

1.018 

-0.271 

0.183 

-0.051 

0.033 

-0.024 

0.026 

: the numerical precision of the 
deviations from this value above. 


Such expansions can be easily incorporated into standard radiative 
transfer code. We have used the above values in Kaufman's published modi- 
fied two stream itethod — the program code was chtained from the Atmospheric 
Science Department at Colorado State University. Figure 1 is a plot of the 
(spline - smoothed) stylized phase funchions at the three wavelengths 
above. 
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iitontg JSqU Reflectance Modeling 

At last year's SRATC some preliminary results from a Monte Carlo model 
for the reflectance of a rough soil were presented. The thrust of our work 
in this area over the past year has been to extend the usefulness of this 
model to surfaces having non-Lambertian elanental reflectance properties. 
In particular, we were concerned with the effect of a distribution of 
Fresnel facets superiirposed on macroscopic surface irregularities . The 
reflectance properties of an elonent of area of surface covered by these 
facets was determined theoretically, and then this reflectance function was 
substituted for the original Lambertian one in the Monte Carlo model. 

In the substitution of reflectance functions into the Monte Carlo 
model it was desirable to allow the nvost general forms possible. For this 
reason two alternatives were investigated: the approximation of the refle<?” 
tance functions in terms of sfiierical harmonics, and the use of a lookup 
table generated from the value of the reflectance functions at discrete 
angles. After some preliminary testing it appears that the former approach 
is slightly better when the reflectance function in question is sinple, 
however for more conplex functions, and in particular for the functions 
arising frcan distributions of Fresnel facets, the lookup table approach is 
the only reasonable one. 

Thus far, we have only considered the case when the Fresnel facets 
have a uniform distribution of inclination angles. An example of the sort 
of results detained for a surface with a siirple cosine undulation is given 
in Fig. 2a. In this case light was incident on tiie surface at a zenith 
angle of 45 degrees and an azimuth perpendicular to the direction of the 
rows. The results are given by plotting contours of equal reflectance fac- 
tor on a polar plot in which distance from the center indicates the zenith 
view angle, and the azimuth of the view corresponds directly to that on the 
plot. Figure 2b shows the bidirectional reflectance distribution function 
for the same surface possessing purely Lambertian elemental scattering pro- 
perties. 


Refle c tan ce From .Ccpplex Ve getation Surfaces 

After reviewing the bidirectional reflectance distribution modeling 
work for vegetation media currently underway in the SRAEC fundamental 
research effort, we concluded that it would be useful to attack the problem 
of reflectance from a rough surface vegetation layer from a multiple 
scattering, radiative transfer perspective. The approach we chose is an 
outgrowth of Norman's work and conplements Strahler's geccnetric optics 
approach utilizing stylized structures for tree crowns. Basically, we are 
adapting a Monte Carlo code by utilizing probability of gap expressions 
which take into account the varying optical path length through vegetation 
layers as a function of canopy density, macro-structure and look angle. 
Geometrical routines are used to compute required parameters from a canopy 
height profile. In practice we anploy Fourier series expansion of canopy 
height profiles. Later we may examine statistical structures inherent in 
the data via power spectra analysis and formulate appropriate models. 
Ttiese ideas are partially generated from some thermal modeling work we are 
performing over a forested site at Oak Ridge National Laboratory and we are 
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capitali 2 ing on the canopy structure data available for this site. 

This work is in preliminary stages, but briefly our approach follows: 

1. An adapted version of our Monte Carlo vegetation reflectance model 
will be used to handle the multiple scattering calculations. 

2. Probability of gap expressions are computed for photon interactions by 
considering both micro (i.e. leaf angle distributions) and macro 
(canopy height profile) structure, i^cifically, 

Po(©) = exp - LAI*g(6)*D(e) /h (3) 

where g( 6 ) is the mean canopy projection in direction 0, LAI is leaf 
area index, and D( 6 ) is the gecnietric path length fr<xn the bottom of 
a canopy later to the top profile in direction 0, and h is the canopy 
layer thickness. 

3. Appropriate geometrical confiutational routines have been developed to 
calculate D(0) for an arbitrary canopy height profile. We have 
applied these routines to saitple data transects frcm the Oak Ridge 
Site. In practice, the canopy height profile data are first 
"detrended" using least squares regression, and the the power spectra 
computed for the detrended data in order to develop Fourier series 
representations for the calculation routines. Figure 3 shews a sairple 
detrended height profile, and schanatically illustrates the varying 
path lengths through this canopy layer as a function of look angle. 
Multiple voids in the canopy are correctly handled. 

SIGNIFICANCE AM2 JSSSE MAJOR STEP 

The main significance in our work is an increased understanding, 
expressed via a variety of computational procedures, of terrain bidirec- 
tional reflectance properties. To our knowledge, there were no macro soil 
reflectance models for complex surfaces incorporating the conplete multiple 
scattering solution and accounting for shadowing effects. Our calculation, 
treatment, and utilization of canopy phase functions also appears to be 
somewhat novel, although we have built upon the efforts of earlier workers. 
Extensions could be made to our work to include, for example, the incor- 
poration of media subelement models to account for physiological or polari- 
zation properties. 

In the eight months ranaining of our third year of funding the final 
step in our research focus will be to complete a tractable treatment of 
complex vegetation surfaces amendable to our Monte Carlo canopy reflectance 
modeling techniques. In addition, we will more gracefully complete several 
loose ends in our completed work including the generation of orthogonal 
expansions of our canopy phase functions and the utilization of the models 
developed under this stu<%^ to illustrate soil/vegetation/atmospheric 
effects in such applications as AVHRR-derived vegetation index calcula- 
tions. 
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Figure 2b. Polar plot of bidirectional reflectance distribution 
function from a rov; soil surface with strictly 
Lambertian elemental reflectance properties. 
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BIDIRECTIONAL REFLECTANCE MODELING OF NON-HCMOGENEOUS PLANT CANOPIES 

John M. Norman 
University of Nebraska 
Lincoln, NE 685fi 

Objective 

The objective of this research is to develop a 3-dimensional radiative 
tran.-fer model for predicting the bidirectional reflectance distribution 
function (BRDF) for heterogeneous vegetation canopies. The model (named 
BIGAR) considers the angular distribution of leaves, leaf area index, the 
location and size of individual subcanopies such as widely spaced rows or 
trees, spectral and directional properties of leaves, multiple scattering, 
solar position and sky condition, and characteristics of the soil. The model 
relates canopy biophysical attributes to down-looking radiation measurements 
for nadir and off -nadir viewing angles. Therefore inversion of this model, 
vdiich is difficult but practical, should provide surface biophysical 
properties from radiation measurements for nearly any kind of vegetation 
pattern; a fundamental goal of remote sensing. Such a model also will help to 
evaluate atmospheric limitations to satellite remote sensing by providing a 
good surface boundary condition for many different kinds of canopies. Further 
this model can relate estimates of nadir reflectance, which is approximated by 
most satellites, to hemispherir?’ reflectance, which is necessary in the 
energy budget of vegetated surfaces. 

Technical Approach 

The approach to this research requires development of the mathematical 
equations and computer coding of the heterogeneous-canopy model, better 
characterization of leaf and soil properties which are a serious limitation at 
the present time, and finally comparison of model predictions with field 
measurements that are obtained from other investigators. The Bidirectional 
General Array Model (BIGAR.) is based on the concept that heterogeneous 
canopies can be described by a combination of many subcanOpies, wliich contain 
all the foliage, and these subcanopy envelopes can be characterized by 
ellipsoids of various sizes and shapes. Tne foliage within an ellipsoid may 
be randomly positioned or non-randonly positioned. Estimates of multiple 
scattering are obtained by transforming each point of interest in a particular 
ellipsoidal canopy to an equivalent one-dimensional canopy and solving the 
radiative transfer equations for this simple case. The BRDF for individual 
leaves has been measured so that appropriate leaf properties can be used in 
the model. The BRDF for several soils has been measured and modeled with a 
simple block-modeling approach to provide a reasonable lower boundary 
condition for the canopy model. Field measurements ot com and soybean canopy 
BRDFs are compared with odel predictions. 
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Research Results 


The results of this study are three fold; 1) A simple soil BRDF model 
tested with measurements, 2) laboratory measurements of leaf BRDFs for live 
corn and soybean plants grown in the greenhouse and the field, and 3) the 
development and validation of a 3-dimensional canopy BRDF model of 
heterogeneous vegetation that combines soil, leaf, and canopy-architectural 
information. 

The simple soil BRDF model approximates a soil aggregate by a single 
rectangular block placed on a fixed lot area vidth a fixed orientation relative 
to the sun azimuth. The relative reflectance factor in any view direction is 
computed from projections of block faces and the shadow of the block on the 
horizontal. Fig. 1 contains a comparison of model results with measurements 
in the solar waveband over a rough surface that was recently plowed from sod 
with furrows approximately 15 cm deep. The RMS difference between the 
surfaces is 2.2% in reflectance. The simple-black-model predictions fit 
measurements from smooth (0.5-1 cm gravel) and medium rough (2 to 5 cm clods 
from multiple tillages) even better with 1% RMS differences in reflectance. 

Leaf BRDF measurements were made at 3 source incidence angles (20, A5, 

70) and about 50 view angles for intact com and soybean leaves attached to 
their respective plants. The distribution visible and near-infrared, 
reflectance and transmittance factors is shown in Fig. 2 for a 45 degree 
source Inciucucc angle. The leaf hemispherical reflectance, which is 
essential to all models of vegetation radiation exchange, was calculated from 
the integral of the BRDF (Table 1). Clearly normal incidence valves which 
currently are being used in other models may not be appropriate since mean 
leaf-sun angles for most canopies are between 45 and 60 degrees and not near 
normal incidence. 

The 3-dimensional bidirectional reflectance general array model (BIGAR) 
predictions have been compared to measurements on com and soybeans that were 
obtained from the Laboratory for Applications to Remote Sensing (LARS) , 

Purdue, University, West Lafayette, Indiana as part of the AGRISTARS program. 
An example of how ellipsoids are used to represent a young com canopy is 
shown in Fig 3 and a comparison of measured and modeled results is in Fig. 4. 
In general the model predictions agree very well with measurements considering 
that soil BRDF measurements cannot be made for the same site as the vegetation 
BRDF measurements. 

A simple 3-term empirical equation has been derived to fit measured and 
modeled, soil and canopy BRDFs over a wide range of wavelengths. The RMS 
reflectance difference between the 3-terra empirical equation and modeled or 
measured distributions (including view zenith angles from 0 to 60 degrees) is 
about 0.2% in the visible and 2% in the near infrared for soybeans and several 
soils; a remarkably good fit. 
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significance of Results 


Significant research results from this study impact on vegetation remote 
sensing in three ways; 1) relation of soil BRDF to roughness and sun angle, 2) 
interactions between soil BRDF and canopy reflectance properties that lead to 
confounding the discrimination of sparse vegetation using angle-of-view 
observations, and 3) interpretation of leaf BRDF measurements to provide more 
suitable leaf spectral properties for models of vegetation radiative transfer. 

The very simple block model, which has beem verified with field 
measurements, relates the soil BRDF to roughness. Unfortunately the roughness 
is a relative roughness and thus not the same kind of "roughness" used by soil 
physicists, Tlius a disked field and a gravel lot can have similar relative 
roughnesses (ratio of size to spacing) and BRDFs but represent very different 
roughnesses to the soil physicist. 

The BRDF of sparse vegetation is not as different from soil as is 
desirable for discrimination based on view angle observations. However the 
structure of sparse canopies may be distinguishable from each other or soil 
backgrounds with zenith view angels of 60 . Unfortunately this is not 
possible from satellites. In addition we now understand how a BRDF is 
different for a soil below a canopy than a soil in the open even though their 
physical characteristics are identical. The 3-term empirical equation fits 
soil and vegetation BRDFs so well that it may greatly simplify extraction of 
canopy properties from off-nadir observations for natural surfaces so that use 
of such dciLci oiay become practical. Perhaps it may even reduce the number o£ 
off-nadir angles needed for obtaining canopy properties from satellite 
observations. 

During this study we have measured leaf BRDFs for com and soybean; a 
rather formidable task that few have attempted. From this work we have 
learned that models based on normal -^incidence integrating sphere leaf spectral 
properties, and to date nearly all leaf spectral property data is for normal 
incidence, are using underestimates of leaf reflectance and overestimates of 
transmittance. Using the more appropriate leaf spectral properties may change 
canopy reflectance by 5 to 15% of the reflectance value. The improved 
estimates of leaf spectral properties will also improve estimates of canopy 
water use and photosynthesis. In fact the reduced water use of an isogenic 
line of dense pubescence (dense hairs) soybeans could only be explained by the 
leaf BRDF and not normal incidence integrating sphere measurements. 

Future Work 


Future research should emphasize the inversion of 3-dimensional 
vegetation- soil bidirectional reflectance factor models and develop simpler 
inversion algorithms from exercising these complex models. Both angular (view 
direction) and wavelength information should be combined in the inversion 
process to extract the maximum aixcunt of information; perhaps inverting the 
distribution of normalized-difference with view angle may hold promise. The 
logical extension of this 3-D model is to include more complex canopoy 
architecture such as conifer trees, and also extend the wavelength band Co the 
thermal . 
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Table 1. Hemispherical reflectance and transmittance of com and soybean 
calculated by integrating the results of leaf BRDF measurements. 

Incidence Visible KIR 


Crop 

Angle 

Refl (%) 

Trans (%) 

Refl (%) 

Trans (%) 

Com 

0 

7.7 

3.9 

- 

- 


20 

9.1 

4.2 

34.6 

39.3 


45 

9.3 

3.4 

33.0 

.4 


70 

14.2 

3.9 

45.2 

45.0 

Soybean 

0 

8.3 

4.9 

- 

- 


20 

9.3 

4.3 

41.7 

30.9 


45 

9.6 

3.7 

44.9 

33.2 


70 

15.3 

3.1 

51.6 

32.0 
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Fig. 1. Soil BRDF predicted from simple model(H=0.4, L=0.85 
W=0.1, =0.67) and measured for a rough plowed soil 

with 15 cm furrov7s anc clods for a solar zenith angle 
of (a) 63° and (b) 28°. The RMS difference between 
measured and modeled for (a) is 2.2% and for (b) is 1.1% 
in reflectance units. Maximum zenith view angle is 60° 
(outer rim), center is nadir and sun azimuth is 0°. 
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Fig. 2. Leaf BRDF and bidirectional transmittance distribution function 

at a source incidence angle of 45° and azimugh of 0° for (a) corn 
at near infrared vavelengths, (b) corn at visible wavelengths, 

(c) soybean at near infrared wavelengths, and (d) soybean at 
visible wavelengths. 
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Fig. 3. Sketch of a corn canopy vitb LAI^ 0.4 including the 
ellipsoids used to approximate the canopy envelope. 
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(b) 


Fig. 4. Comparison of measurements and model predictions for a 
corn canopy of = 0.4 having north-south rows. (a) 
solar azimuth «• 7°, solar zenith = 18°, (b) solar azimuth= 
276°, solar zenith = 47°. South Is 6°. Maximum view 
zenith angle is 60° (outer rim) and nadir view is at 
the center. 
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UNDERSTANDING THE RADIANT SCATTERING BEHAVIOR OF VEGETATED SCENES 

Daniel S. Kimes 
NASA/GSFC 


Fundamental knowledge of the physics of the scattering behavior of 
vegetation as gained in this study is important. This growing body of 
knowledge will ultimately serve the remote sensing and earth science 
community in many ways. For example, it will provide (1) insight and 
guidance in developing new extraction techniques of canopy characteristics, 

(2) a basis for better interpretation of off-nadir satellite and aircraft 
data, (3) a basis for defining specifications of future earth observing 
sensor systems, and (4) a basis for defining important aspects of physical 
and biological processes of the plant system. Such fundamental knowledge 
is very important to the long term advancement of remote sensing and earth 
science programs. The overall objective of the three-year study is to 
improve our fundamental understanding of the dynamics of directional 
scattering properties of vegetation canopies through analysis of field 
data and model simulation data. The specific objectives were to (1) collect 
directional reflectance data covering the entire exitance hemisphere for 
several common vegetation canopies with various geometric structure (both 
homogeneous and row crop structures), (2) develop a scene radiation model 
with a general mathematical framework which will treat 3-D variability in 
heterogeneous scenes and account for 3-D radiant interactions within the 
scene, (3) conduct validations of the model on collected data sets, and 
(4) test and expand proposed physical scattering mechanisms involved in 
reflectance distribution dynamics by anaizying both field and modeling data. 

Using a Mark III 3-band radiometer, directional reflectance distributions 
spanning the entire exitance hemisphere were measured of natural vegetation 
canopies ranging from bare soils to complete with 100% ground cover. NOAA 
7/8 AVHRR bands 1 (0.58-0.68 ym) and 2 (0.73-1.1 ym) were used in data 
collection. The cover types reported were corn, soybeans, grass lawn, 
orchard grass, alfalfa, cotton row crops, pine forests, hardwood forests, 
plowed fields, annual grassland, steppe grassland, hard wheat, salt plain, 
and irrigated wheat. All measurements were taken from the ground except 
the forest canopy measurements which were taken from a helicopter. Leaf 
area index, leaf angle orientation distributions, probability of gap functions, 
leaf reflectance and transmittance, and soil reflectance were measured on 
some of the canopies. 

A three-dimensional radiative transfer model was developed (Fig. 1) 
and is unique in that it predicts (1) the directional spectral reflectance 
factors as a function of the sensor's azimuth and zenith angles and the 
sensor's position above the canopy, (2) the spectral absorption as a function 
of location within the scene, and (3) the directional spectral radiance as 
a function of the sensor's location within the scene. It is one of two 
existing models with such general three-dimensional capabilities. This 3-D 
model was expanded to incTude the anisotropic scattering properties of 
leaves as a function of the leaf orientation distribution in both the zenith 
and azimuth angle modes. The model was applied to complete vegetation 
canopies of various leaf orientations—erectrophi le (mostly erect leaves). 
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planophile (mostly horizontal leaves), spherical (equal probability of all 
leaf orientations), and heliotropic (sun tracking leaves). The model was 
also validated and analyzed for agricultural crop canopies, natural 
grassland canopies and forest canopies. The model was used to explore the 
radiative transfers that take place in these various scenes. 

The analysis of the field data and model simulations yielded significant 
information on the radiant scattering behavior of the vegetation canopies. 

The results showed unique reflectance distributions ranging from bare soil 
to complete vegetation canopies for agricultural crops and natural grasslands. 
Physical mechanisms causing these trends were proposed based on scattering 
properties of soil and vegetation. Soil exhibited a strong backscattering 
peak toward the sun (Fig. 2). Complete vegetation exhibited a "bowl" 
distribution^with the minimum reflectance near nadi r (Fig. 3). Sparse 
vegetation canopies showed shifting of the minimum rehectance off of 
nadir in the forward scattering direction because both the scattering 
properties of the vegetation and soil were being observed. In addition, 
sparse canopies exhibited the strong backscatter peak of the soil at 
relatively small solar zenith angles. The largest variations in reflectance 
with changing view angle occurred at large solar zenith angles for canopies 
with very low vegetation covers. As vegetation densities increased and 
the solar zenith angle decreased, reflectance variations decreased. The 
dynamics of the directional reflectance distributions were analyzed and 
physical principles responsible for the observed dynamics were proposed. 

Past studies have demonstrated that the normalized difference transformation 
[AVHRR (Band 2 - Band l)/(Band 1 + Band 2)] is useful in monitoring green 
vegetation biomass. A difficulty in utilizing AVHRR data that scans out 
to 56° is that the signal for any particular target can change significantly 
with various view angles. It was demonstrated that this transformation 
generally decreased the directional variation of the signal. However, 
there were exceptions. For each remote sensing application the user should 
be aware of these variations for the specific cover types being studied, 
solar zenith angle and scanning direction of the sensor with respect to 
the solar azimuth. It was found that complete canopies with different 
leaf orientation distributions (erectophil e, spherical, planophile, and 
heliotropic) had unique reflectance distribution characteristics. The 
dynamics of these distributions were physically explained by directional 
scattering effects of two mechanisms. The first mechanism causes the 
characteristics "bowl" shape of complete canopies that is— increasi ng 
reflectance with increasing off -nadir view angle for azimuth directions. 

It is caused by shadowing gradients and view projection gradients within 
the canopy. The second mechanism is the primary directional scattering of 
the leaves (phase function) due to leaf orientation, source direction, and 
leaf transmittance and reflectance values. The combination of these two 
mechanisms in the various canopies are responsible for the dynamics of the 
overall shape of the directional reflectance distributions. 

The results showed that the scattering behavior of relatively dense 
forest canopies is very similar to the scattering behavior of agricultural 
crops and natural grasslands. Only in more sparse forest canopies with 
significant spacing between the tree crowns (or clumps of tree crowns) does 
the scattering behavior deviate from homogeneous agricultural and natural 
grassland canopies. This clumping of vegetation material has two effects 
on the radiant transfers within the canopy. (A) It increases the probability 
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of gap to the understory and/or soil layers which increases the influence 
of the scattering properties of these lower layers and (B) It increases 
dark shadows within the scene causing increased backscatter and decreased 
forward scatter to occur relative to the homogeneous case. Both phenomenon 
tend to increase backscatter relative to forward scatter. For typical 
forest canopies, the peak backscatter reflectance can be increased as much 
as 30% relative to an equivalent homogeneous canopy due to phenomenon A and 
35% due to phenomenon B. The combined effect of phenomenon A and B can 
cause typical increases of 65% or higher. It is hypothesized that these 
phenomenon are especially important in sparse conifer forests such as the 
boreal forest which account for 50% of the world's forest area. 

The understanding of the radiant scattering behavior of vegetated 
scenes gained in this studied can lead to intelligent techniques for 
extracting canopy parameters from remotely sensed data. For example, 
estimating the hemispherical reflectance (albedo) of terrestrial surfaces 
is of great importance in studying biospheric and atmospheric processes. 

It is proposed that satellite borne instruments represent the only practical 
means of obtaining global estimates of surface albedo data at reasonable 
time resolution; the problem being how to relate the nadir or directional 
reflectance observations from such sensors to the integrated hemispherical 
reflectance. This study investigated (1) the relationships between 
directional reflectances and hemispherical reflectance and (2) the effect 
of solar zenith angle and cover type on these relationships. The results 
showed that errors in inferring hemispherical reflectance from nadir 
reflectance can be as high as 45% for all cover types and solar zenith 
angles. By choosing a solar zenith angle between 30°-40° the same error is 
reduced to less than 20% in both bands. For both bands a view angle of 60° 
off-nadir and j^90° from the solar azimuth reduces this error to less than 
11% for all sun angles and cover types. A technique using two specific view 
angles reduces this error to less than 6% for both bands and for all sun 
angles and cover types. These techniques may yield considerable dividends 
in terms of more reliable estimation of hemispherical reflectance of natural 
surfaces. 

These fundamental research efforts have improved our physical 
understanding of radiant scattering in vegetation canopies. This has been 
accomplished through the analysis of both model simulation data and field 
directional reflectance data in the visible and near infrared wavelengths. 
The directional scattering behavior of vegetation has been described and 
explained for various solar zenith angles, and various canopies with 
different densities, leaf orientation distributions, and optical properties. 
The work has been in response to providing an intelligent basis for defining 
specifica+'ions of earth observing sensor systems and for inferring important 
aspects of physical and biological processes of the plant system. 
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coordinate system, diffuse solar sources (only one is illustrated) , and 
the direct solar source. The diffuse and direct solar sources are ex- 
tended down to the surface of each cell on the top surface of the cell 
matrix. Any 3-D scene can be represented by "filling" the cells with 
various materials with different densities and scattering properties. 
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MEASUREMENT AND MODELING OF THE OPTICAL SCATTERING 
PROPERTIES OF CROP CANOPIES 

V. C. Vanderbilt 

Laboratory for Applications of Remote Sensing 
Purdue University 


I. Justification of Research 

The amount of sunlight specularly reflected by such plants as sunflower, 
corn, sorghum, wheat, and grass is often so large that canopies of these plants 
appear white instead of green when viewed obliquely toward the sun. 

Specular reflections from the shiny leaves of plants originate at the 
interface between air and the cuticle wax layer. Unlike the diffuse portion of 
the light reflected by a leaf, the specular portion of the reflected light is 
reflected at the first surface it encounters; it never enters the leafc From 
the Fresnel equations of optics, the light reflected by such a shiny surface is 
polarized. The polarized portion of the reflectance may be species dependent 
and related potentially to the physiological status and development stage of the 
canopy - to such botanical variables as leaf age and plant water status and 
temperature regime. 

All of this suggests that remotely sensed polarization measurements of a 
plant canopy will contain information about the leaf surfaces — information 
independent of that already identified in the light reflected from the interior 
of the leaves measured in the various spectral regions. The information is 
independent because the polarized portion of the reflected light does not enter 
the leaf to interact with cell pigments, walls or water. 

Polarization measurements, although certainly a potential source of useful 
information about the earth, have never been acquired routinely by satellite- 
borne Imaging sensors as part of earth observation research. The reason is 
simple; the research area is virgin. Hard evidence in the form of physically- 
based theories supported by actual data has only recently begun to demonstrate 
the actual — not merely potential — information in such data. More evidence 
is needed before a dedicated satellite-borne sensor system would be justifiable. 

II. Research Objectives 

The overall objective is to investigate the potential information in the 
polarization data of both single leaves and plant canopies. This research is 
measuring, analyzing, and mathematically modeling the specular, polarized, and 
diffuse light scattering properties of several plant canopies and their 
component parts (leaves) as a function of view angle and illumination angle. 

The potential of these bidirectional radiation properties for ground cover 
discrimination and condition assessment is being evaluated. 

III. Approach 

The research approach has proceeded in both empirical and theoretical 
directions. Measurements performed at our laboratory demonstrate the 
relationship between polarization data and various optical and botanical 
properties of both pieces of foliage and plant canopies. The data provide a 
basis for gaining fundamental understanding of how light is scattered and 
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polarized by a plant canopy. A mathematical model has been developed for 
predicting the specular and polarized light scattering properties of plant 
canopies. By exercising the model we have developed better understanding of the 
potential information in polarization data. 

The empirical part of the approach has involved demonstrating our new 
technique for determining specular, polarized, and diffuse components of the 
reflectance factor of both leaves and plant canopies. Applying this new 
technique to our field measurements, we have determined the specularly scattered 
and polarized light from a plant canopy as a function of view and illumination 
directions. To these reflectance data, we have appended the ancillary data of 
the plant canopies for use by ourselves and other investigators developing and 
testing light-canopy interaction models, such as our specular 
reflection/polarization model. 

A polarization photometer was developed to investigate the potential 
information in the mean and standard deviation of the polarized/diffuse 
components of the reflectance factor of leaves- ..lea^ured Individually in Vivo at 
six wavelengths at the Brewster angle. Measur^mert s wore made (1) in a survey 
of plant species and varieties representing crops, forests, "weeds," and 
horticulture, including as factors (when appropriate) xaaf pigmentation, 
development stage, position of leaf on plant, and position of instrument on 
leaf, (2) of a corn canopy as a function of its moisture stress, and (3) of a 
greenhouse-grown wheat innoculated with wheat rust. 

IV. Results 

We present three results of our research, (1) a theoretical model of the 
specular-polarized light scactered by a plant canopy, (2) data demonstrating a 
linear relationship between the relative water content of a leaf and its non- 
polarized reflectance, and (3) data demonstrating the importance of the 
variable, angle of incidence, in explaining the polarization of light scattered 
by plant canopies. 

A. Theoretical Model 

A model was developed for predicting the amount of light specularly 
reflected and linearly polarized by the leaves of such plant canopies as wheat, 
corn, and sorghum. The model is based upon the morphological and phenological 
characteristics of the canopy and upon the Fresnel equations which describe the 
light reflection process at an optically smooth boundary separating two 
dielectrics. 

The theory demonstrates that, potentially, measurements of the linearly 
polarized light from a plant canopy may be used as an additional feature 
discrimination. Examination of the model suggest that, potentially, satellite 
polarization measurements may be used to monitor plant development stage, leaf 
water content, leaf area index, hail damage, and certain plant diseases. The 
modeling results show that the angles of the polarization analyzer on a 
radiometer or satellite-borne sensor measuring a ground scene may be predicted 
from the view and illumination directions. 

Applicability of the model of the canopy specular reflectance should extend 
to many species because leaves - which specularly reflect sunlight - are 
ubiquitous, unconfined by geography or climate. 

The modeling results. Fig. 1, show that for the predictions from the model 
the single variable, angle of incidence (on the leaf) of the sunlight specularly 
reflected to the sensor, explains much of the variation, as a function of view 
direction (both zenith and azimuth), for a plant canopy with a uniform 
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(spherical) leaf angle probability density function. For examr>le, at any 
specific angle of incidence, 30 degrees, for example, the predicted R^, Fig. 1, 
changes little for zenith view angles less than about 40 degrees, only changing 
significantly for angles approaching 90 degrees. The physical interpretation of 
this prediction. Fig. 1, is that at any given angle of incidence there are more 
specularly reflecting leaves in the field of view of the sensor at large zenith 
view angles near 90 degrees than near 0 degrees. For example, for the sun on 
the horizon in the east and for view directions varying from 0 degrees up to the 
horizon to the north, the angle of incidence will be constant at 45 degrees in 
all these view directions; however, the number of specularly reflecting leaves 
in the field of view of the sensor will be least at 0 degrees zenith view angle 
and greatest at the horizon. 

B. Leaf Moisture vs. Reflectance 

To investigate if relationships exist between reflectance and moisture- 
stressed vegetation, measurements were acquired on the leaves of corn, grown in 
a field under moisture-stressed conditions. (A physiological disturbance of a 
plant usually results in an increase of reflectance in the visible portion of 
the spectrum.) Twenty-four hours prior to reflectance measurements, a portion 
of the field was flood irrigated. Leaves for reflectance measurements came from 
(a) this irrigated treatment, (b) the field-grown, stressed conditions, and (c) 
leaves excised from plants subjected to rapid desiccation. The water status of 
each individual leaf was quantified by measuring the relative water content of 
each leaf sampled for leaf reflectance measurements. 

Fig. 2 depicts the reflectance in the red wavelength band (650 nm) as a 
function of relative water con* .i, '*C). The results. Fig. 2a, show that when 

the RWC of these corn leaves i. . .. , the reflectance factor, R, tended to 

decrease. But within this range of relative water contents sampled (between 50 % 
and 100^), the variaton in the reflectance factor is great - making the 
usefulness of the reflectance measurements for predicting the relative water 
content slight. Fig. 2b shows there is no relationship between the polarized 
component of the reflectance factor, R_, and relative water content. 

But there is a relationship between the non-polarized component of the 
reflectance factor and relative water content, Fig. 2c, which increases linearly 
(R = 0.77) with decreasing relative water content. This relationship appeared 
valid even for relative water contents greater th’,n 80^, a moisture regime for 
which other investigators have found hemispherical leaf reflectance * ■> be a poor 
estimator of relative water content. The ron-polarized component of t’ j 
reflectance factor thus appears to be a hatter predictor of relative water 
content than the reflectance factor. 

C. Plant Canopy Angle of Incidence 

In field experiments designed to provide comparisons between the model 
predictions and polarization data, we measured two plant canopies, each in a 
variety of view directions as the sun moved, providing illumination in a 
continuum of directions. More than 200 spectra were acquired on two wheat 
canopies in the boot (preheaded) and dough (headed) stages of development, and 
on each date in each of 33 view directions. 

The degree of linear polarization (Fig. 3) at a wavelength of 0.66 um is 
plotted for 19 June and 17 July for four view zenith angles and angle of 
incidence, gamma. Regardless of zenith or azimuth view angles, the data points 
for June 19 fall within a narrow region defining an arc. On July 17 the scatter 
in the data^is generally greater, although those data acquired at a view zenith 
angle of 60^ define a fishhook-shaped curve. 
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Fig. 3 shows that most of the variation on June 19 in linear polarization as 
function of the two variables, view zenith and view azimuth angles, is explained 
by the single variable, angle of incidence, a prediction of the 
specular/polarization model. The angle of incidence is computed knowing that a 
small area of shiny leaf must be uniquely directioned to speculcj.’ly reflect 
sunlight to an observer. For the headed wheat of JUly 17, the angle of 
incidence explains the variation in the linear polarization at a zenith view 
angle of 60 (albeit the arc-shaped relationship of June 19 is a fishhook on 
July 17) but less well at smaller zenith view angles. The decrease in the 
degree of linear polarization at large angles of incidence for 6o” zenith view 
angles on July 17 (the hook of the fishhook) is due to the heads (poor specular 
reflectors) decreasing the visibility of the flag leaves to the radiometer. 

V. Conclusions 

The specular reflection process has been shown to be a key aspect of 
radiation transfer by plant canopies. Polarization measurements have been 
demonstrated as the tool for determining the specular and diffuse portions of 
the canopy radiance. The magnitude of the specular fraction of the reflectance 
is significant compared to the magnitude of the diffuse fraction. Therefore, it 
is necessary to consider specularly reflected light in developing and evaluating 
light-canopy interaction models for these two wheat canopies. Models which 
assume leaves are diffuse reflectors correctly predict only the diffuse fraction 
of the canopy reflectance factor. The specular reflectance model, described 
here, when coupled with a diffuse leaf model, would predict both the specular 
and diffuse portions of the reflectance factor. The specular model predicts and 
the data analysis confirms that the single variable, angle of incidence of 
specularly reflected sunlight on the leaf, explains much of variation in the 
polarization data as a function of view-illumination directions. 

Design of hardware to remotely sense the polarization of the light reflected 
by a canopy under a clear sky is simplified by the results of this research. 
First, the lack of fine structure in wavelength in the polarization spectra 
suggests that a design with a single wavelength band covering the entire visible 
wavelength is a possibility. Second, the angle of the polarization analyzer on 
the polarization sensor does not depend on the data but solely on 
view/illumination directions and can be set prior to data acquisition. 
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GEOMETRIC-OPTICAL MODELING OF A CONIFER FOREST CANOPY 


Alan H. Strahler 
Principal Investigator 

Hunter College 
City University of New York 


The primary objective of this research is to explore 
how the geometry of trees in forest stands influences the 
reflectance of the forest as imaged from space. Most plant 
canopy modeling that has been carried out thus far views the 
canopy as an assemblage of plane-parallel layers on top of a 
soil surface. For these models, leaf angle distribution, 
leaf area index, and the angular transmittance and reflec- 
tance of leaves are the primary optical and geometric param- 
eters. Such models, and their derivatives, are now suffi- 
ciently well developed to explain most of the variance in 
angular reflectance measurements observed from homogeneous 
plant canopies. However, forest canopies as imaged by air- 
borne and spaceborne scanners exhibit considerable variance 
at quite a different scale. Brightness values vary strongly 
from one pixel to the next primarily as a function of the 
number of trees they contain. At this scale, the forest 
canopy is nonuniform and discontinuous. 

This research focuses on a discrete-element, geomet- 
ric-optical view of the forest canopy. Trees are considered 
to be solid objects casting shadows on a contrasting plane. 
The brightness of a pixel is then a function of the back- 
ground brightness, the number of trees in the pixel, the 
size and shape of the trees and the individual shadows they 
cast, and the chance overlapping of shadows and tree crowns 
that occurs when trees fall close together. This model of 
the forest canopy presents an alternative approach to canopy 
modeling that is well suited to applications of multispec- 
tral imaging by aircraft and spacecraft. 

The primary objectives of this research are: (1) 
mathematical formulation of a discrete-element, geometric- 
optical model to explain the mean and pixel-to-pixel vari- 
ance in reflectance of sparse to dense forest canopies; (2) 
testing of the formulation through Monte Carlo simulations 
and real image data; and (3) development of an inversion 
procedure to allow remote estimation of size, shape and 
spacing of trees from imagery when the imagery is too coarse 
to permit identification and >neasurement of individual 
trees . 

As developed, the model considers conifers as cones 
(Fig. 1). Parallel-ray geometry is used to describe the 
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Geometric-optical Modeling 


illumination of the three-dimensional cone and the shadow it 
casts on the background. Cones are assumed to be randomly 
placed and freely overlapping; their heights are assumed to 
be ognormally distributed. Both these assumptions are sup- 
ported by literature studies and field v;ork. For the ini- 
tial model, the cones are assumed to be of a fixed shape 
determined by the apex angle; however, in application, the 
apex angle is allowed to vary and an effective apex angle is 
determined for the stand as a whole. 

The research has been carried out in four overlapping 
phases, all of which are nearly complete. The first phase 
involved development of the basic mathematical model 
explaining the distribution of reflectance values expected 
for stands with given size, shape, and spacing parameters 
under fixed conditions of illumination. This phase included 
the development of an important general theorem in geometric 
probability describing the expected variance in random over- 
lapping of shapes in the plane. 

The second phase utilized Monte Carlo simulations of 
forest stands both to confirm the general theorem and empir- 
ically calibrate some parts of the model. As one part of 
the Monte Carlo modeling effort, the bidirectional reflec- 
tance of a pixel consisting of Lambertian cones on a back- 
ground was simulated. The result (Fig. 2) showed many of 
the features of bidirectional reflectance distribution func- 
tions that have been observed for real forests. 

The third phase involved assembling Landsat 80-m and 
SPOT simulation lC‘”m digital imagery for two test sites in 
northern California, and visiting the sites to collect 
ground data for model calibration and verification. The two 
sites contrasted a rather open stand of smaller, broader 
trees (pines) with a dense stand of taller, more narrow 
trees (red firs) . Field measurements included recording the 
diameter and location of every tree within 12-16 circular 
subplots firranged on a grid at each site, as well as the 
collection of accurate heights and apex angles for a random 
subsample of trees within each subplot. 

The fourth phase included the development of an inver- 
sion procedure to estimate stand parameters of tree size, 
shape, and spacing remotely through the observed distribu- 
tion of brightness values for a stand, as well as the check- 
ing of the results against the field observations. For the 
two stands tested, the inversion procedure yielded good 
results. The 80-m imagery was less accurate for the denser 
of the two stands, whereas the 10-ra imagery was generally 
better. The contrasts between the two stands emerged 
strongly at both resolutions in spite of some differences in 
observed and calculated parameter values. 
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The significance of this research is twofold. First, 
it establishes firmly that the three-dimensional geometry of 
ground scenes alone can explain much of the variance in 
digital imagery seen at resolutions typical of present and 
planned sensing instruments. This observation is widely 
applicable beyond forests to many other types of natural 
plant communities, as well as to any other scenes imaging 
three-dimensional objects, such as urban areas, suburban 
developments, etc. Second, it provides a method to directly 
parameterize the shape, size, and spacing of ground objects 
remotely even when such objects are below the level of reso- 
lution of the sensing instrument. The most obvious example 
is the remote estimation of forest biomass, which is a 
direct function of the shape, height, and density of the 
trees. Used in conjunction with quantitative relationships 
between spectral reflectance and leaf area index/standing 
biomass now being developed by other researchers, as well as 
with a stratified sampling procedure, this technique could 
greatly aid global understanding of carbon flows and cycles. 

The future development of this research will have two 
main thrusts. First, more mathematical modeling is neces- 
sary to incorporate new crown shapes and mixtures of shapes, 
as well as to account for terrain effects. Second, more 
field testing is necessary to confirm the utility of the 
inversion procedure and extend the technique to vegetation 
types beyond conifer stands. These developments will not 
only improve our basic understanding of remotely sensed 
scenes, but also allow the direct extension of that improved 
understanding to relevant problems in remote sensing and 
global habitability. 
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CROP CHARACTERISTIC RESEARCH: G. D, Badhwar 

GROWTH AND REFLECTANCE ANALYSIS Johnson Space Center 

Houston, TX 77058 


Justification, Goals, and Objectives 

Much of the early research in remote sensing follows along developing 
"spectral signatures" of cover types. It was, however, found that a 
signature from an unknown cover class could not be matched to a catalog 
value of known cover class. This approach was abandoned and "supervised 
classification" schemes followed. These were not efficient and required 
extensive training. It has been patently clear that data acquired at a 
single time could not separate cover types. 

A large portion of the proposed research has concentrated on modeling the 
temporal behavior of agricultural crops and on removing the need for any 
training data in remote sensing surveys — the key to which is the solution 
of the so-called "signature extension" problem. 

A clear need to develop spectral estimaters of crop ontogenic stages and 
yield has existed even though various correlations have been developed. 
Considerable effort in developing techniques to estimate these variables 
was devoted to this work. 

The need to accurately evaluate existing canopy reflectance model (s), 
improve these models, use them to understand the "crop signatures," and 
estimate leaf area index was the third objective of the proposed work. 

The next section gives a synposis of this research effort. 


Technical Approach 


The technical approach consisted of first developing an accurate model that 
would describe the temporal development of various spectral transforms with 
time, henceforth called a profile, that would depend primarily on crop 
characteristics. It would thus permit features to be extracted from real 
spectral data that describes a specific crop and would not depend on external 
variables such as row direction, sun zenith angle, the atmospheric state, 
etc. Having extracted a very small set of features use the canopy reflectance 
model (s) to gain (a) better understanding and limitations of existing canopy 
models, (b) use the model (s) to accure a deeper physical understanding of 
why these features permit crop separation and develop method(s) of deciding, 
a priori, what features would permit crop separation, and (c) explorate the 
applicability of these models to natural forest communities. 
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Results 


The current research effort has shown that the Kauth-Thomas (K-T) greenness, 
in the spectral space of both the multi spectral spectJ'al scanner or the 
thematic mapper, can be described by a model of the form, 

a/ 2 

/ 23e\ a 

p(t) = pQ + (pm ■ Pq) y — ) " tg) Exp [~3(t - tg)*^] (1) 

where p(t) the K-T greenness as a function of time, p^j the maximum value of 
greenness reached at time of peak greenness 



Pg the value of soil greenness at times at and before emergence, tg, and a 
and 3 are two crop and condition specifics constant. These two constant are 
related to the inflection points of the profile 


ti - tg 


(2a + 1) - y 8a + 1 
4 3 


1/2 


t2 - tg = 


P(2a + 1 ) + \/ 8a + 


, 1/2 
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This research effort has established that the peak greenness above the soil 
line, = Pm " Po> separation 



1/2 


and the time of peak greenness, tp,‘are three characteristics or features 
that carry 95% of all information (Fisher information criteria) available 
in both spectral and temporal data and has led to a drastic reduction in 
the number of variables and a simplification of classifier design. 

Figure 1 shows the power of these features in sep?rating two summer crops, 
corn/soybeans, based on data extracted from the tnematic mapper. The axis 
out of the plane of paper is the number of pixels, the other two axes being 
G|j,ax 3nd a in days. The distributions are more or less gaussian and provide 
excellent separability. Not only has it been shown that these features (G^ja^j 
a, and tp) are applicable to different sensor systems but are truly 
"signature extendable" over vast areas in the United States and, for the 
first time ever in remote sensing, to areas in Argentina. In particular, 
it has been found that, (i) G^g^ (soybeans) > G^g^ (corn), (ii) a (soybeans) < 
a (corn), and (iii) tp (soybeans) > tp (corn). 
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It has been shown that the cardinal points t-|, tp, and t2 are related to 
three specific ontogenetic stages of a crop. This has been shown to be 
case for wheat, barley, corn, and soybeans and appears to be true for 
rice also. This, then permits, a priori, calculation of these stages from 
an agrometeorological model and in complete automation of crop classification. 


It has also been shown that integral, 

t2 

/ P(t)dt, 

tl 

which acts in a manner very similar to the leaf area duration, is strongly 
correlated to yield in the case of both corn and soybeans. More research 
in this direction is called for; ho'-'ever, the potential of a true integrated, 
automatic, and objective crop production system applicable to global corn 
and soybeans areas seem within grasp. 

In order to understand the reasoning behind why G^ax (soybeans) > Gf^a^ (corn) 
almost universely, the existing canopy reflectance models were used. Figure 
2, shows a two dimensional histogram of the G,^ax ^®^f area index, 
calculated using the SAIL model. The distribution for each crop was obtained 
by varying soil types (12 covering the world soil reflectances) and fourteen 
different leaf angle distribution for a total of 168 different calculations 
in each distribution. These calculations show that if the leaf area index 
of these crops is greater than about 4, the primary reason of their 
separability is because corn is C4 plant and soybeans is C3 plant and their 
leaf reflectance and transmittance are intrinsically different. Differences 
in soil type and leaf angle distribution simply require a somewhat larger 
difference in the two leaf area indices to get the same separation. These 
results have also established that if the input to the canopy model such as 
leaf optical properties, soil reflectance, etc. for two crops is known, a 
transform and optimum observation period to provide a given level of 
separability can be calculated, a priori. 

The existing canopy reflectance models had not been subjected to a good 
test, particularly for off -nadir view angle. In order to acquire n degree 
of confidence in the models and their predictions, the homogenous SUITS, 

SAIL, and CUPID canopy reflectance models w-'-'e very carefully evaluated for 
canopies of corn and soybeans. All of these models captured the salient 
features of the canopy reflectance, with the SAIL model showing the best 
overall performance. However, it was found that it is necessary to include 
the specular directional characteristics of leaves into these models. 

Figure 3 shows the performance of the SAIL model before and after the 
inclusion of leaf specular reflectance systematic angular deficiencies of 
this model are substantially removed. This should improve usefulness 
of these models, to evaluate key biophysical parameters, s as '^eaf area 
index. Considerable effort has been devoted in evaluating models for 

coniferous canopies such as Black Spruce, Jack Pine, Red Pine, White Pine, 
and Ponderosa Pine with little success, partly due to lack of input data on 
key input parameters; however, their performance on Aspen and Birch canopies 
appear to be much better. We have also found that the temporal profile of 
an Aspen canopy can be adequately represented by an equation of the form 
(1). The meaning of the t-|, tp, and a would of course be very different. 
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Si qni fica.'c e 

(1 ) For the first time in remote sensing, crop features have been found 
that are truly "signature extendable." In case of corn and soybeans 
they have been shown to be applicable to vast geographic regions of 
the United States for four years 1978, 1979, 1980, and 1982 and are 
extendable to Argentina. Based on SAIL canopy reflectance model, the 
reasons for this applicability have been understood. 

(2) It has been shown that critical ontogenetic crop stages can be estimated 
from spectral data. Based on this work and more detailed work, it is 
suggested that it may be more accurate to estimate crop phenology 
using spectral data than current methods. 

(3) Preliminary evidence suggests that the area of greenness profile from 

ti (when new leaf development stops) to t 2 (dent) is strongly correlated 
to yield. 

The above three results have made it possible to seriously consider an 
automatic and objective crop production system. 

(4) An improved canopy reflectance model has been developed that includes 
the leaf specular component. 

(5) The effectiveness of these models in estimation of leaf area index 
of wheat, corn, and soybeans and recently in study of forest species 
separation and aspen leaf area estimation has been demonstrated. 

Future 


Most of the current work on feature extraction, ontogenetic stage and yield 
has been done on corn and soybean. Some start was made on wheat, barley, 
and oats. The work on spring grains should be intensified. Currently, 
no technique exists to separate the three crops. Additional work on 
ontogenetic stage and yield of corn and soybeans still needs to be done. 

The performance of existing canopy models on forest canopies has been found 
to be sadly lacking, much more so for coniferous forest than for decidous 
forest. Major improvements in these models are called for and a corresponding 
adequate input data set must be collected. Techniques to estimate leaf 
area index or phytomass of vegetation cannot be developed real i ably without 
such an effort. 


75 






MkiShi' 


ii«n*^rtum»»« 




NUMBER OF OBSERVATIONS 







RELATIVE AZIMUTH 
(DEGREES) 


n 


FIGURE 3 


N86- 13801 


GROUND REFLECTANCE MEASUREMENTS FROM THEMATIC MAPPER DATA 


J^ aes E. Conel 
Jet Propulsion Laboratory 
Pasadena, California 91109 

1. Background - Satellite radiance data are measures of solar radiation 
that has been reflected by the Earth's surface and scattered and absorbed by 
atmospheric gases and aerosols. Of concern to geologists are the surface 
reflectance and the degrading effects on surface resolution and albedo 
contrast introduced by atmospheric phenomena. The objects of the presently- 
described research have been to: (1) provide an empirical relationship 

between scanner radiance and ground reflectance allowing interpretation of the 
satellite data in terms of the surface parameter, (2) assess the precision 
with which surface spectral reflectance may be recovered from Landsat-4 TM 
data in the presence of perturbing atmospheric and instrumental factors. Our 
approach is field-oriented, and utilizes ground observations of surface 
spectral reflectance with portable spectrometers and radiometers to develop 
the required empirical relationships. 


2. Method - For a satellite scanner system over locally uniform ground 
with homogeneous atmosphere above the upward directed radiance is given by 
(Chandrasekhar, 1960; Gordon, 1976; Pearce, 1977). 


I = 


I + -- — 
^1 1-Rs 


( 1 ) 


where the first term (Ip represents radiation from the atmosphere alone, and 
the second term radiation that Interacts both with the atmosphere and the 
surface. R is the uniform Lambertian surface reflectance, s is a parameter 
describing the overall probability of backscatter from the atmosphere after 
reflection by the ground, and I2 the overall probability of transmission by 
the atmosphere after reflection by the ground. All quantities in (1) are 
functions of wavelength. The determination of R from equation (1) requires 
measurement of the functions I„, and s. In principle this may be carried 
out by theoretical study of model atmospheres, provided information is 
available on concentration and distribution of aerosols as well as the 
absorbing gas species present. These data are not however generally available 
for the times and places of satellite data acquisition. 

We illustrate here an alternative method wherein the required parameters 
of equation (1) are determined empirically using ground-based measurements of 
spectral reflectance. The image data set used is the TM scene for Wind River 
Basin, Wyoming (acquired November, 1982). The scene covers an area of 32,400 
km^ of central Wyoming. The atmosphere is largely free of clouds and haze, 
which is typical of scenes ordinarily used for geologic remote sensing. Snow 
blankets higher parts of mountains surrounding the basin. 

3. Results for Surface Spectral Reflectance - Field spectral 
measurements of surface bidirectional reflectance were made for selected 
targets throughout the Wind River TM scene using a portable reflectance 
spectrometer (Conel, et al., 1985). Image radiance values, expressed in 
digital number (DN) vrere determined for each of the field sites. Scatter 
plots of image DN versus surface reflectance averaged over bandpasses of the 
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TM were prepared and are shown in figure 1. After conversion of the DN to 
energy units the term is obtained as the intercept value at R = 0, and the 
value of I 2 from the slope. These plots are found to be linear with 
correlation coefficients of 0.96 or greater for all channels of the TM. It is 
also possible to obtain good fits using parabolic equations. Quasi-parabolic 
functional behavior would Indicate a contribution to the radiance from the 
multiple reflection factor 1/(1-Rs) and would provide a means for 
determination of s directly from the curvature term. Based on statistical 
analysis of these various possible fits to the data, the presence of nonlinear 
terms in equation (1) cannot be established. For the atmospheric conditions 
prevailing during time of satellite data acquisition we conclude that the term 
s la negligible compared to unity, and that the linearized form of equation 
(1) is appropriate. 

To Illustrate application of the method for recovery of ground 
reflectance outside the calibration areas we present a comparison between 
image-derived spectra and measured field reflectance spectra for five areas 
throughout the Wind River TM scene (Figure 2). The chronology of events for 
gathering the data was as follows: (1) the original TM scene was obtained in 
November, 1982, (2) field measurements were obtained in November of 1983 and 
used to construct the calibration lines of Figure 1, (3) field measurements 
using portable spectrometers and a hand-held broadband radiometer of the 
additional "unknown" areas were obtained in July, 1984. Thus the comparisons 
Involve measurements spanning a period of approximately two years. Agreem(int 
between the two sets of observations is excellent to good. It is expected 
that the comparisons could be improved if the time interval between satellite 
and field measurements could be shortened. TWo additional factors undoubtedly 
contribute to the differences observed. First, field sites were chosei that 
were homogeneous and (except for Riverton High School lawn) unvegetated. 
Despite precautions the best natural targets are always inhomogeneous at the 
level of a few pjercent reflectance, and we must rely upon a limited sampling 
to provide representative reflectance data. Second, some targets are small 
and it is often difficult to locate them accurately in the images. Factors 

that are likely to contribute second order effects Include those of so-called 

adjacency effects, phatrs effects in the surface reflectance, and the 
comparison of bi-directional versus Lambertian surface reflectance 
functions. Adjacency effects reduce contrasts between contiguous areas of 
differing surface reflectance. Kaufman and Joseph (1982) have provided some 
numerical examples. These problems can be minimized by making measurements 
far from the edges of large areas of differing albedo, but this is not usually 
possible in practice. For optically thin conditions (optical depth on the 
order of 0.2 or less) the effects may amount to a relative change in 

brightness of a percent (see Kaufman and Joseph, 1982) . Relative effects of 

this magnitude will be difficult to detect (see below) let aLone correct for 
in image data. These effects may influence both the calibration relationships 
as well as the determined brightness values for isolated pixels in the 
scene. Phase effects in the surface reflectance can be minimized by taking 
field measurements at the same phase angle and solar elevation as characterize 
the satellite observations. The third problem of interpreting satellite- 
determined reflectance as bi-directional is not considered serious for the 
present observations since the atmosphere, especially at longer wavelengths, 
was optically thin. 
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Considering the possible complications, agreement between measured and 
satellite-determined reflectance properties is remarkably good. This lends 
support to the assumption, implicit in Equation 1, of homogeneity of the 
atmosphere over large areas of the scene during the time of satellite data 
acquisition. The measured surfaces also appeared to have been stable over the 
time span of satellite and field data acquisition. 

4. Sensitivity of the TM-derived measurements of reflectance - By 
differentiating Equation (1) with respect to R (Gordon, 1976) a basis is 
provided for determining what differences in spectral reflectance AR may be 
obtained from IM scanner data, including effects of both atmosphere and 
precision of scanner irradiance measurements AI. The TM system slgnal-to- 
noise ratio (i.e., I/AIjj for all visible and near infrared bands can be 
expressed as a linear function of the irradiance I (Barker, et al., 1983), 
i.e., 1/ Alfj = A + BI where the A and B are known constants, and AI^ is the 
noise-equivalent Irradiance of the scanner. An expression for AR Including 
the signal/noise ratio can then be written down. Figure 3 shows the results 
of calculations for such noise-limited reflectance differences for all bands 
of the TM scene of Wind River Basin. The values of and I 2 ate derived from 
the data given in Figure 1. Based on these calculations, the intrinsic 
sensitivity of the TM system to detection of changes in reflectance everywhere 
in the scene less than about 1%. 

Tfie actual detectable reflectance differences may be limited by 
instrumental gain. These differences AR' are given by AR' = l/(dDN/dR), which 
is the reciprocal slope of a calibration line in Figure 1. The values of AR' 
obtained in this way for the calibration data in Figure 1 are comparable to 
those obtained accounting for the measured pre-flight signal/noise ratio of 
the TM system together with the atmospheric attenuation factors (Figure 3). 
The present gain settings of the TM system are thus consistent with the actual 
noise characteristics of the TM. 

5. Significance of Results - We have employed an approximate analysis of 

the real terrestrial atmosphere and surface (embodied in Equation 1) to 
investigate the question of recovering surface reflectance from TM radiance 
data, and to estimate the uncertainties present in such determinations. It 
was shown that the surface reflectance could in the best examples be recovered 
to within a few percent, and it is expected that these estimates could in 
principle be Improved by better sampling procedures. The theoretical limit on 
such determinations was shown to be about a percent for the TM system. 

Equation 1 and its empirically determined constants can be used in conjunction 
with sharp edges in the image to estimate the atmospheric portion of the 

system modulation transfer function (MTF) . 1^ and I 2 in Equation 1 can be 

thought of as functions of the optical depth and phase function for a 
homogeneous model, and the numerical values obtained as slope and intercept 
used to estimate these parameters for an optically thin atmosphere. These 

determinations have been made for the present data set, but space has not 

permitted a discussion of the results. They along with other results will be 
set out in forthcoming papers. It also appears possible to provide 
independent estimates of optical depth and scattering albedo using the so- 
called "two-halves" method of Kaufman (1982), which will provide an 
interesting comparison with the present results. 
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6. Future Research - Despite erapliasis placed on interpretation of 
satellite-acquired scanner data, aircraft scanner systems like TMS and AIS 
continue to provide essential data sets for the study of many geologic and 
geobotanical problems. Ihe aircraft observations provide useful high spatial 
resolution and (for TMS) encom?^:;:.- a vjide range of surface viewing 
directions. We will study the problem of correcting aircraft data for 
atmospheric effects, and provide metho'’. ; for looking at directional properties 
of the surface independent of the atiuosphere for one or two test sites. The 
approacn will emphasize "field" detemlnation of atmospheric parameters along 
lines indicated above, together with observation of directional atmospheric 
scattering from the aircraft data themselves as inputs to modeling programs. 


7. References 

Baker, J.L., D.L. Ball, K.C. Leung, and J.A. Walker, 1983, Abstracts of Early 
Results Symposium, NASA/Goddard Space Flight Center, Greenbelt, MD. 

Chandrasekhar, S., 1960, Radiative Transfer , Dover, New York. 

Conel, J.E., H.R. Lang, E.D. Paylor, and R.E. Alley, 1985, IEEE Trans. Geoscl. 
and Remote Sensing , in press. 

Gordon, H.R., 1976, Appl. Opt. 15(8), 1974-1979. 

Kaufman, Y.J., 1982, J. Geophys. Res., 87(C6) , 1437-1447. 

Kaufman, Y.J., and J .H. Joseph, 1982, J. Geophys. Res., 87(C2), 1287-1299. 

Pearce, W.A., 1977, Rep. 004-77 , EG&G/ Washington Anal. Serv. Center, 

Rlverdale, Maryland, 136 pp. 


8. Figure Captions 

Figure 1: DN vs. Reflectance Calibration Curves. TM-4, Wind River Basin, 

Wyoming, Nov., 1982. 

Figure 2: Comparison of TM-4 and Field Spectra, Wind River Basin, Wyoming. 

Figure 3: Noise-Limited Uncertainty 4R Derived from m-4. Wind River Basin, 

Wyoming . 
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Thermal Infrared Geologic Remote Sensing Research 
at the Jet Propulsion Laboratory 

Mary Jane Bartholomew and Anne B. Kahle 
January 29, 1985 


JUSTIFICATION 

Remotely sensed thermal Infrared spectral data have great potential to 
improve rock type discrimination if the factors that control thermal infrared 
spectral reflection and emission can be better understood. Improved rock type 
discrimination in turn leads to better and more efficient geologic mapping 
which is carried out in the exploration, assessment and documentation of 
mineral resources and geologic hazards. Geologic maps also contribute 
significantly to the understanding of the natural history of the earth and 
neighboring planetary bodies. 

OBJECTIVES 

The particular goal is to study the spectral emission and reflection 
behavior of rock materials (rocks, rock weathering products and soils) as a 
function of earth-like environment and occurance. Some of the important 
environmental variables that influence the earth's near surface thermal regime 
and hence its thermal emission are sky temperature, sky field of view, 
microclimatology, time of day and season of the year. In addition, chemical 
and {Aiysical factors, mineralogy, {article size, surface roughness, presence 
of desert varnish and rock weathering, play a role. 

APPROACH 

The approach involves both field and laboratory investigations. The 
field work is divided between airborne and ground based studies. Remotely 
sensed mid-infrared spectral data of the earth can be acquired by the airborne 
Thermal Infrared Multis pectral Scanner (TIMS) (Kahle and Goetz 1983; Palluconi 
and Meeks, 1985). TIMS measures normal spectral radiance of the earth's 

surface. It has six channels from 8 to 12 and is flown in a NASA Lear 

Jet. (Within this range the transmission of the earth's atmosphere is good 
and intense absorption features are present in surface silicate materials.) 
Computer image enhancement research is a significant part of TIMS 

investigations . 

The ground based field work is carried out with the Portable Field 

Emission Spectrometer (PFES) which is used to make normal spectral radiance 
measurements from 3 to 14 ijra. This JPL-built instrument is mounted and 

carried by backpack. The PFES has a spectral resolution of approximately 
.2 pm which makes it an important tool that can be used in interpr. ting TIMS 
data . 


However, some scientific questions in this subject require the kind of 
experimental control that can only be provided in the laboratory. For 
example, the influence of earth-like near surface thermal gradients, desert 
varnish, rock we hering, particle size and surface roughness on the spectral 
emission of rock materials needs to be understood and can only be studied in 
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an environment that is free from variations in sky temperature, 
mlcrocllmotology , and sky field of view and uneffected by the season of the 
year or the time of day. 

Iherefore, we are developing th*’ capability to measure in the laboratory 
the normal ipectral radiance or rock materials at ambient temperatures (•' 
20°C) . Ihis involves the careful design and construction of a sample 
accessory, called the emissivity cold box, for a commercial Fourier Transform 
infrared spectrometer that allows control and stabilization of sample 
temperature while shielding it from background radiation, (Brown and Young, 
1975; Aronson and Ensile, 1973; Logan and Hunt, 1970; Conel, 1969; Low and 
Coleman, 1966; Lyon, 1965). A schematic diagram of the sample accessory is 
shown in Figure 1 . 

Liquid nitrogen will be used to control the temperature of the walls of 
the cold box and dry gaseous nitrogen will control the sample temperature. 
The temperature of the gaseous nitrogen will be controlled by a liquid bath. 
A separate liquid bath will be used to control the temperature of gaseous 
nitrogen that purges the spectrometer's optic bench. 

The cold box will take the place of the spectrometer's source mirror 
during the measurement of sample radiance. Modifications to the spectrometer 
that allow this are already complete. Furthermore, the spectrometer's 

standard pyroelectric detector has been replaced by a more sensitive liquid 
nitrogen cooled raercury-cadmium-telluride detector. 

With this configuration the sample radiance at the detector can he 
characterized as 

L(X,T^) = K(X) [e(X) W(X,Tj) + R(X) + S(X)] (1) 

where L(X,T^) = radiation at detector 

Tj = temperature of sample 

K( X) = instrument resonse function 

e( X) = emissivity of sample 

W(X,T|) = black body radiation at temperature Tj 

R(X) = radiation reflected by sample 
S(X) = background radiation 

then by making these additional measurements, 


L(X,T2) = 

K(X) 

[e(X) W(X,T2) + R(X) + S(X)] 

(2) 

B(X,T^) = 

K( X) 

[W( X,Tp + S( X)] 

(3) 

B(^.TlN2) = 

K(A) 

S(X) 

(4) 
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where 


^2 - temperature of sample different from Tj 

^LN2 ~ temperature of liquid nitrogen 

B(X,T]^) = radiation at detector from blackbody at 

temperature Tj 

the sample eralssivity can be calculated as follows. 

L(X,T^) - L(X,T2> W(X,T^) 

"" B(X,T ) - B(X,T_ ') ’ ' W(X,T ) - W(X,T )' (5) 

1 LN2 1 ^ 

All values in equation 5 are either measured or calculated from Plank's Law 
once the temperature of the sample is known. 

Two techniques will be used to determine the sample's temperature. 
First, thermocouples will be located on or near the top and bottom surfaces of 
the sample. Second, the surface temperature will be computed from the 
radiance at the sample's Christiansen frequency (Conel, 1969; Aronson et al., 
1969). 


RESULTS 

An example of TIMS data from Death Valley, California is given in Figure 

2. These images show variation in emittance across part of the Panamint 

Mountains, valley floor and part of the Funeral Mountains. The light areas 
indicate high emittance and the dark low emittance. The values upon which the 
images are based were calculated from Plank's Law and TIMS data. These data 
have been used to make lithologic maps of bedrock and alluvial fans in Death 
Valley (Kahle and Goetz, 1983; Gillespie, Kahle and Palluconi, 1984). 

Figure 3 shows PFES data from Death Valley. The silicate absorption 
doublet from 8 to 10 pm can be clearly seen in the quartzite spectra. This 

feature is only hinted at in the spectra of the volcanic rocks. The presence 

of the absorption features at 6.5 and 11.5 pm clearly Indicate the carbonate 
nature of the dolomite example. 


FUTURE PROSPECTS 

We will evaluate the geologic utility of multls pectral thermal infrared 
data using TIMS, and we will determine the lithologies that can be effectively 
discriminated on the basis of TIMS measurements, either alone or in 
combination with other remotely sensed data sets. Our laboratory studies will 
be directed towards understanding the spectral emission behavior of Earth 
surface materials so as to optimize the interpretation of TIMS data. 
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1. Justification and Objectives 

Among the tools used in passive remote sensing of Earth resources in the visi- 
ble and near-infrared spectral regions are measurements of spectral signature and 
bidirectional reflectance functions (BDRFs). Determination of surface properties 
using these observables is complicated by a number of factors. Including (a) mixing 
of surface components, such as soil and vegetation, (b) multiple reflections of 
radiation due to complex geometry, such as in crop canopies, and (c) atmospheric 
effects. Differences in surface conditions, spectral coverage, viewing and illumi- 
nation geometry, and atmospheric conditions make comparison of observations obtained 
in the laboratory, in the field, and from aircraft and spacecraft extremely diffi- 
cult. In order to bridge the diversity in these different approaches, there is a 
need for a fundamental physical understanding of the influence of the various effects 
and a quantitative measure of their relative Importance. In particular, we consider 
scene complexity effects using the example of reflection by vegetative surfaces. 
The interaction of sunlight with a crop canopy and interpretation of the spectral 
and angular dependence of the emergent radiation is basically a raultidimensionai. 
radiative transfer problem. In the first step of the radiation history, sunlight 
is directly and diffusely transmitted downward through the Earth's atmosphere. The 
downward radiation field then interacts with the plant canopy. Each element of the 
canopy (leaves, stems) reflects (and transmits) light according to its intrinsic 
bidirectional reflectance (and transmittance) properties. The complex canopy geo- 
metry, underlying soil cover, and presence of diffuse as well as collimated illumi- 
nation will modify the reflectance characteristics of the canopy relative to those 
of the individual elements. Finally, upwelling radiation emergent from the canopy 
will be directly and diffusely transmitted by the atmosphere to space, and some 
contamination of the spectral signature with radiation reflected from neighboring 
surface regions will occur due to scattering associated with the adjacency effect. 
The intent of this research program is to develop some of the tools needed in. order 
to achieve a quantitative understanding of such phenomena, using a combination of 
experimental and theoretical methods. 

2 . Technical Approach 


The tools which are needed in order to study this radiative transfer problem 
include (a) knowledge of the bidirectional scattering properties of individual 
canopy elements, primarily leaves, as well as of the canopy as a whole, (b) a 
model describing the interaction of radiation with the canopy, and (c) a radiative 
transfer algorithm for describing the interaction of a nonuniform radiation field 
with the atmosphere. In order to obtain realistic leaf bidirectional reflectance 
functions for incorporation into canopy and scene models, we are using laboratory 
spectrogoniometry to .study the intrinsic scattering properties of individual leaves. 
These measurements provide a data set which obviates the need for nonphysical model 
assumptions, such as Lambertian leaf BDRFs. We have been using a goniometer located 
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at the University of California, Los Angeles to obtain spectra of individual leaves 
for a variety of Illumination and viewing geometries. The instrument operates in 
the spectral range 0,4 to 2.5 ym. The illuminating source may be placed at an 
inclination of 0 to 90“ relative to the surface normal and essentially the entire 
range of emission and azimuth angles may be measured, including near-zero, grazing, 
and high phase angles. An important attribute is the ability obtain data out o 
the principal plane of scattering. The leaf BDRPs obtained in this manner will be 
made available to other researchers for incorporation into canopy models. With 
regard to atmospheric effects, including extinction (absorption and scattering); 
addition of path radiance; multiple reflections between atmosphere and ground; and 
adjacency effects (whereby radiation reflected from surface regions not in the 
field of view are scattered into the line of sight), we have been developing compre- 
hensive computational tools for solving multidimensional radiative transfer prob- 
lems. Our technique incorporates a two dimensional spatial Fourier transform of 
the radiative transfer equation, and the resulting expressions for each Fourier 
component of the radiation field are computed. The full 3-D intensity field is 
then reconstructed using the inverse Fourier transform. The mathematical descrip- 
tion is of general form so as to include consideration of non-Lambertian BDRFs. 
Our technique has been specifically designed to overcome some of the approximations 
employed in previous methods and has been implemented on a minicomputer. The method 
is nearly as general as Monte Carlo but has the advantage of providing deterministic 
solutions for the upwelling and downwelling radiations fields at many altitudes, 
zenith angles, and azimuths from a single computer run. 

3. Research Results 


Experimental data have been acquired using the UCLA goniometer on camellia, 
cucumber, and cauliflower leaves. In order to monitor the effect of using cut leaves 
versus leaves still attached to the plant, the reflectance of a cucumber leaf was 
monitored for several hours following clipping. Reflectivity at 0.98 ym increased 
during this period, indicating water stress and a decrease in the strength of the 
water band. To study spectral dependence of leaf BDRF, we obtained spectra of a 
camellia leaf at view angles of 0 and 30 degrees for illumination at normal inci- 
dence. The ratio of these spectra in the 0.4 - 0.9 ym region is displayed in Fig. 
1. The increased noise level at the short wavelength end is due to lower detector 
response and decreased source light levels. This plot shows that the leaf BDRF 
deviates markedly from Lambertian behavior in a spectrally dependent fashion. The 
greatest departure from Lambertian reflectance is observed in the 0.68 ym chloro- 
phyll absorption feature, possibly as a result of a decrease in the amount of 
multiple scattering within the leaf. 

As an example of the success ,vf our multidimensional radiative transfer algo- 
rlthta, we consider the albedo step function problem, in which two semi-infinite 
fields of differing albedo are adjacent to one another. This model is useful in 
that it provides (a) an evaluation of the diffusion of radiation across the sharp 
albedo boundary separating the two regions, and (b) the ability to verify that the 
solution asymptotically approaches the correct result for a homogeneous surface at 
large distances from the boundary. This surface model may be considered to repre- 
sent a coastal boundary, in which a large expanse of low albedo water borders a 
broad land area of higher reflectivity. We compare results obtained using the 
Fourier transform method with Monte Carlo results for an atmospheric model consist- 
ing of a clear atmosphere with a scale height of 8 km and an optical depth of 0.189, 
corresponding to Rayleigh scattering at 0.47 ym. The two half-planes are Lambertian 
with albedos of 0.0 and 0.5, and the solar zenith angle is 40°. This set of condl- 


94 



tions Is referred to as Model 1. Figure 2 shows the ''jults of calculations for 
Model 1 using the Fourier transform technique, along wit i selected points from the 
Monte Carlo results. The total upwelling intensities are shown as the sum of the 
diffuse (multiply scattered) radiation field and the direct field, which consists 
of photons transmitted to space without being scattered. The agreement between 
the two methods appears to be excellent. We also show at large distances from the 
boundary the intensities calculated using the one-dimensional matrix operator 
method for uniform albedos of 0.0 and 0.6, demonstrating that the Fourier transform 
results approach the expected values away from the boundary. In the vicinity of 
the boundary, the adjacency effect is apparent. Due to scattering within the atmo- 
sphere, the upwellitsg intensities on either "ide of boundary are affected by photons 
reflected from the. ^jurface on the opposite side of the boundary and then scattered 
into the line of sight. In effect, the atmosphere "blurs" the boundary with a point 
spread function whose width is a few kilometers. The Fourier transform technique 
is useful in that its underlying mathematical formalism demonstrates how the atmo- 
sphere filters high spatial frequency information from the upwelling radiation 
field. This filtering can be described by the computation of atmospheric modulation 
transfer functions (MTFs). Figure 3 is a comparison of MTFs calculated using our 
technique for a Rayleigh atmosphere and an aerosol laden (hazy) atmosphere, each 
of total optical depth 0.3 and 5 km scale height. These MTFs indicate that the 
more highly forward scattering haze phase function results in a greater observabil- 
ity of the adjacency effect. 

Another area in which the 3-D radiative transfer algorithm is proving useful 
is in the development of a technique for determining atmospheric transmittance 
from surface images acquired at off-nadir view angles. The basis of the method is 
as follows: As shown in Fig. 2, the total upwelling radiance is the sum of a direct 
and diffuse component. Radiation directly transmitted to space is proportional to 
the surface albedo distribution, and a sharp discontinuity in the surface albedo 
gives rise to a similar discontinuity in the direct field. The diffuse field, on 
the other hand, while exhibiting a qualitative similarity to large scale variations 
in surface albedo, exhibits a much smoother character, owing to the filtering effects 
of the atmosphere, as discussed above. The direct field, which is also proportional 
to atmospheric transmittance exp(-x/)i), where t is the optical depth and p is the 
cosine of the view angle, thus has a very different spatial character from the 
diffuse field. If we compute the spatial Fourier transform of upwelling radiances 
within an image, the spatial "smoothness" of the diffuse field means that it will 
make a negligible contribution at high spatial frequencies. Thus, if we eliminate 
low frequency Fourier components we axe left with a measurement that is directly 
proportional to the atmospheric transmission factor exp(-T/p). By obtaining images 
of the same surface region at several view angles, it should in principle be possible 
to solve for x, provided non-Lambertian surface effects can be minimized. Our cal- 
culations suggest that acquisition of images in the vicinity of 60° off-nadir should 
be sufficient to accomplish this and recover atmospheric transmittance with an accu- 
racy of 5-10%. 

4 . Significance of Results and Future Directions 


Many investigators have studied the spectral properties of a wide variety of 
crop samples. The spectral reflectance of the foliage is modified by changes in 
the internal structure of the leaves associated with various forms of stress, such 
as water or mineral stress. The specularity of the reflected radiation increases 
at wavelengths where increased absorption occurs, presumably due to reduced depth 
of penetration of the incident radiation into the leaf resulting in fewer multiple 
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reflections between layers. Therefore, we expect changes in leaf structure result- 
ing from stress to affect the angular as well as spectral reflectance properties. 
Thus, it is important to incorporate realistic (non-Lambertian) leaf bidirectional 
reflectance functions into canopy and scene models. The spectrogoniometric leaf 
measurements will provide such data. 

The Fourier transform radiative transfer algorithm which we have developed is 
a useful tool for studying the effect of the atmosphere on high spatial resolution 
imaging of the Earth's surface. As discussed above, the sharpness of the adjacency 
effect in the vicinity of albedo boundaries is dependent upon the vertical distribu- 
tion of the atmospheric scatterers. We are using our radiative transfer method to 
perform a systematic study of the variation of atmospheric point spread functions 
with aerosol vertical distribution, phase function, single scattering albedo, and 
optical thickness. The results of such calculations will permit evaluation of the 
possibility of observing atmospheric blur directly in high resolution Images of 
the Earth's surface, with the goal of determining aerosol properties. Conversely, 
these calcuilations will lead to techniques for removing spectral contamination 
associated with the adjacency effect. 

Further work is necessary in order to demonstrate tbs feasibility of implement- 
ing the multiple view angle optical depth retrieval technique on a space platform. 
The major uncertainty at present is the manner in which the unknown angular reflec- 
tance properties of the surface should be taken into account. For the model BDRF we 
used in our simulations, off-nadir imaging near 60° appears appropriate to compensate 
for lack of a priori knowledge of the shape of the surface BDRF. Whether this range 
is appropriate for natural surfaces can only be determined by experimentation in the 
field. A number of portable field spectrometers are available for obtaining BDRF 
measurements. The success of this technique would be significant because it provides 
a means of monitoring tropospheric opacity, which is important not only for atmo- 
spheric studies but also because it is one of the most important parameters govern- 
ing the effect of the atmosphere on surface images. An extension of the optical 
depth retrieval technique is to near infrared spectral regions where water vapor 
absorption occurs. In its current form, the technique relies on the assumption that 
atmospheric transmittance variesias exp(-T/p). In the water vapor bands, any finite 
bandwidth spectral channel will Incorporate many individual lines and this law may 
not be followed. In order to study the view angle dependence of transmittance in 
such regions, we are developing a radiative transfer algorithm which permits calcu- 
lation of broadband radiances for channels containing many absorption lines in the 
presence of aerosol scattering. 

Finally, we plan to investigate the applicability of our 3-D radiative trans- 
fer method to scattering in plant canopies. In principle, our algorithm is appli- 
cable to heterogeneous media and has the advantage of providing a true solution 
to Che multidimensional radiative transfer equation. In addition, solutions are 
obtained at many view and illumination angles and azimuths, thus providing the 
full canopy BDRF. Variation of the BDRF of Che individual leaves will permit 
study of the effect of changing reflectance properties of Che individual canopy 
elements on the canopy angular and spectral reflectance as a whole. Such calcula- 
tions will mesh nicely with Che goals of Che experimental portion of this investi- 
gation. 
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Fig. 1. Ratio of reflected Intensity at 30° view angle to 0° view angle for camellia leaf. 
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1. INTRODUCTION 

The Earth’s atmosphere has a significant bearing on the quality of images of the Earth’s 
surface as observed from space by reflected sunlight. This atmospheric effect is caused by absorp- 
tion and scattering of sunlight by m.olecules and aerosols. The molecular effects can be accounted 
for, but the aerosol effects are variable and poorly known during satellite observations. The op- 
tical parameters used to compute the atmospheric effects are; its optical depth to determine the ex- 
tinction of light; the single scattering albedo, which is the fraction of light scattered from total ex- 
tinction; and the scattering phase function giving the probability of scattering in a particular 
direction. 

The atmosphere tends to make dark surfaces appear brighter and bright surfaces darker, but 
the resultant effect depends on the surface reflectance pattern itself. Atmospheric scattering dif- 
fuses the boundaries between adjacent fields. This effect, called the adjacency effect, is due to 
light reflected from the ground surrounding a sensor’s field-of-view and then scattered into the 
field-of-view. The adjacency effect depends on the above mentioned optical parameters and on the 
depth of the aerosol layer. 

Theoretical analyses of atmospheric effects, ground-based measurements of atmospheric 
parameters, and analysis of satellite data are reported on here. This work occurred during the 
period 1982-4. 


2. JUSTIFICATION, OBJECTIVES, AND APPROACH 
2.1. Justification 

Recently developed theoretical models of radiative transfer in the atmosphere can be used to 
perform atmospheric corrections for known aerosol characteristics. The development of correction 
algorithms requires testing the models against measurements to investigate the main atmospheric 
characteristics that affect certain remote sensing techniques, and to develop methods to measure 
these characteristics from ground and from space. 
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2.2. Objectives 


To study the atmospheric effect on remote sensing by theoretical modeling and to measure 

this effect. To relate the effects of the atmosphere to the aerosol characteristics in order to deter- 
mine which are the main characteristics that affect a given remote sensing function. The specific 

goals are: 

a. Laboratory simulation and analysis of field measurement of the atmospheric effect to study 
the effect and to validate theoretical models. 

b. Development of a theoretical model that can simulate the atmospheric effect over any variable 
surface, 

c. Study the atmospheric effect on the spatial resolution of satellite images of the surface and on 
classification of surface features. 

d. Analyze the relative importance of aerosol scattering and absorption on remote sensing of the 
surface in order to determine if atmospheric corrections based on the optical thickness alone 
are feasible. 

e. To measure the atmospheric characteristics that are important for atmospheric corrections: op- 
tical thickness, path radiance, scattering phase function, and single scattering albedo, and find 
the relations between them in order to establish the minimum information required for at- 
mospheric corrections, 

2.3 Approach 

a. In order to measure the atmospheric effect and test theoretical models, a field and controlled 
laboratory experiments are performed in the presence of nonuniform surface reflectance. 

b. The simulation of the atmospheric effect over a general varying surface reflectance is perform- 
ed by the calculation of the atmospheric Modulation Transfer Function (MTF) and application 
of two-dimensional Fourier Transforms for the calculation. 

c. Analyze the aerosol effect on spatial reduction and classification by means of the atmospheric 
MTF and application of the 2-D Fourier Transform technique. 

d. Radiative transfer computations are used to derive the relative importance of aerosol scattering 
and absorption on remote sensing of the surface albedo, vegetation index, and classification. 

e. A system is developed to measure simultaneously from the ground the optical thickness, the 
path radiance, the scattering phase function, and the single scattering albedo. The covariance 
between these parameters is studied. Inversion procedures are developed to derive from the 
observations the aerosol optical-thickness, scattering phase function, and absorption. 
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3. RESEARCH RESULTS 


3.1. Measuremenl of the atmospheric effect 

A laboratory e.xperiment [6] was conducted in which a half white, half black plate repre- 
sented the nonuniforni surface and a dish with hydrosol of latex spheres represented the hazy at- 
mosphere. A solar simulator irradiated the dish and a detector scanned the upward radiance across a 
line perpendicular to the border between the white and black area. Fig. lA shows an example of 
the measured radiance and the theoretical fit [1] as a function of the distance from the boundary 
X (normalized by the hydrosol height H). In this experiment the adjacency effect (the effect of a 
bright object on the radiance above a nearby dark object) was measured as a function of hydrosol 
optical thicknesses for nadir and off nadir observation directions. Good agreement was found be- 
tween theory and experiment. The measured adjacency effect had an amplitude 20% higher than 
the theoretical one, and a range 25% shorter than in theory. 

A field experiment was also conducted [8] in which the upward radiance above and below the 
haze layer was measured simultaneously with ground and airborne measurements of the at- 
mospheric characteristics. Fig. IB shows an example of the upward radiance from the nadir below 
and above the base. The radiances are measured over the dark water as a function of the distance 
from the bright vegetated land (the wavelength is 773 nm). The reflectance of the vegetated surface 
was approximately 0.35. A comparison with theory shows similar good agreement as in the 
laboratory study, with the adjacency effect somewhat stronger and narrower than in theory (as 
was found in the laboratory study). This difference can be due to the approximations in the theory 
and due to errors in accounting for very large particles. 

3.2. Resolution and classification 

The experimental verification of the radiative transfer theory and the application of Fourier 
Transform [7] to generalize the theory so that the upward radiance can be simulated for any 2-D 
varying surface, enables the application of the theory to study the atmospheric effect on spatial 
resolution and classification of surface features [3], [5J and [7]. It was found (see Fig. IQ that the 
atmospheric effect reduces the spatial resolution substantially. For the Thematic Mapper (TM) the 
spatial resolution of 30m (defined for Modulation transfer value of 0.35) is reduced to 100m by a 
hazy atmosphere (aerosol optical thickness = 0.5, X= 550nm). An atmospheric correction that ig- 
nores the adjacency effect can correct the resolution to 50m. Similarly, significant effects were 
found on classification of surface features [3] and on separability between field classes [7]. 

3.3 Scattering versus absorption 

A theoretical study was conducted [9, 12] to test the common hypothesis that the atmospheric 
effect depends mainly on one atmospheric parameter — the optical thickness (t^). An increase in 
the amount of aerosol tends to brighten dark surfaces and darken bright surfaces (see Fig. 2A). In 
the brightening process aerosol scattering has a dominant effect and absorption is of secondary 
importance; but in the darkening process the absorption has an important role. For a given bright 
surface the atmospheric effect can be of brightening or darkening depending on the ratio of 
aerosol scattering and extinction optical thicknesses (w^,). A reliable estimate of such effects can 
not be made until more is known about atmospheric absorption and its variability. However, Fig. 
2B shows that remote sensing of contrast such as vegetation index is weakly affected by absorp- 
tion, but still affected by changes in the optical thickness. The vegetation index is a function of the 
slope of a line in this figure and thus does not change appreciably with absorption (from model I 
to 2 and from model 3 to 4). Figure 2B also shows how the reflectance of the earth-atmosphere 
system is altered by the atmosphere. 
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3.4 Ground-based experiment 


In order to learn about the optical characteristics of the atmopshere as they affect remote 
sensing, a ground experiment is being performed in order .to collect and analyze a substantial data 
base. Simultaneous measurements in the visible and near infrared spectrum are made of the solar 
transmission, the sky radiance in the direction of the North Star and also in the almucantar 
through the sun and the degree of polarization of the skylight. The aerosol optical thickness is 
derived from the solar transmission. The path radiance is measured in the direction of the North 
Star in order to minimize geometric effects and observe changes associated with aerosols. During a 
day the angle between the observer, sun, and North Star is essentially constant. 

Attempts have been made to relate atmospheric corrections to only one parameter — namely 
the aerosol optical thickness, in only one spectral band. Figure 3A shows the strong correlation 
between path radiance (at 870nm) and optical thickness (at 610nm). 94 percent of the radiance 
variance is associated with optical thickness variance for these wavelengths. The correlation bet- 
ween the radiance at 870nm and the optical thickness at 610nm was bigger than the correlation 
with optical thickness at 870nm. This is due to a different size dependence of the radiance and the 
optical thickness. The correlation between path radiances in two spectral bands (Fig. 3B) shows 
the accuracy of estimating the path radiance in one band from that ' ^ another band. The radiance 
at 870nm can explain 74 percent of the radiance variance at 480ns ,sid 62 percent of the radiance 
variance at 1670nm. 

The aerosol scattering phase function is derived from the radiance of the skylight in the solar 
almucantar. The inversion proceeds by an iterative method to account for molecular and multiple 
scattering. An inversion of radiances computed for a model atmosphere is shown in Fig. 3C. The 
true aerosol phase function of scattering angle is given by P^. P would be the estimated phase 
function derived from the radiance, if only molecular scattering was accounted for. P„ is the first 
guess of the aerosol phase function, which is used to estimate the radiance of multiple scattering. 
This value is subtracted from the measured value to obtain the P, -estimate of the aerosol phase 
function. The procedure is repeated with this new value of the phase function (P,). After two 
more such iterations, the curve Pj is obtained, and it is nearly identical to the true value (P 4 ). 
Although the simulated radiances are free of errors, an analysis of the error budget has been 
made. This inversion method has been applied to measured radiances. Eventually, these functions 
will be compared with those derived from solar transmission data. 


5. CONCLUSION 

Radiative transfer theory (RT) for an atmosphere with a nonuniform surface is the basis for 
understanding and correcting for the atmospheric effect on remote sensing of surface properties. 

In the present work the theory was generalized and tested successfully against laboratory and field 
measurements. There is still a need to generalize the RT approximation for off-nadir directions 
and to take into account anisotropic reflectance at the surface. 

The adjacency effect results in a significant modification of spectral signatures of the surface, 
and therefore results in modification of classifications, of separability of field classes, and of 
spatial resolution. For example, the 30m resolution of the Thematic Mapper is reduced to 100m by 
a hazy atmosphere. The adjacency effect depends on several optical parameters of aerosols: optical 
thickness, depth of aerosol layer, scattering phase function, and absorption. Remote sensing in 
general depends on these parameter, not just adjacency effects, but they are not known well 
enough for making accurate atmospheric corrections. It is important to establish methods for 
estimating these parameters in order to develop correction methods for atmospheric effects. Such 
estimations can be based on climatological data, which are not available yet, correlations between 
the optical parameters and meteorological data, and the same satellite measurements of radiances 
that are used for estimating surface properties. Our knowledge about the atmospheric parameters 
important for remote sensing is being enlarged with current measurements of them. 
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formity in the surface reflectance. spatial frequency (c;km) 

Fq is the solar spectral irradiance 
and 0o=arccosine(po). the solar 
zenith angle. 

A - Laboratory simulation for optical thickness 0.3, nadir observation. The 
theoretical fit (....) and the surface reflectance (- - -) are indicated. X is 
the distance from the bright area and H the hight of the scattering medium, 
b - Field measurements for aerosol optical thickness 0.5. The aircraft measurements 
X ) and theory (....) are indicated. 

£ - The atmospheric effect on spatial resolution of the Landsat Thematic Mapper (TM). 
The modulation transfer functions for the TM( — ) alone, with atmospheric effect (....) 
and with correction that ignores the adjacency effect {- - -) are shown. The original 
resolution (30m) is reduced to 50m by the adjacency effect and to 100m by the total 
atmospheric effect. 
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RADIANCE (L) - HAZY DAY 



Fig. 'l\ A - Scatter diagram of the radiance from a 
Landsat Tinage on a hazy day as a function of the 
radiance on a clear day taken over Washington, DC; 
for 700-800nm band. The radiance is normalized to 
reflectance units. 


radiance !660) 



RADIANCE 1820) 


B - Radiance of models of the earth-atmosphere system. 
The surface^consists of dense alfalfa. The solar zenith 
angle is bO ; the direction of observation is given by 
the polar angle 0 ,and azimuth (ji from the direction of 
propagation of incident sunlight. The number besides 
the p lotted points refer to the aerosol models: 

IMODEL co r, TA^Onm) t A(820nm) ! 

I 1 0.% 0.287 0.200 ■) 

I 2 0.88 0.287 0.200 j 

3 0.y6 0.861 0.600 I 






SCATTERING PHASE FUNCTION 




RADIANCE/Fo (870 nm) 



RADIANCE/Fo (1670nm) 


Fig. 3: Example of analysis of ground measurements of the sky radiance. 

^ - The relation between the optical thickness at(610nm) and the sky radiance (at R7nnm) 
in the direction of the North star. Each symbol represents a different observation day. 
iB - The correlation between the sky radiance at 480nm and R70nm; and 87nnm and IRTOnm. 

£ An example of inversion of the almucantar simulated radiances for the aerosol scatterinn 
phase function. 


106 







N86- 13805 


THE CHARACTERIZATION OF SURFACE VARIEGATION EFFECTS 
ON REMOTE SENSING 

William A. Pearce 

EG&G Washington Analytical Services Center, Inc. 
5000 Philadelphia Way 
Lanhara, Maryland 20706 


Improvements in remote sensing capabilities hinge very directly 
upon attaining an understanding of Che physical processes contributing 
to the measurements. In order to devise new measurement strategies and 
to learn better techniques for processing remotely gathered data, we 
need to understand and to characterize the complex radiative 
interactions of the atmosphere-surface system. In particular, it is 
important to understand the role of atmospheric structure, ground 
reflectance inhomogeneity and ground bidirectional reflectance type. 
Our goals, then, are to model, analyze, and parameterize the observable 
effects of three dimensional atmospheric structure and composition and 
two dimensional variations in ground albedo and bidirectional 
reflectance. 

To achieve these goals, we employ a Monte Carlo radiative 
transfer code to model and analyze the effects of many of the 
complications which are present in nature. We can, for example, treat 
finite clouds as well as vertical atmospheric structure. We can model 
the effects of two dimensional structure in ground albedo and 
directional reflectance. We can examine alternate data acquisition 
strategies perhaps involving variations in detector viewing direction, 
solar illumination direction, or polarization sensitivity. Most 
Importantly, for the purposes of analysis and characterization, we can 
divide the detected radiance into components. The radiance transmitted 
directly through the atmosphere from the ground to the detector is the 
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componeat which carries the most significant Information about the 
character of the surface. The radiance arising solely from atmospheric 
scattering is a term which, in effect, should be subtracted from, the 
data as it has no bearing on the surface character. Lastly, there is a 
ground diffuse component comprised of radiation which, having scattered 
from the ground, is atmospherically scattered prior to detection. This 
last terra is responsible for non-local adjacency effects. Its influence 
on the data usually may be described by an atmospheric spread function. 

Specifically, the objectives of the present study are to 
investigate the sensitivity of the radiance components and the 
atmospheric spread function to: variations in the vertical aerosol 
profile; variations in the aerosol size distribution (and, hence, its 
scattering distribution); variations in the detector look direction and 
solar Illumination direction; variations in ground reflectance type and 
bidirectional reflectivity; the presence of clouds in the vicinity of 
the point of view. 


RESULTS 


As the modal radius of the atmospheric aerosol size distribution 
increases, the associated scattering phase function is enhanced both in 
the forward and backward directions. For remote sensing, it is the 

changes in the forward scattering peak which are more important. Our 
simulations show significant sensitivity in the inner regions of the 
spread function (within 250 meters of the point of view). Figure 1 
Illustrates this, showing the atmospheric modulation transfer function, 
or MTF (the Fourier transform of the atmospheric line spread function) 
for three log-normal size distributions with modal radii of 0.13, 0.26, 
and 1.04 micrometers (bottom curve, middle curve, and top curve, 
respectively). 
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For Lambertian ground reflectance, there is relatively little 
sensitivity to detector nadir angle for nadir angles less than about 
40 degrees. For larger nadir angles, however, the atmospheric line 
spread function becomes asymmetric and depends upon the illumination 
angle. Figure 2 shows the atmospheric line spread function for a 
detector nadir angle of 50 degrees (detector pointing roughly into the 
sun), while the dashed curve is for 180 degree receiver azimuth. The 
solar zenith angle for both cases is 47 degrees. The line spread 

function here is computed from the gradient of the intensity in a scan 
across a. fixed albedo boundary. The asymmetry is due primarily to the 
forward— backward asymmetry of the aerosol scattering phase function. 

Variations in the form of the vertical aerosol density profile 
tend to scale the range of the atmospheric spread function. This 

scaling appears to be roughly linear in the average aerosol altitude. 

Clouds in the vicinity of the point of view can influence the 
detected radiance. In particular, the shadow of a cloud produces much 
the same intensity gradient when scanning across its edge as is obtained 
for a high contrast ground albedo boundary. Figure 3 shows a detailed 
comparison of the radiance components for two cases. The solid lines 
show (from top to bottom on the right side of the figure) total 
radiance, ground direct radiance and atmospheric radiance as detected 
from nadir while scanning across a Lambertian albedo boundary (albedos 
0.4 and 0.07). The symbols correspond to a similar scan across a cloud 
shadow boundary (the ground albedo is 0.4, the solar zenith angle Is 
47 degrees, and the altitude of the cloud top is 7 km.). The pluses 

indicate total radiance; the solid squares show the ground direct 

radiance; the open squares show the ground diffuse component; and the 
crosses mark the atmospheric radiance. Here, the atmospheric component 
varies because the cloud shades the scattering atmosphere as well as the 
ground. The ground direct radiance varies less abruptly on the dark 
side of the shadow due to cloud transparency and translucence. The 


109 



for t:h« two Is re«arfeafel:f stmU-Mt though the 

€0apo««jits contrihoting: to th« total: feahav^ rather dif fere«t-;i,y<. 
t'ltaiiC the tattge of the effect of^ the shadow qs% the defeectad ithtehS'lty is 
islgo ll: tea n.t ly- 'large * 


mmtrichfiCE 4is!i,> wti - \jm mm 

ikxt sirBulatlon.s have showtt that the form of the aeroaol scatter'tng 
phase fiiactlOTi is a« tapartant determinant of the: atmospheric spread 
fiwicttoh* For nadir abaervatios,,- the caatral peak of the Spread 

funct'ioa is priiaarlly affected:. For o.hservotioe:®: well away from naditj 
the spread ftmcflots can beco?ae: msymmtrtc even though ground ref lectaoce 
Is haaii&ertlafj*: The shape of die spread funetlon is not very sensitise 

to the vertical proflie of thfe aercfsol. We: .have foumt t.hat elie shadows 
of htfh, dense cloyds: produce- chaages 1 b the' total detected ;radla«e:e 
guite sitsllac In range and laagnlfeo.de with those; as;sociated with a high 
eoatrast ground albedo houndar.y. Hm*; radiance ccaponeats * however, axe 
not so slailar* 

The nex:£ major varlahl.e ta: be s.tt»4:le<i is the sen8itivl:ty of 
observed radiances to varlatio.n8. in ground tiidlrec:t,lc«'ial reflectance and 
to two’~dime;nt loaai .patteegs of varying reflectance type* We wish 
especially to undersfeand. to,- what extent va.riat.loits lo bl;direct;iofial 
refle,ctaa;« can he used to l:Mp:rove class If icat. ion and lderittf|cation' of 
surface areas. Wa wish :to study .the ucillcy of pola.r lEacion senait-lve 
detectors* We will extend o.ui study of the effects .of c'io.uds, 'lo;aklng 
bo'th at high thin elowds oJnd at- the* detatled effects asepciated with .low 
thick :Cl:Oudp,* 


m. 






FIGURE 1. 


ATMOSPHERIC SPREAD FUNCTION 




INTENSITY 


R«Dl«NCeS: CLOUD SHADOW AND ALBEDO BOUNDARY 


0. 080H 


0 . 060 - 


i- 


I 

I 

1 


— -r" 


=E 


0. 0A0! — 


h- 


0 . 0 £ 0 |- 


I 


0. 000L- 


m 


+ 


L I 

' II 

II 




i 




i ;>-‘r ' 




:fc 


A. 


X 


□ 


X 


-£. 00 - 1 . 50 - 1 . 00 - 0.50 0.00 0.50 1.00 1.50 £.00 
DISTANCE FRON SHADOW EDGE ( KN ) 


Cl CLOUD 


113 


FIGURE 3 



N86” 13806 


MULTIDIMENSIONAL MODELING OF ATMOSPHERIC EFFECTS 
AND SURFACE HETEROGENEITIES ON REMOTE SENSING 


S.A.W. Gerstl, C. Simmer, and A. Zardecki 

Theoretical Division 
Los Alamos National Laboratory 

The overall goal of this project is to establish a modeling capability that 
allows a quantitative determination of atmospheric effects on remote sensing 
including the effects of surface heterogeneities. This includes an improved 
understanding of aerosol and haze effects in connection with structural, 
angular, and spatial surface heterogeneities. One important objective of our 
research is the possible identification of intrinsic surface or canopy charac- 
teristics that might be invariant to atmospheric perturbations so that they 
could be used for scene identification. Conversely, an equally important 
objective is to find a correction algorithm for atmospheric effects in satel- 
lite-sensed surface reflectances. 

Our technical approach is centered around a systematic model and code 
development effort based on existing, highly advanced computer codes that w.?re 
originally developed for nuclear radiation shielding applications. 
Computational techniques for the numerical solution of the radiative transfer 
equation were adapted on the basis of the discrete-ordinates finite-element 
method which proved highly successful for one- and two-dimensional radiative 
transfer problems with fully resolved angular representation of the radiation 
field . 

Results of the initial phases of our research project, namely, the 
implementation and verification of computer codes and modern atmospheric data 
bases, have been published and proved highly successful. In the next step, a 
coupled atmosphere/canopy model was developed that allowed us to perform a 
sensitivity analysis and to quantify how satellite-sensed spectral radiances are 
affected by increased atmospheric aerosols, by varying leaf area index, by 
anisotropic leaf scattering, and by non-Lamber tian soil boundary conditions’ 
Our numerical results allow the following conclusions: 

1. Atmospheric perturbations are the major contributors to 
Landsat MSS bands 1 and 2 (visible), while they play only a 
minor role in MSS bands 3 and 4 (near-IR) for remotely 
sensed images of vegetative surfaces. The Kauth-Thomas 
greenness parameter is found nearly proportional to the 
atmospheric optical depth over uniform vegetative surfaces. 

2. The inclusion of anisotropic leaf scattering character- 
istics in our canopy model is significant to explain 
reflectance measurements directly above the canopy but 
produces only negligible effects on satellite-measured data 
above the atmosphere at nadir. 
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3. Non-Lambertian and specular surface reflectance charac- 
teristics (BRDF) affect satellite measurements at nadir only 
insignificantly but can make large differences for off-nadir 
observations . 

The details of these results are also published and led us to our most recent 
investigations of how typical non-Lambertian angular canopy reflectance charac- 
teristics may be used for scene identification. 

Identifiers for horizontally homogeneous vegetated surfaces that may be 
used in remote sensing are the spectral distribution and the angular 
distribution of radiances above the canopy. The spectral distribution is 
determined by the optical properties of the constituents of the particular 
canopy, e.g., leaves and stems, while the angular distribution is created mainly 
by the canopy architecture. While the spectral characteristics of most canopy 
components are very similar for most plants, the canopy architecture .lay be most 
dissimilar. This fact leads to intrinsic canopy angular reflectance patterns 
that may be typical for certain classes of plant canopies. Thus, the 
atmospheric perturbations of the angular reflectance pattern of a Lambertian, a 
mixed Lambert ian/ specular surface, and of the measured BRDFs of savannah and 
coniferous forest canopies were studies using aerosol-free and polluted 
atmospheres . 

Two typical results of our calculations are shown in Figs. 1 and 2, where 
measured angular reflectance patterns above real canopies are computationally 
"transported" through a mildly aerosol-loaded atmosphere to a hypothetical off- 
nadir satellite sensor. These results show that maxima, like the canopy "hot 
spot" (a dominant peak when viewing in the direction of the sun with the sun at 
observer's back), are still detectable above the atmosphere even if observed 
through heavily aerosol- loaded atmospheres. It appears that especially those 
patterns that vary with the azimuth angle are better preserved than variations 
with view zenith angles. Note, in particular, how the shape of the reflectance 
pattern along the 180° azimuth changes when transported through the atmosphere, 
while variations with the azimuth along the 30° zenith remains basically the 
same. The underlying cause for this effect is obviously the constant path 
length through the atmosphere for all directions with the same view zenith 
angle. Thus, strong azimuthal variations typical for any surface or canopy 
reflectance pattern undergo only minor perturbations by the atmosphere, while 
the angular patterns for varying view zenith angles are much more perturbed by 
atmospheric effects. This finding suggests that off-nadir satellite 
observations with varying view azimuth angles may contribute significant new 
data that is expected to be valuable for crop identification. 

In our search for an effective new atmospheric correction algorithm, we are 
presently performing exploratory calculations with the objective of recon- 
structing, from radiance distributions above the atmosphere, a typical angular 
reflectance pattern as measured directly above a canopy. Figs. 1 and 2 
demonstrate clearly the severe atmospheric perturbations that are cumulatively 
added to the upwelling radiance distributions, especially in the visible 
wavelength region and for large view zenith angles. However, we show in Fig. 3 
the result of a simple correction algorithm applied to the coniferous forest 
data (see also Fig. 1) for an aerosol-polluted, rural atmosphere of horizontal 
visual range of 23 km, equivalent to vertical optical depth at 0.55 pm of 0.45 
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for the entire atmosphere. The top half of the iso-radiance contour plots give 
Kriebel's measured data directly above the forest canopy while the bottom half 
shows our reconstructed radiance distributions from values above the atmosphere 
but corrected for atmospheric perturbations. This correction is based on a 
simple subtraction of the atmospheric path radiance component as determined from 
an independent additional radiative transfer calculation assuming purely 
Lambertian surface reflectance but with the same albedo as that of the real 
canopy. For our forest canopy, we have Ot = 11% at 0.85 pm. The comparison in 
Fig. 3 shows an agreement within approximately 10% between the measured and re- 
constructed radiance distributions for view zenith angles lower than 60° and all 
view azimuth angles. The hot spot is almost perfectly retrieved after this 
simple atmospheric correction. These results are both surprising and 
encouraging, because the correction apparently works very well even at the 
visible wavelength where atmospheric effects dominate the canopy signal by far. 

The next logical step in our research project must be a two-pronged 
approach. First, the analysis of our (as yet) simple atmospheric correction 
algorithm must be broadened and studied in depth. Secondly, our physical 
understanding of how certain typical canopy reflection patterns are correlated 
with measurable biophysical canopy parameters must also be deepened. To achieve 
this latter goal, we propose to build a simulated vegetative canopy structure 
using only artificial (non-living) structure elements for leaves, stalks, 
branches, etc., that provides a well-defined and fully reproducible testbed for 
verification and testing of instruments and theoretical models in remote 
sensing. Since every natural vegetation canopy changes through growth, and its 
canopy components and architecture are intrinsically variable, reproducibility 
is seriously limited and measurements of biophysical parameters, such as leaf 
area index and leaf angle distributions, are prone to large uncertainties. In 
contrast, the artificial canopy provides a facility for "ground truth" measure- 
ments, testing and calibration of instruments, verification, and sensitivity 
testing of theoretical models as well as specific instruments. A synopsis of 
this proposal is attached. 

If our preliminary findings about the invariance of typical local extremes 
in surface reflectance patterns to atmospheric perturbations can be further 
substantiated, then important conclusions should be drawn for the operation of 
future earth observing satellites. This result leads to the recommendation that 
off-nadir satellite observations in the near-infrared may contribute additional 
valuable information to crop identification as well as other scene identifica- 
tion when enough non-Lambert ian (e.g. specular) reflectance is present. For 
crop identification it might be sufficient to swing an MSS or TM type scanner 
into the downward sun direction to observe the hot spot, and then rotate the 
instrument around a 360° azimuth for the constant solar zenith angle. 
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Figure 3. Upwelling radiance distribution of coniferous forest at surface (upper half) compared with 
radiance at top of the atmosphere, corrected for the atmospheric perturbation (lower half). 
Vq = 23 km. 
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SCATTERING MODELS AND BASIC EXPERIMENTS 
IN THE MICROWAVE I’^GIME 

Pung# A.K, and A.J. Blanchard 
Wave Scattering Research Center 
Electrical Engineering Department 
University of Texas at Arlington 
Arlington, Texas 76019 


1. Introduction 

The objectives of the research over the next three years are; 

(1) To develop a randomly rough surface scattering model which is 
applicable ove^r the entire frequency band, 

(2) To development computer simulation method and algorithm to 
simulate scattering from known randomly rough surfaces, Z(x,y) . 

(3) To design and perform laboratory experiments to study 
geometric and physical target parameters of an inhomogeneous layer. 

(4) To develop scattering models for an inhomogeneous layer which 
accounts for near field interaction and multiple scattering in 
both the coherent and the incoherent scattering components. 

(5) Comparison between theoretical models and measurements or 
numerical simulation. 

Although various new approaches in rough surface scattering 
have appeared in the literature, they are generally restricted to 
either low or high frequency regions because' of the simplifying 
assumptions used. General theories which are not restricted by 
such simplifying assumptions are so complex that either no 
explicit solutions are available or only the scattered field is 
available because the average scattered power is too involved to 
compute. In order to obtain a general scattering model for 
randomly rough surface which is explicit and simple to use it is 
necessary to determine a sufficiently accurate surface current 
including just enough multiple scattering terms to have the 
correct behavior over the entire frequency band but not too 
involved to deny an explicit solution. This is our goal for a 
rough surface scattering model. In volume scattering problems 
most theories are restricted to either a sparsely populated medium 
or coherent scattering in a dense medium. Incoherent scattering 
which is of major interest in remote sensing has been investigated 
mainly by either the Born or the distorted Born approximation 
in the field approach and the radiative, transfer formulation in the 
intensity approach. Other methods are needed to include multiple 
incoherent scattering in the field approach and techniques are 
needed to determine the autocorrelation function in the coherent 
field calculation in dense inhomogeneous medium. Currently, the 
use of Percus Yevick approximation and the associated pair 
distribution function for hard spheres in coherent scattering 
computation leads to nonphysical results when the volume fraction 
of the medium is around 0.6. 

All theoretical scattering models involve some sim.plifying 
assumptions. To verify these models it is best to have 
controlled experiments so that target parameters can be varied and 
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any given experiment can be repeated. In the case of random 
surfaces it is possible to use computer simulation. This is not 
only cost effective compared with full scale field measurements 
but it is also more convenient to change target parameters than 
in any controlled experiment. 

2. Approach 

The approach to be used to develop the proposed random surface 
model is to use the integral equation method. We have used this 
method successfully to develop a scattering model for a perfectly 
conducting surface in the past year and we intend to generalize it 
to dielectric surfaces. To perform computer simulation for the 
rough surface problem we intend to investigate first the 
generation of a two- dimensional random surface with specified 
autocorrelation function. This is a major effort because 
algorithms and techniques for this purpose have been restricted to 
one-dimensional surfaces. From the numerical point of view this 
is also a large job requiring large capacity computer. Past 
experience on a one-dimensional surface indicates that it is 
desirable to have at least a 10,000 point surface to insure 
convergence of statistical parameters and it is desirable to have 
a 50,000 point surface to allow scattering computation at 
different frequencies. Once the surface is available the standard 
moment method will be used for scattering coefficient 
computations. Here again it is necessary to generalize from one to 
two dimensions and numerically this is a problem because we cannot 
allow the computation-size to increase like a square. 

The measurement program will be divided into two major 
tasks, (1) the design and construction of a versatile volumetric 
target structure, and (2) the measurement of radar backscatter 
cross section as a function of system and target parameters. The 
target is designed and constructed to simulate vegetation type 
structures. The target can be used to provide a variety of both 
geometric and electrical characteristics to include, leaf density 
in three dimensions, leaf orientation and location, leaf size and 
shape, and permittivity characteristics of the volumetric 
structure including relative permittivity and conductivity. The 
target is designed to allow a number of different parameters to 
be evaluated while maintaining the basic target shape. 
Measurements will include the effects of leaf distribution, 
shape, orientation and location on radar backscatter. Initial 
experiments wil 1 minimize the effects of the stalk structuresi 
however, future measurement programs could include these 
interaction effects as well. The measurements will be made in 
an indoor Radar Cross Section Range. The measurements will 
include amplitude and phase measurements at X band^K band, 
incident angles from 0 to 50 degrees, and the four linear 
transmit/ receive polarization states. 

In view of the difficulties with the current theoretical 
approaches as stated in the introduction to the volume scattering 
problem, we propose to use the eikonal approximation for incoherent 
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scattering coiaputation anA th» aolutton o( tbe roaster eg«ateion for 
deterrolfilng a roore realistic aafcocor relation function* l:t should 
be noted that .the pair Aiiifcr.lbiifcion function was derived for use 
in a dynaalc medium sucli as; a iiuld* Hence, for a weAluij; with 
random but rigid distr ifcdtions off particles It is not appropriate 
and does lead to nonphyalca.1. .results, 

3, Resulta to Date 

Over the past year we have d»v«:l0:p#d a scattering &odel for 
an inhomogeneous layer containing needle*- shaped leaves, We have 
made comparlsona between scatterinf fr©!® such a layer and another 
layer containing disc-ahaped leaves, The significant difference 
between the two types of vegetation is that the back scattering; 
angulat cu.cve for the disc-type has a hoiip in the mid angular 
region while the needle-typ®: has a dip in the same region when the 
sire of the leaves it comparable to the incident wavelength { see 
Plgares 1 and 2), In addiction, study has also been, carried out 
to detemine the effects of fresnel phase in the scattering phase 
matrices on s;catterers that are coHipatable tO: the incident 
wavelength bat are much smaller than tie wavelength in. at least 
one dimension. The conclusion la that c.he Inclusion of the 
Eresnel phase gives a more aceurate fregueaey trend dependence 
Csee lteerl984), Also studied is the effect of closer spacings 
between scatter ers accounted for In the scafetet-ing phase matrix, 
Results Indicate that s more accurate target ■parAin.eter can be 
inferred from such a model than Is possible using a standard 
radiative transfer pefchod (Fung and Eomi June, 1984 The inversion 
fros! Bieasuremersts for albedo and optical depth; was studied by Puncj 
and Bo.m {August rl:§84i with the conclusion that at. l.e.as;t a 
first-order theory meat be mmmS In the radiative transfe.r 
formulation for emiaaion ftpift an inhoiBogeneous layer. 

A scattering model fpr a perfectly conducting, rantdofflly rough 
surface has been developsa whic;h is valid over the entire 
freguency band CPanil984>, It has been demonstrated, analytically 
that the general soliition Bbtained reduces correctly fc.o the 
geometrical optics solution In the high frequency limit and to the: 
first-order small perturbation Bblufelon In the low frequency 
Halt, In the intermediate regies it has a behavior which can be 
very different from either limiting solution but is liable to be 
.mis. interpreted as one or the othe.c. The reason is that it could 
fb»:eipble the low frequency solution in polarisation behavior but. 
follow the general angular trend ■of the ;Klrchhoff solution {see 
ylgute 3). Also developed was the :SC;afct.e:ring model for a perturbed 
sinusoidal. ;au:rface { Eom and ;Pnng,l:9B4) , This was done to account 
for soil surfaces with row d.;»re,ct, tonal dependence* 

A plastic; model for a layer of leaves of either circular or 
elliptic shape has been designed which pe.rmits the Iteaf angular 
distribution, leaf size, leaf density, leaf dielectric value, and 
row directio-n. effects to be studied under laboratory conditlo.ns. 

It is expected that measurements conducted with this target will 
be able to provide reliable data for the v€;rl£lcafcion of existing 
volume scattering siodels. 
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REMOTE SENSING OF EARTH TERRAIN 


Jin Au Kong 

Massachusetts Institute of Technology 


1. Objective 


The objective of this research is to develop theoretical models that are useful and 
practical in the remote sensing of the earth environment including the earth terrain, the 
lower and the upper atmospheres. We have been very successful in developing various mod- 
els applicable to the microw’ave remote sensing of vegetation, snow-ice, and atmospheric 
precipitation. We shall extend such studies to the higher frequency range to unify the 
optical band and the microwave theoretical foundations. We shall also extend our study, 
which had an emphasis on vegetation canopy to include all terrain media, and the whole 
earth environment. A data base will be developed to generate scene radiation characteris- 
tics which will benefit the stxidies of global inhabitability, meteorological applications, and 
crop yield. 


2. Technical Approach 


The radiative transfer theory has been fully developed to treat the volume and rough 
surface scattering effects. The combined random rough surface and volume scattering ef- 
fects have been studied by employing a Gaussian random surface and applying the small 
perturbation methods which are modified with the use of cumulant techniques. The rough 
surface effects are incorporated into the radiative transfer equations by modifying the 
boundary conditions for the intensities. The T-matrix method, which utilizes vector spher- 
ical wave functions for the expansion of incident, scattered, and surface fields is used to 
derive the extinction matrix and the phase matrix for the radiative transfer equations. 
Densely distributed spherical dielectric scatterers is studied with the quantum mechanical 
potential approach. The quasi-crystalline approximation is applied to truncate the hier- 
archy of multiple scattering equations and the Percus-Ycvick result is used to represent 
the pair distribution function. The modified radiative transfer theory is studied with the 
Feynman diagrammatic representation. Atmospheric precipitation is also studied with the 
vector radiative transfer equations by making use of the Mie scattering phase functions and 
incorporating the rain-drop size distributions. To study scattering by anisotropic medium 
modelling sea ice and vegetation field with row structures, the dyadic Green’s function is 
first derived and used in conjunction with the first-order Born approximation to calculate 
the backscattering coefficients. 
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3. Research Results 


In active microwave remote sensing, a two-layer anisotropic random medium model 
has been developed to investigate the anisotropic effect of sea ice in which, due to the 
flow of sea water, the brine inclusions are tilted and to study the azimuthally anisotropic 
behavior of vegetation canopy which is attributed to the orientations of the vegetation 
stalks. In this scheme, the dyadic Green’s function for a two-layer anisotropic medium is 
developed and used in conjunction with the first order Born approximation to calculate the 
backscattering coefficients. It is shown that strong cross-polarization occurs in the single 
scattering process and is indispensable in the interpretation of radar measurements of sea 
ice at different frequencies, polarizations, and viewing angles. The effects of anisotropy on 
angle responses of backscattering coefficients are also illustrated. For passive microwave 
remote sensing of a two-layer anisotropic random medium, the principle of reciprocity is 
invoked to compute the brightness temperatures. The bistatic scattering coefficients are 
first calculated with Born approximation and then integrated over the upper hemisphere 
to be subtracted from unity, in order to obtain the emissivity for the random medium 
layer. The theoretical results are illustrated by plotting the emissivities as functions of 
viewing angles and polarizations. They are used to interpret remote sensing data obtained 
from vegetation canopy where the aniostropic random medium model applies. Field mea- 
surements with corn stalks arranged in various configurations with preferred azimuthal 
directions cire successfully interpreted with this model. 

The problem of microwave scattering from periodic surfaces is solved by using a rig- 
orous modal technique which conserves energy, obeys reciprocity, and takes into account 
the multiple scattering and shadowing effects. The theoretical results have been applied 
to the calculation of the radar backscattering cross sections in active remote sensing and 
the brightness temperatures in passive remote sensing, and used to match field data from 
soil moisture measurements. The angular behavior of the brightness temperatures has 
been explained with the threshold phenomenon by considering the appearance and dis- 
appearance of the various Floquet modes. The Kirchhoff approximation is used to study 
the scattering of electromagnetic waves from randomly perturbed quasi-periodic plowed 
fields. The narrow-band Gaussian random variation around the spatial frequency of the 
sinusoidal variation is used to introduce the quasiperiodicity. The physical optics integral 
is evaluated to obtain closed form solutions for coherent and incoherent bistatic scatter- 
ing coefficients. In the geometrical optics limit, it is shown that the bistatic scattering 
coefficients are proportional to the probability of the occurrence of the slopes which will 
specularly reflect the incident wave into the observation direction. The theoretical results 
are illustrated for the various cases by plotting backscattering cross sections as a func- 
tion of the incident angle. It is shown that there is a large difference between the cases 
where the incident wave vector is parallel or perpendicular to the row direction. When 
the incident wave vector is perpendicular to the row direction, the maximum value of the 
backscattering cross section does not necessarily occur at normal incidence. The scattering 
coefficients can be interpreted as a convolution of the scattering patterns for the sinusoidal 
and for the random rough surfaces. For the backscattering cross sections we observe the 
occurrence of peaks the relative magnitudes and locations of which are explained in terms 
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of the scattering patterns for sinusoidal surfaces. The combined random rough surface 
and volume scattering effects have been studied by employing a Gaussian random surface 
and applying the small perturbation methods which are modified with the use of cumulant 
techniques. The rough surface effects are incorporated into the radiative transfer equations 
by modifying the boundary conditions for the intensities. The radiative transfer equations 
are then solved numerically using the Gaussian quadrature method and the results axe 
illustrated cind compared with experimental data. 

The radiative transfer theory is applied to calculate the scattering by a layer of ran- 
domly positioned cind oriented nonspherical particles. The scattering amplitude functions 
of each individual particle are calculated with Waterman’s T-matrix method, which uti- 
lizes vector spherical wave functions for the expansion of incident, scattered, and surface 
fields. The orientation of the particles is described by a probability density function of the 
Euler ian angles of rotation and a rotation matrix is used to relate the T-matrix of the prin- 
cipal frame to that of the natural frame of the particle. The extinction and phase matrices 
for the radiative transfer equations cire expressed in, terms of the T-matrix elements. The 
extinction matrix for nonspherical particles is generally nondiagonal. It is shown that there 
are only two attenuation rates in a specified direction of propagation. The scattering of a 
plane wave obliquely incident on a half space of densely distributed spherical dielectric scat- 
terers is studied with the quantum mechanical potential approach. The quasi-crystalline 
approximation is applied to truncate the hierarchy of multiple scattering equations and 
Percus-Yevick result is used to represent the pair distribution function. While results at 
high frequencies are calculated numerically, closed form solutions are obtained in the low 
frequency limit for the effective propagation constants, the coherent reflected wave cind the 
bistatic scattering coefficients. 

In the strong fluctuation theory for a bounded layer of random discrete scatterers, the 
second moments of the fields in the second order distorted Born approximation axe obtained 
for co-polarized fields. The backscattering cross sections per unit area axe calculated 
by including the mutual coherence of the fields due to the coincident ray paths, and 
those due to the opposite ray paths, corresponding to the ladder and cross terms in the 
Feynman diagrammatic representation. It is proved that the contributions from ladder 
and cross terms for the co-polarized backscattering cross sections axe the same, while the 
contributions for the cross-polarized backscattering cross sections are of the same order. 
The bistatic scattering coefficients in the second order approximation for both the ladder 
and cross terms are also obtained. The contributions from the cross terms explain the 
enhancement in the backscattering direction. 

Both passive and active remote sensing of atmospheric precipitation are studied with 
the vector radiative transfer equations by making use of the Mie scattering phase functions 
and incorporating the rain drop size distributions. For passive remote sensing we employ 
the Gaussian quadrature method to solve for the brightness temperatures, and for active 
remote sensing an iterative approach carryied out to second order in albedo is used to 
calculate the bistatic scattering coefficients and the backscattering cross sections per unit 
volume. 
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4. Significance of Results 


The anisotropic medium model that has been developed is now able to account for 
the like as well as cross polarization measurements of the backscattering cross sections in 
active remote sensing of sea ice. The same model applies to the interpretation of data 
from the passive remote sensing of vegetation field. Extending this two-layer medium 
model to multilayered configuration will facilitate the interpretation for more complicated 
scene radiation characteristics. 

The strong fluctuation theory has been successfully developed to account for the 
volume scattering effects, not only for earth terrain media such as vegetation and snow-ice 
fields, but also for the remote sensing of atmospheric precipitation. The periodic rough 
surface theory provides a rigorous approach that satisfies both the principle of reciprocity 
and the principle of energy conservation. 


5. Next Major Steps 


There are three major steps: 1) To extend our study to include the effects of the upper 
atmosphere such as the troposphere, the ionosphere, and the magnetosphere. Our initial 
success with the usage of the random medium model applied to atmospheric precipitation 
serves as a solid starting point. 2) To extend our study to other earth terrain medium 
configurations such as canopy-covered soil moisture investigation, snow-covered fields, and 
layered-medium model including both volume scattering and rough surface scattering ef- 
fects. 3) To extend our study to unify the theoretical foundation for the microwave and 
the optical spectra. The starting point is the application of the radiative transfer theory, 
which has been developed in both frequency ranges rather independently. 
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1 . Introduction 

Microwave remote sensing of vegetated terrain has been studied. 
Vegetation is modeled so that backscattered radar signals can be used to 
infer parameters which characterize the vegetation and underlying 
ground. The vegetation is modeled by discrete lossy dielectric 
scatterers with prescribed characteristics. The goal of the modeling 
effort is to remotely sense vegetation type (classification), growth 
stage, and plant/qround moisture. This information can then be used as 
input into agricultural, forestry and global circulation models. The 
microwave frequency spectrum, particularly L and C bands, are especially 
appropriate for this purpose since the wavelength is comparable to plant 
leaf and stem size. The resulting resonant interaction leads to 
backscattered data highly depend on plant shape and orientation. In 
addition, the transparent nature of the atmosphere in this frequency 
regime allows for algorithm development which requires no atmospheric 
correction. 

2. Technical Approach 

The modeling approach characterizes the vegetated terrain as a 
layer of plant growth over a possibly rough surface. The plants are 
viewed as composites of leaves and stems which are modeled by dielectric 
discs (not necessarily circular) and rods. The discs and rods are 
placed at random locations throughout the layer; their size and 
orientation statistics are chosen to correspond to the vegetation being 
modeled. 

Once the discrete stochastic model has been specified, the average 
backscattering coefficients are computed for both like and cross 
polarizations. The calculation is made amenable to analysis by the 
presence of a small parameter - the fractional volume of vegetation. It 
is then possible to employ the distorted Born approximation or transport 
theory (depending on albedo size) to calculate the backscattering 
coefficients. This constitutes the formulation of the discrete problem, 
i.e., given the distribution of scatterer sizes and orientations, the 
backscattered power can be computed. The final step in the procedure to 
obtain a remote sensing algorithm is to invert this direct problem. The 
performance of the inversion assumes that the backscattering coefficient 
is given as a function of angle of incidence and/or frequency. From 
this data, the statistical characteristics of the scatterers are 
computed. This calculation is often numerically ill-conditioned and 
needs special treatment. 
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3. Research Results 

Over the three year period, results have been obtained in all areas 
needed for algorithm development. Electromagnetic modeling of plant 
elements has been addressed, backscattering coefficients have been 
calculated and inversion algorithms have been developed. 

Simple plant canopies such as soybeans can be modeled by replacing 
the plant's leaves and stems by dielectric discs and rods respectively 
as is shown in Fig. 1. The relative dielectric constant used for the 
discs and rods is the equivalent dielectric constant of the plant 
material. The use of this equivalent dielectric neglects scattering 
effects caused by internal variations within the plant. This is small 
however, because of the long wavelength involved. 

Scattering amplitudes have been computed for both the disc and stem 
elements. The resulting formulas for the scattering amplitudes are 
particularly simple because they take advantage of the small ratio of 
thickness to radius in the case of the discs and radius to length in the 
case of the rods. The simple formulation is particularly important 
since the results must be averaged over both inclination and size 
distributions. 

Backscattering from the modeled plant canopy has been treated by 
two techniques: the distorted Born approximation and the transport 

theory. The two methods are equivalent when the albedo of the 
scatterers is small and no coherent (planar) boundaries are present. It 
has been shown for the case of a flat ground that the transport theory 
neglects certain coherent interference terms which are taken into 
account by the distorted Born approximation. It should be noted that 
this effect is a low frequency (L-band) phenomenon because at higher 
frequencies the surface appears rougher, and as a result, coherence 
effects disappear. Since the distorted Born method is applicable in the 
L band frequency regime, since it contains the interference terms and 
since it is relatively simple in formulation as compared to vector 
transport theory, it has been used almost exclusively, to calculate the 
backscattering coefficients. 

Application of the distorted Born theory to the modeling of crops 
such as soybeans has yielded interesting results. Fig. 2 shows some of 
these results for a frequency of 1.5 GHz. The figure is a plot of 
versus angle of incidence. An examination of the figure shows that for 
vertical-like polarized returns the leaves are dominant at low angles 
and the stems are most important at large angles of incidence. This is 
not surprizing since the electric field is aligned with most stems at 
large angles of incidence. 

The model whose results are shown in Fig. 2 consists of leaves 
having radii of 2.5, 3.5 and 5.0 cm. Each of these leaf types has a 
density of 333/m . The stems are all taken to be 20 cm long and with a 
density of 1000/m . The leaf and stem size characteristics were 
obtained by the principal investigator who made on-site measurements at 
Beltsville . 

Although Fig. 2 just shows the vertiral-vertical backscattering 
coefficient, simple expressions were obtained for the 
horizontal-horizontal case as well as the cross polarized returns. 
These appear to be valid for frequencies below 4 GHz as long as the 
correct dielectric constants for the discs and rods are used. The de- 
Loor formula has been used for the calculation of disc and stem 
dielectric constants however, at lower frequency (= 1-2 GHz) a 

conductive contribution should be included. 
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The final step in any remote sensing problem is the inversion 
process. Results in this area have beeri obtained this year. Attention 
has been focused on the relationship between the backscattering 
coefficient and the joint probability density of disc radii and 

inclination angles. The expression derived by the distorted Born 
approximation has been used. An examination shows that it is nonlinear 
in nature. The relationship has been linearized in the 1-2 GHz region 
where the skin depth is large. The linearized expression (a Born 
approximation) is a Fredholm integral equation of the first kind. 
Inversion problems of this type are usually ill-conditioned and must be 
regularized. 

To simplify the inversion procedure further, it has been assumed 
that the radii and inclination angle densities are independent. If it 
is assumed that the density of one variable is known, then the other can 
be found through the integral equation with the knowledge of the 
backscattering coefficient for various angles of incidence. Such an 
inversion is shown in Fig. 3. Here the leaves have fixed radius of 4 cm. 
The inclination angles density is shown by the solid line in Fig. 3. 
The combined joint density is then used to generate 0 ^ 1 -, for various 
angles of incidence. The calculated values of the backscattering 
coefficient are corrupted by noise and a Phillips-Twomey inversion 
algorithm is used to compute the inclination angle density at certain 
discrete points. The results, which are good, are shown by small 
circles, squares and stars in Fig. 3. 

The Phillips-Twomey method has also been applied to the inversion 
of the complete two-dimensional density of radii and inclination angles. 
Here it has been assumed tht backscattering data is known at different 
angles of incidence and for different frequencies. 

4. Significance of Results 

The results just presented provide a direct relationship between 
the backscattering coefficient and the detailed parameters that 
characterize the scatterers. This means that by remote L band radar 
measurements at different angles of incidence and for different 
frequencies, such quantities as leaf inclination angle distribution and 
leaf area distribution can be determined. These distributions can then 
be directly used to estimate growth stage by knowledge of the area 
density of leaves and to estimate stress, due to lack of water, by 
knowledge of the inclination angle density. 

5. Future Work 

The model developed thus far has been directed toward leafy 
agricultural crops such as soybeans. By incorporating other scattering 
types of a more complex nature, other crops and forests can be modeled 
and inversion algorithms developed. Future scatterer development should 
include the modeling of non-planar leaves, branches and stalks with 
periodic nodes (corn). 

An experiment should be planned so that the theoretical 
developments can be tested. The experimental approach is particularly 
important in the inversion algorithm development since equipment 
limitations can force consideration of different inversion techniques 
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A vegetation canopy is a highly inhonjogencoos raedium at 
fflicrowava :f re-quoacie# , and bacaqae the- seatt-ering elements 
Cleaves* stalks* fruits* and branches) have a nonuniform 
distributidn in orientation*, the canopy is likely to exhibit 
non isotropic attenuat ion properties . In some canop ies , the 
stalk may contain the overwhelniing .«aj.ority of the plant’s 
bioiaass, which suggests that an iecideirit r;- Car wave would be 
differentially attenuated by the: canopy: depending on the 
d:irectiof} of the incident eieetrlc field relMiive to the stalks' 
or ie:ntation , The propagatian properties of a vegetation canopy 
play a cwntT a I role in njodeling. bo:th the backscattering behavior 
observed by an 5 mag 1 ng radar a.fid the .emi.ssian observed fey a 
radiometer. These propagation properties are in turn governed by 
the dielectric properties and the siEe, shape, and 
s loped i St r j hut i ons of tha s.cstterers * In spite of the crit ical 
need for canopy propagatioh irpdeXa and experimental data* very 
few invest tgat ions had been .cond.acted (prior to this study) to 
determine the extinct ton propercies: of vegetation canopies, 
eittier by constituent type (leaves* sta.Iks , etc.) or as a whole* 

A three- faceted approa.feh waS: undertake.n in this :stud:y 
consisting of comp:lifflentary iTtyest.i.gfetions with each providing 
answers to specifi.c aspects of 'this overall question. The first 
facet c-onsisted of diel;ectric measureircents that were cooduetee! 
over the 1-10 GHz band to d.etermine the va:ri.atio.n.s of the 
dielectric permittivity and loss factor as a function of water 
content and salinit.y. Dielectric fnixing fnodeis were developed 
and evaluated against the experijn-entai data for individual plant 
parts flea:ves* stalks, fruit) of several types: of plants (wheat*: 
corn* and soybeans). An example is shown in F.ig, 1. 

The- .second facet of the study con:sisted: of direct 
mea:S:«:r:: 0 mersts of eanopy attenuation. U:sing truck-mounted radars 
a:S fcransait.titig sources positioned above the. canopy, an.d: small 
receive antennas laounted on fiberglass ra.ils lunderTOat.h the 
canopy Con the so:il surface),: attenuation measurerfients were 
conducted as a function of incidence angle (relative to. nadir), 
mierO:wava frequency, and polarization c.onf i.g.uration » In 
addition, special experiments were conducted to evaluate the 
relative contributions, of the various canopy constituents to the 
total loss for attenuation) factor of the canopy. These included 
wheat deGapitation experiments., ip which mea.surcmerit s of the 
ea.nopy loss factor wefe. made before- .and after cutting off .and 
removing the wheat heads, and soybean def ol.lation axperlraents 
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which, provided coriaparifOH oC th# e'Ktinctlon loss due tO: the 
stalks alone with the total, loss tfua to the stalks and leaves 
together « ikcr:.p'rdi:ng to Pig.* 2., the' st-alk-s recount For the 
aiajority of the propagation los.s in a soytoeana canopy for 
vertically polariz-ed waves, .w.heTe.as the ieave..s: are the doroi-na-nt 
absorbers for horizontal polariaa.ti.ofi.* 

k canopy pr.opsgat ion mo-del was developed consisting of three 
terms! la) leaf terra that aecaunts lor afeaorpt.ion by randomly 
oriented; dlsc-'ihaped leay.es, (b) a .stalk tetin that accounts far 
absorption by a uniaxial crystal consisting of thin, vertical, 
lossy cylinders in an air background, and (c) a branch term that 
ac'coonts for the loss by a medium containing ra.ndomly -oriented 
need.l.e”S:hsped branches . The first and third terms are 
polarization-independent and vary wjth iricidencs .angle as sec €, 
The stalk term has a stmng doptotdnnce un both p?olariKatipn and 
Incidence .a.ngle. The mod*> ! wav foimd t --i be in good agreeme;n.t 
with experiinental obs‘»r v<it i'>nv. ''bjo of the key parameters in the 
raodel is the dieieccri. const anr ,it vegetation matarial, 

which was cnlculatof* using t hi- die)tc"rrtr r*>i>fing models developed 
in the first facet of f fi i s lovt-sf Jgat son. The combination of the 
dieleetric mixing rnodt' 1 an! c ho iU!u*py prof.sgat ion model provixies 
the tseans to cowimt s hi* ! art or uf t hr canopy if given the 

water coatent.® and voiMnic l ra^ticno f>i i h* canopy eofis:t.i.t.u.ents 
and the canopy hei^'ld . 1 hr srude! ^ far the depandenee ph 

iocidance a-ngle, plat i/nt ars J wav^' 1 £-na,t h , 

In the next facet of this study, we attempted to relate the 
r©sy;lts obtained above to the- remote serssing problem. In order 
to focus the study on the canopy itself and avoid problems 
associated wi th v«r i-*t iona jft f he so j I L<dckgrtnmd , it was. decided 
to cover the soil surface with s<*reen wire and make passive 
microwave, obsarvat i on.s a funcfittn of time over the growing 

cycle. The screen wire, while it ha.*; no discernible effect on 
the growth of the plants, acts .is; a conducting surfac:©. with a 
near-zero emis.s.lvity . Thus, soil tn«i!>ture variations exercise nop 
influence oh the canopy emiss^ion. This technique was :used. to 
investigate the var iat i on u! the t>r i >;ht ness temperature of the 
canopy with l.ook direction (relative to the row diraGtion ) , 
incidence an.g,le, and polarization conf iguratipn -for wheat, 
soybeans, and corn cAnopies:. A radiative transf.er model was used 
to relate the. brightness te.raperatuTe to biophysical parameters of 
the canopy through the csno'p.y lef:SS' faCt-or. Figure 3 ;shows that 
the experimental observations a.r.e in good agreement with the 
model predict ions . A second set of radiometer obser vafe.ions were 
made for canopies under nstura.l conditions (no .screen used:). The: 
radia:tive transfer raode.l w.as ext:ended' to iii:e-i.ude' a soil-mo:4s.ture 
term that accounts for bhe emiss ion: .and r.qf l;ect ton by the soil 
surface. Again good agreement was 'foy:nd between theory and 
exper irneat ., 

The maj:or:i.ty of miccow.ave sciaitc ring, and e.im;i.ss ion models: 
reported in the literature have .treated the canopy as an 
isotropic medium. This impl ies: t hat, the at tenuat ion. coef f i.c ient 
is indeperwient of polarization, inc. idesree angle, and look 
directiori: (in azimuth). The? present study has shown that in most 
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cases, none of the above assumptions are valid. Consequently, 
modelers will have to incorporate the anisotropic character of 
the canopy in their models, v.'hich is likely to lead to more 
applicable models Cor the backscattering and emission from 
vegetation than now exists. 

The work conducted to data has established the overall 
functional relations between the propagation coefficients of a 
vegetation canopy and its biophysical and geometrical properties 
The results, however, are based on a number of specific 
experiments that were conducted for a specific set of crops at 
specific stages of growth. The goal of the next phase of this 
study is to develop a general propagation model that may be 
applied to a wide range of canopy types and conditions. The 
experimental component of the study will be designed accordingly 
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Figure 1. Measured Dielectric Behavior of corn 
leaves. 
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Figure 2. Comparison of attenuation measurements 
made for soybeans canopy before and 
after defoliation for vertical and 
horizontal polarizations. 
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OBJECTIVES 

The objectives of this research are to identify the primary contributors 
to lO-GHz radar backscatter from various natural and man-made surfaces and 
objects, and to use this information in developing new and better models for 
the scatter. Vfhen the true sources are known for the scattering that leads to 
variation in intensity on radar images, the images (and sets of them) may be 
interpreted more meaningfully in terms of the variation of parameters of 
interest for science or application. For example, better interpretation of 
vegetation images may be possible for yield forecasting and stress detection. 

APPROACH 

Backscattering coefficient (o°) measurements from crops, soil, snow, etc. 
have been conducted for over three decades, but only a minimal effort has been 
focused on determining the sources of backscattering in those targets. Using 
a pulse radar, Graf and Rode f 1982] in Germany used a solitary fir tree to 
determine its major sources of backscattering. Ulaby et al . [1982], using a 
defoliation technique, conducted scattering measurements on several types of 
crops. The only previously reported extensive direct measurements of sources 
of scatter have been conducted by the Remote Sensing Laboratory at the Univer- 
sity of Kansas through this research in 1983 [Zoughi et al., 1985; Wu et al., 
1985] . The approach and the results of these experiments provided the primary 
guidelines for the 1984 experiments. 

To measure the relative backscattered energy from various constituents in 
a target (i.e., com cob, stalk, leaves), a radar system with a fine range 
resolution and a small footprint at the target range is required. Measure- 
ments conducted during 1983 utilized an available FM-CW radar system adapted 
so that its transmitted frequency swept over a wide band to produce a range 
resolution of about 11.0 cm. A focused parabolic-reflector antenna system was 
used to provide narrow bearawidths in both the azimuth and the elevation direc- 
tions; this resulted in a footprint of 15 x 18 cm at the range of 4 m. Asso- 
ciated with this system were relatively high range-sldelobe levels which made 
target detection ambiguous at times. Therefore, funding for a specialized 
radar system was requested and approved. Consequently, a short-range very- 
fine resolution FM-CW radar system was designed and built. This radar system 
(SOURCESCAT) also provides a range resolution of 11 cm, but its range-sidelobe 
levels are appreciably reduced. The same focused antenna system used in 1983 
was used with the SOURCESCAT. 

Targets examined during the 1984 growing season included corn, soybeans, 
wheat, alfalfa, short grass, tall prairie bluegrass and several types of 
trees. To determine the main sources of backscattering in these targets, 
constituent defoliation (i.e., for com and trees) and layer-by-layer defoli- 
ation (i.e., soybeans and tall grass) were implemented. All the measurements 
were conducted at incidence angles of 30° and 50°. 
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RESULTS 


For an individual corn plant, backscattered energy is mainly due to the 
leaves; scattering from the cobs and the stalk is negligible, as shown in 
Figure 1. For a canopy, echoes from the upper 1 - 1.5 m portion of the plants 
dominate the total echo. Energy returned from the ground is insignificant due 
to the two-way attenuation of the radar signal through the canopy. 

Soybeans are very lossy volume scatterers. The two-way loss through one 
row of plants ranges from 20 to 15 dB at 30° and 50°, respectively. Backscat- 
tered energy from the upper 30 cm of the plant dominates the returned sig- 
nal. The energy returned from this portion is about the same at both inci- 
dence angles. Layer-by-layer defoliation Indicates that the unattenuated 
backscattering from various leaves in the canopy is about the same. 

The backscatter from wheat is dominated by the heads in both the early 
and late stages of the growth. Backscatter from heads remains almost constant 
at different stages, but attenuation through the heads is considerably more at 
the early stages of growth due to its higher moisture content at this stage. 
This shows backscatter Insensitivity to moisture content but direct loss due 
to moisture content of the heads. 

Individual alfalfa plants of about 40 cm tall were observed with no 
defoliation performed. In the presence of the plants, the returns due to the 
entire plants dominated the radar signal (individual parts of the alfalfa 
plants are too small to resolve). The returns in this case were slightly 
higher at 30° than at 50° incidence angles. With the plants removed, back- 
scatter from the roots was observed at the depth of 5 to 12.5 cm beneath the 
ground. 

Two different tall prairie grasses (50 cm high) were examined. The first 
site was an undisturbed and natural site which contained the dead grass m.ater- 
ial from the previous year on top of the soil . The second site had this dead 
material intentionally burnt off the soil before the growing season. There- 
fore, the latter site contained less soil moisture than the former site. The 
results of our measurements were also in agreement with this fact. In both 
cases, a volume of grass about 15 cm thick and about 25 cm above the soil gave 
the strongest backscattered energy. This experiment was conducted in colla- 
boration with investigators from Kansas State University. 

For pine trees, the needles showed the strongest backscatter, and caused 
the strongest attenuation in the radar signal. Cones, although insignificant 
contributors to the total backscatter, exhibited more backscattering than 
attenuation properties. Figure 2 illustrates these results. 

Four types of deciduous trees were examined. Leaves were a strong cause 
of backscattering and attenuation. However, removing the leaves and keeping 
the leaftails and small twigs and branches reduced the backscattered energy 
very little. Therefore, their contribution to the total backscatter cannot be 
ignored. 

These results are in great agreement with the results of the 1983 experi- 
ments, but there is more confidence in these results because a greater number 
of independent samples were observed and the ambiguities caused by range 
sidelobes were much less. 
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niRcnssiON 


To obtain a complete understanding of the wave-target interaction pro- 
cess, experiments need to be conducted and models developed in each of the 
following areas ; 

(1) Target Dielectric Properties : To model a vegetation canopy, the 

dielectric constant of canopy constituents (leaves, stalks, fruits, 
etc.) as a function of temperature and water-volume fraction needs 
to be known. 

(2) Target Attenuation Properties ; For media such as snow, forests, and 
vegetation canopies, attenuation properties of the medium are essen- 
tial in modeling. 

(3) Target Voli’.me Geometry : Methods need to be developed to measure the 

three-dimensional statistical distribution of orientation and den- 
sity of scatterers within the volume. 

(4) Target Scattering Properties ; Extensive measurements need to be 
conducted on the scattering behavior of (a) targets under natural 
conditions, and (b) target constituents. 

Information from all of these four areas enter into a scattering model. 
The existing scattering models in many cases assume or estimate the contri- 
bution of some of these parameters to the whole model. To develop complete 
models, these assumptions and estimations must very closely coincide with the 
true state of these parameters. This can only be achieved through sound 
experiments in the above areas. 

By modeling the scattering properties of target constituents and using 
our experiment results to cross-check the models, we may be able to model 
correctly the scattering properties of the target as a whole. 

A method-of-raoments solution is proposed to compute the scattering from 
the Individual plant parts, which may be modeled as dielectric cylinders, 
ellipsoids or smooth-curved dielectric slabs of arbitrary shapes. The 
analysis is based on the LeVine- Schwinger integro-dif ferential equation, 
derived from the equivalent polarization current. The solution for the total 
field is obtained numerically. The far-zone approximation is applied to the 
integral equation. Given the total field inside the dielectric scatterer, we 
can calculate the scattered field and, therefore, the scattering cross- 
section. 

More concentrated experiments and research in this area are needed to 
answer all the related questions raised and those that may be introduced as 
the research progresses. Experiments on more natural targets and man-made 
targets (such as buildings) will continue in 1985, along with extensive 
efforts on developing scattering models. 
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ACTIVE MlCROSaVg PROPERItES OF WETATI08 CASOPIES 


Jack F, Farls 

Callforiila Inatltute of Technologjf/ 
Jet Propalslon Laboratory 


Judtlf feat ton . Potential users of radar Imagery need a better fundamental 
understanding of the capabl I it te« of radar systems for vegetation studies than 
past studies provide. One approach Is the use of theoretical mode 1 a to predict 
observable active microwovt properties of vegetation. This In turn requires 
accurate observations of backscatterlng coefficients and other active microwave 
properties In field research studies . The background document for the SRARC 
program eaphaslzus the need to relate electromagnetic parameter* to classical 
biophysical descriptors and to onderstand the role of polar l*at ion » especially 
ctoas-polarlitatlon* Gfoala « The broad goal of sy study is to Increase the 
understanding of the effects of canopy structure on the active tnlcorwave 
properties of vegetation canopies, with particular attention to polarlaatlon. 

Approach . In the first year, I studied aircraft radar 8cattero»eter 
.^asureisents of .corn and soybean; .f ields In Iowa to detemlne the propertlea of 
theae aature flelda at iL-band ahdi C-bapd fb .C») a:n,d at .sensdr Ipak” 

angle.B :fro« 5 to §0 d«|ress,. In the .second; year, :1. used the ^GlOud .Mo;<iel 
(Ref . 1) with, a :s«t of lC~ba.nd., fS c«) .and KAband (2 i») trUck-based. radar 
.scat ter OBeter llke^polarleatlon data for corn In .Kansas to .estend, the oodel .and 
to relate the model parameters to biophysical eharaetertstlcs of the corn* for 
the third and current year, I used a C~band dUal-polartaatlOB truck scatter- 
oaeter at the University of California Kearney Agrtcultura:! Center to measure 
the active microwave propert;tea of trees in orchards. Using the Cloud Model, I 
was able to: estimate not only the backs catte ring coefficient, but also the 
attenwation coefficient, the density of baekscatterlng cross sections, and the 
transiitttance. The research acitiyity lead to a new procedure In wslttg such 
truck-based radar scatteroaeter data for thick, volume scattering aedia. 


Results . 

First Year . In th*> first year, 1 found that mature corn and soybean fields 
could be best separated with the use of C-hand cross -polar! zat Ion (HV) data at 
TM degrees . No separation existed for C-hand 1 ike-pol.) rizat ion (HH I . At 
L-band, the separation of corn and soybeans was best for 1 Ikc-poI ari ration (HH), 
but was not as good as It was for C-bano HV data, Mo separation existed in 
L~band cross-polarl rat ion data, 1 observed significant row direction effects 
f or sensor look-angle near I h degrees when like-polarization oas used. These 
effec's were strongest at. L-band, No row direction effects existed for cross- 
polarization for either hand. Finally, I observed that the presence of wet soil 
and/or vet canopy due to rain reduced the separations of corn and soybeans and 
prodticed a significant increase in backsc.atter , especially near 1 b degrees . 

Second Year . In examining tlte Cloud Model , I found tb he model form was 
significant 1 y different for cfoss-polarizat ion than If, .s for its usual like- 
polarization foriaulation. Also, I found that the use of the ratios of back- 
scattering coefficients for the various polarization combinations (VV, HH, HV 
and VH) could lead to the Isolation of canopy angular orientation , ropercies . 



Th 1 B al Ir>wed a new interpretation nf the corn and soybean data studied in the 
first year. -Similar separattonn existed In the dcpolarlzst ! on ratio (cross” 
polarization divided by 1 Ike-polar } aat ten ) such as existed Independently in the 
Setter of the two chaniieLs, Finally, by applying the Cloud Model to through-- 
the-season, truck -based radar scat teronet er nteasureiitetits , 1 could estimate the 
baekscattering and ext Inc.t ion cross sections of average loaves, 1 found that 
the former was a simple power-law funet 1 on of corn leaf area {f*" of D.9?) 
tbronghont the season , 

Third Year . J used the Microwave Scat terometer C~band (MSC), on loan from the 

Hf SA Johnson Space Center, mounted on a JFL t ruck (height of II m ) to make field 
aeasureaents. This tnstrument obtains calibrated radar backacattering coeffl- 
e.lent data in the field, however , It samples a volume that is limited in size 
due to the finite range-resolution and the finite beam-width. It obtains data 
at both like-polarization and croas-po lari ration. Through the use of the rang- 
ing capability of the MSC and the Cloud Model, T was able to estltnate the active 
a-tcrowave properties of vegetation. :cap,<sples, m follows! (. 1 ) the troe, 
eorrected backseattering eocf f Iclent, S”,|.i ( 2 ) the volume, extinction coefficient 
K <»“**')} (3) the hackscattering cross section density, E ( 4 ) the two- 

way transtnlttance of the canopy, t, t. {1 and f stand for H Or V to denote the 
polarlratloh :Coinl»inatlon «s,ed):|, 0 j tne canopy” backseat ter log coeffl- 

.clent;, .S®ipgj 'and '{€.)■ the, dtstr.lbat,.lon 'Of tottr<!.e:a' with, range,, FOt' a, glveft- 
vegetation. c:a,no:py' these, aetive: a.ic,rowav* chafacteristlcs: ■vary with :poia,rl:*atlon, 
or polartzatton oomTsioatioa.* 

Illustration, of the Ifew Proceditre . the pracedore to ohEal,n these pataaeters. Is 
illwstrafed below for an. Ideal case where the canopy^ properties are homogen- 
eous. One view of the Cloud Model la that the ohaerved backseat teeing eoeffl- 
ctent Is given by 


R3 


H 

*'3 

dS 

m 

If j dE ‘ 

f: /fg dR 

«I 


»i 

«2 


for the one-layer case, where; the canopy elements are found between the radges, 
8j and 8j (a), with the reflected (mirror Image) canopy extending to raoge 
83 Co), 


fj = cos (T) exp I-CiCj + K^) (R - Rj)I (2) 

^2 ” ^'IJ *"lj ^*2 ~ 

. exp hCK^ + Kj) (8 - 8j)l 

where T 1,8 the sensor look-angle (degrees), r| . is the surface reflectance 
Cassuaed to specular In baste character) for polarlxatlon IJ, and R Is the raage 

(m). 
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The Integration of Kq , (1) over the entire range frets 8| to 83 with the 
expressions In (2) and (3) leads to the true canopy backseat taring coefficient 


S, , cos (T) (1 + r . . t, t . ) (1 
ij Ji 1 J 

+ t, t . s , , 

i J Ijs 


ti t ) 


m 


where !s the surface backscatterlng coefficient . 

However , for the MSC (or any scatteroneter having a small, finite range resolu- 
tion ), the sensor measures only a fraction of the total backscatterlng. I call 
this the partial backscatterlng coefficient (measured ) which is related to the 

fuactloB f by 


«* |f dll Coy«r the range resotutinS, M) t:5)' 


Since AR Is small it i» possible to estittate the function f by AS®/AR. 

Figure 1: shows the comparison between the given, (true) wine of f and the 
estimated value of f as obtained from a modeled prediction of the "tteasured" 
piTtisl badkscatf.etini :cocfflcle,it, (see 'Fiiure S). The estimated f la close to 
the given ,p for some ran|eS' but .not .for all 'ranges:. Sow,: It (Is 'posslhie 'tO' 
,pe.r:(:or« a numerleal, Integeatton. fco «s:tl»a,:te the t:rae' valne: of teckscatterlng 
co.e'f:fl'C:4e»t, as l.ndlcat«d la, .Eg.. {::l) 'abov*:, To,r 'thick 'vege't'atlon :ca'ttdples, the 
true vailue :ls slgnlflcanfcty higher than the largest :partia:i ■vald«:,. 'Thus, this 
procad«:re Ss required .for sach #itu«t:tons. 

Since the variation of f with range is a decaying exponential function, the use 
of the logarithm of f (or alter natively, expression of f :in deClBels, dB) leads 
t :0 a straight line poirtion of the f (dB) versus R curve (see Figure 3). Massing 
a straight line through points In. this segment allows estimatlloh of both' E and :K 

by 


fi “ fjj set (T) and K « slope/ (■•'8,686) (6) and (7) 

1 

kith the >’<«t Imated va Ine of K and the known ilant range thickness at the canopy , 

t| f , can he est Itnated. Using rough estimates of r and surface bac'k.scatterl ng 
roef f 1 rlent , one can estitBat» lndepend'»nt ly the value for E, ITiIs value should 
be close to that estimated for Eq, (6), With the new values for K and K It is 

{'■osstble to estlisate the "lull caiuipy" backscatterlng coeffieleot from 

S°yg - E cos (T)/(Kj - Kj) (8) 

If the second (reflected ) portion of the f (dB ) versus R curve can be detected , 

it l« possib e to extrapolate the two straight lines (one f roa ranges less than 

the range to the surface, and one for ranges greater than the range to the sur- 
face) to the surface range values of ( left side) and (right side) and 



calculate r by %/ft « Since the estlaiated value of f at R 2 «ay be higher than 
the average value of and f g due to the contribution of attentuated surface 
bacfcecattettng, It is possible to estimate the surface backscatterlng coeffl"* 
eient,.. 

Results with: Orchards . Figure 4 shows the partial backscatterlng coefficients 
versus range for a young peach orchard (see Figure 5 which Indicates the tree 
structure). Note that near the top of the canopy the radar return Is doalnated 
by the W where HU dominates the return near the base of the canopy, After the 
noise floor Is estimated and subtracted, the data were used to estimate f CdB) 
as show In Figure 6. The active radar parameters shown In Table 1 are obtained 
from these data using the above procedure. 

Significance of Results. The three year study shows chat, in many cases, cross** 
pol.'irlEatlon or a combination of H8, W and/or cross-polar iaat ton yellds infor- 
mation on the angular orientation of vegetation canopies. Ratios of channels 
should be used ttj Isolate such lafomsatioii from other vegetation characteriat f cs 
stjci! as total areal biomass or water content. Furthermore, with the proper siso 
ot sinpie microwave scattering models, such ;as the Cloud Model, and acourat*-' 
radar scot tc rometer data, It Is possible to eslt*ate the total backscatt or J ng 
coefficient, the attentuatlon coeffieleKt , the transalttance , and average near- 
terlng and ext S notion radar cross sections of scattering elements. These darn 
should nr.'ive ascrisl for wertficetton of aodellng approaches and for the develop- 
ment nf edu.jtiate iinderstanding of the info r motion content of ra d-ir serujor data 
such as From the svntbettc aperture radar (SAR) laagers, 

Se Fe.'ences . 

i. Attcma, K ,P, W. and F»T, Hlaby, 197St Vegetation aodeied as a water 
.'I'ltid. Radio Scl . 13C2)t357-364, 


Table 1, 


Example of Measured Active Microwave Properties 
of Orchards Through Analysis of Truck Radar 
Scatterometer Measurements at C-Band and 60 Degrees 


Field' 

Pol 

^ max 

(dB) 

E 

(m~b 

K 

(m~^ ) 

tj 

,CdB:) 

Mature 

?V 


-10.2 

0.077 

0.199 

0.021 

-10,1 

Peach 

HH 


-11.6 

0.050 

0.178 

0.031 

-11.5 

Orchard 

VH 


-17.3 

0.014 

0.188 

0,026 

-17.2 

Young 

VV 


-10.4 

0.036 

0.082 

0.203 

MM 

Peach 

HH 


-12.1 

0.028 

0.J87 



Orchard 

VH 


-15,8 

0.011 

0.084 
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LIST OF PUBLICATIONS AND PRESENTATIONS 


The following pages list the refereed journal articles and 
other publications and presentations that report on the results 
of the research performed under the SRAEC Project. The listings 
are presented as provided by the principal investigators and are 
arranged in the same order as the summary reports presented in 
Sect ion 1 1 . B . 

SRAEC Project scientists have produced approximately 80 
refereed scientific journal papers during the first three years 
of the program as well as dozens of additional writings in the 
form of theses, dissertations, symposium and conference proceed- 
ings, and special documents. 
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tio. f submitted. 29 ms. pp. t 7 fig. + 3 tab. 
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Geoscience and Remote Sensing , submitted. 54 ms. pp. + 
16 fig. + 4 tab. 

Strahler, Alan H., C. B. Noodoock, and J. A. Smith, 1985. 
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on Geoscience and Remote Sensing , submitted. 42 ms. 
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Proceedings 

Smith, J. A., K. D. Cooper, and A. H. Strahler, 1984. 

Electro-optical terrain reflectance modeling: A per- 
spective. Proc . SPIB Tech . Symp . . Spr . 1984 . pp. 12- 
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Strahler, Alan H., X, Li and J. Franklin, 1985. 
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t 
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scenes, Proc . Eighteenth Int . S ymp . on Remote Sensing 
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press, 10 pp. 
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Georaetric~optical Modeling 
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Pat .te rn , o£. Coniferous Trees . M. A. Thesis, Dniv. of 
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fornia, Santa Barbara, March 1985, 78 pp. 
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CROP CHARACTERISTIC RESEARCH: 6. D. Badhwar 

GROWTH AND REFLECTANCE ANALYSIS Johnson Space Center 
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PUBLICATIONS 

Journal Articles 

1. G. D. Badhwar, W. Verhoef, and N. J. J. Bunnik: Comparative Study of 
Suits and SAIL Canopy Reflectance Models, Remote Sensing of Environment, 
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A member of AgRISTARS documents describing the data produced. 
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ATMOSPHERIC AMD SCENE COMPLEXITY EFFECTS ON SURFACE BIDIRECTIONAL REFLECTANCE 


David J. Diner (Principal Investigator) 
Refereed Journal Publications 


1. Diner, D, J. and Martonchik, J. V» (1984). Atmospheric transfer of radia- 
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J. Quant. Spectr. Rad. Transf. 31 , 97-125 • 
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MULTIDIMENSIONAL MODELING OF ATMOSPHERIC EFFECTS 
AND SURFACE HETEROGENEITIES ON REMOTE SENSING 
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4. A. Zardecki, S.A.W. Gerstl, and J.F. Embury, "Multiple scattering effects in 
spatial frequency filtering," Appl. Opt. 4124-4131 (1984). 
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A8. Radiative transfer tbeory for active remote sensing of a layer of nonspberical par- 
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Lee and J. A. Kong), IBBE Transcations on Geoscience and Remote Sensing, sub- 
mitted for publication. 

C Q B,fgr ?ncs .Ariidg8 

CL Tbeoretical models for microwave remote sensing of vegetation and soil moisture 
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DISCRETE RANDOM MEDIA TECHNIQUES FOR 
MICROWAVE MODELING OF VEGETATED TERRAIN 
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Approach, with O.S. Sidhu, IEEF-GE-21, p. 62-71, 1983. 
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Particles, with Seker and LeVine, to be submitted to Radio Science, 1985 
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with LeVine, presented at URSI, Boulder, CO, 1983. 

Mean Wave Propagation in a Trunk Dominated Forest, with A. Schneider and 
S. Seker, presented at lEEE/URSI Symposium, Houston, TX, 1983. 

Backscattering from a Vegetation Layer-Coherent and Incoherent Models, 
presented at URSI Commission F Symposium on Microwave Signatures in 
Remote-Sensing, Toulouse, France, 1984, 
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Scattering from Thin Dielectric Disks with D. LeVine - presented at 
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Area Distributions, with H. Saleh, to be presented at lEEE/APS URSI 
Meeting, Vancouver, Canada, 1985. 
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MEASURING AND MODELING OF THE DIELECTRIC 
PROPERTIES AND ATTENUATION OF VEGETATION 

Fawwaz T. Ulaby 
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IEEE Trans, on Geoscience and Remote Sensing , GE-21, N.l, pp. 
51-61, January 1983. 


Brunfeldt, D. R., and F. T. Ulaby, "Microwave Emission from Row 
Crops," submitted to IEEE Trans, on Geoscience and Remote Sensing . 


Ulaby, F. T., and E. Wilson, "Microwave Attenuation Properties of 
Vegetation Canopies," submitted to IEEE Trans. Geoscience and 
Remote Sensing SRAEC Special Issue. 


Allen, C. T., and F. T. Ulaby, "Modeling the Polarization 
Dependence of the Attenuation in Vegetation Canopies," submitted 
to IEEE Trans, on Geoscience and Remote Sensing. 


Symposia 


Jedlicka, R. P,, and F. T. Ulaby, "Measurements of the Microwave 
Dielectric and Attenuation Properties of Plants," National Radio 
Science Meeting (URSI Commission F), University of Colorado, 
Boulder, January 5-7, 1983. 


Brunfeldt, D. R., and F. T. Ulaby, "Measured Microwave Emission 
and Scatter in Vegetation Canopies," IEEE International 
Geoscience and Remote Sensing Symposium (IGARSS '83), San 
Francisco, California, August 31 - September 2, 1983. 
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Ulaby, F. T., "Microwave Properties of Vegetation Canopies; An 
Overview," International Colloquium on Spectral Signatures of 
Objects in Remote Sensing, Bordeaux, France, September 13-16, 
1983. 


Ulaby, F. T., R. P. Jedlicka, and C, T, Allen, "Measurements of 
the Sources of Scatter from Vegetation," Microwave Signatures in 
..aote Sensing (URSI Commission F), Toulouse, France, January 
3-20, 1984. 


Brunfeldt, D. R., F.T. Ulaby, and M. C. Dobson, "Active and 
Passive Observations of Vegetation Canopies," Microwave 
Signatures in Remote Sensing (URSI Commission F), Toulouse, 
France, January 16-20, 1984. 


Allen, C. T., and F. T. Ulaby, "Modeling the Polarization 
Dependence of the Attenuation in Vegetation Canopies," IEEE 
International Geoscience and Remote Sensing Symposium (IGARSS 
*84) , Strausbourg, France, 27-30 August, 1984. 


Brunfeldt, D. R., and F. T. Ulaby, "The Effects of Row Direction 
on Microwave Emission from Vegetation Canopies," IEEE 
International Geoscience and Remote Sensing Symposium (IGARSS 
*84), Strasbourg, France, 27-30 August, 1984. 
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PDBLICATIONS AND PRESENTATIONS 
Determination of Backscattering Sources in Various Targets 
R.K. Moore, Principal Investigator 


1 « Refereed Journal Articles 


"Determination of Backscatter Sources in Surface Type Targets," R. Zoughi, 
R.K. Moore, L.K. Wu, F.T. Dlaby and A. Afifi, to be published in 
International Journal of Remote Sensing, 1985. 

"Preliminary Results on the Determination of the Sources of Scattering from 
Vegetation Canopies at 10 GHz," L.K. Wu, R.K. Moore, R. Zoughi, F.T. 
Dlaby and A. Afifi, to be published in International Journal of Remote 
Sensing, 1985 . 

"SOURCESCAT; A Very-Fine-Resolution Radar Scatterometer ," R. Zoughi, L.K. Wu 
and R.K. Moore, submitted to Microwave Journal, 1984. 


2. Conference and Symposia Presentations 

"Measurements of the Sources of Scatter from Vegetation," R.K. Moore, R, 
Zoughi, L.K. Wu, A. Afifi and F.T. Dlaby, presented at the IIRST 
Specialists Symposium, Toulouse, France, January 1984. 
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APPENDIX A. THIRD ANNUAL PRINCIPAL INVESTIGATOR WORKSHOP — 
INVITATION LETTER AND AGENDA 
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National Aeronautics and 
Space Administration 


Goddard Space Flight Center 
Greenbelt, Maryland 
20771 


Reply to Attn of: 623 



Dear Colleague: 

Trivia question No. 1: "Which state of the original thirteen was 

the sixth admitted to statehood in 1788 and has a land area of 8,257 
square miles?" Trivia question No. 2: "Which city in that state 

has well-known institutions of higher learning that were established 
In 1636 and 1861?" Trivia question No. 3: "When and where is the 

third annual ‘^RAEC Pi's meeting going to be held?" The annual 
ga^-hering of the NASA Fundamental Remote Sensing Science Research 
Program's Scene Radiation and Atmospheric Effects Characterization 
Project's principal investigators is not a trivial event, this 
year especially, as each of you will be reporting the results and 
significance of three years of research I This letter is to 1) 
confirm the time and place for the meeting, 2) advise you concerning 
follow-up communications concerning the meeting, and 3) inform you 
of the preparations that you should make for this meeting. 

As you have been informed via telecon, the third annual SRAEC Pi's 
meeting will be held on the campus of the Massachusetts Institute 
of Technology in Cambridge during the last three days of January 
(29, 30, and 31) 1985. Dr. Jin Kong has graciously offered to host 
this meeting and is making the necessary local arrangements and 
finalizing the agenda for the meeting. Dr. Kong will be sending 
you further details concerning these items within the next few 
days. It might be useful to mention at this time however that a 
special hotel room rate has been negotiated and a block of rooms 
have been reserved for our group at the nearby Sonesta hotel. 

I'm very pleased that at this time it appears that we may have a 
100% P.I, participation at this meeting. I want to encourage you 
to make every effort to attend the entire three days of the meeting. 
We are planning to take a group photo at this meeting and I would 
like everyone in the picture. More important, however, is that we 
want full participation in the results reporting and research issues 
discussions. 
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This meeting, and your individual reports, will be focused on summarizing 
the significant contributions that these past three years of fundamental 
research have made to the remote sensing scientific community. In addition 
to learning about the more technical details and specific progress that has 
been made in your more focused research activity, we all together need to 
be able to summarize the new insights that have been gained, the remote 
sensing problems that have been solved and, in general, what benefit this 
project has provided to NASA and the remote sensing and broader scientific 
communities. We need to collectively address a few key research issues and 
what are the new research priorities, based on what we have learned. Thus, 
we need to spend some "quality time" that will help both direct and defend 
this program. 

In addition to bringing your body and nrind to this meeting, I need you to 
bring with you a brief written summary of your research that will be compiled 
with the summaries from all of the Pis into an annual report. This year I 
am asking for a single written report that falls somewhere between the 
"Investigation Summary Results Abstract" and the "Investigation Synopsis" 
that you prepared last year. I need a succinct report that is introduced 
by a paragraph on the justification, goals and objectives for your research 
(the broader two or three year plan). This should be followed by a paragraph 
describing your technical or experimental approach, A discussion of the 
research results should then follow and can use up to three key figures to 
illustrate your findings. The final one or two paragraphs should highlight 
significance of your results and (briefly) what the next major step 
should be in the logical progression of your research focus. This report 
should be a minimum of two and a maximum of four single-spaced typewritten 
pages, excluding figures. Figures should be suitable for reproduction 
(clearly readable) at one-half of full size as submitted to me on 8 1/2" x 
11" paper. As the heading at the top of the first page you should include 
the descriptive title of your study (all capital letters), your name (first) 
and then any co-investigator's names, and your institutional affiliation 
beneath the names. 

On a separate page (typed single-spaced with a double space between entries) 
you should include your list of publications and presentations that have 
resulted from your funding from the SRAEC Project. You should give your 
study title and your name (as PI) only at the top of this page followed by 
a separate listing of 1) refereed journal publications and 2) symposia and 
conference presentations. 

Please adhere to the limits given above in preparing your report. Remember 
that this report is not a journal paper but rather is intended for a less 
technical readership (i.e., scientists from other disciplines and research 
program managers) and should be written to provide an interesting overview 
of your work. Your journal papers are available for anyone v/ho needs or 
wants all the details and highly specific/technical results. Help me put 
together a report that will enable us to keep this program alive and well 
within NASA during this era of significant government budget cutting. 

Your oral report should be prepared as last year, allowing 15-20 minutes 
for presentation and 10 minutes for discussion. One final topic concerns 
the IEEE Transactions on Geoscience and Remote Sensing "Special Issue on 
Scene Radiation and Atmospheric Effects Characterization". For various 
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reasons the communication concerning the deadlines for this special transactions 
has not been what It should have been. Nevertheless, the current "deadline" 
for you to submit five copies of your paper (to Bob Murphy) for the September 
1985 specical Issue is December 15, 1984 (moved back from the original 
September deadline). However, I have spoken with the editor and he has 
Indicated that they can be submitted later than this date but the probability 
of their being reviewed, revised and re-submitted in final form in time 
for this special issue decreases with each week past this deadline. Probably 
early January will be sufficient time. Please submit your paper as soon 
as possible, as we would like a paper from each one of the Pis In the 
SRAEC project. 

Thank you for your eagerness to participate and future help in making this 
year's meeting a success. If you have further questions please feel free 
to call me at (301) 344-9186 or Jin Kong at (617) 253-5625, as appropriate. 

I look forv/ard to seeing you again In January, 

Sincerely, 


/ 



Donald W, Deerlng 
Earth Resources Branch 
Laboratory for Terrestrial 


X 


\ 


\ 


Physics 


P.S. Answers to Trivia Questions: 

1. Massachusetts 

2. Cambridge (Harvard, MIT) 

3. January 29-31, 1985; Cambridge, MA 
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Ms.gtmg„^g9nda 


Third Annual Meeting of the NASA Fundmiental Research Scene Radiation 
and Atmospheric Ejects Characterisation 


January 29-31, 1985 


Tuesday, January 29, 1985 


09:00-09:20 

09:20-09:50 

09:50-10:30 

10:30-11:00 

10 : 00 - 12:00 

Lunch 

13:00-15:00 

15:00-15:30 

15:30-17:00 

18:00 


MIT Welcome {Jon Allen) 

NASA Welcome and Meeting Objectives (Don Deering) 

Invited Speaker Ken Hardy on applications of satellite data for me- 
teorological analysis 
Coffee Break 

PI Reports (Gerstl, Diner) 

PI Reports (Fraser,Pee>’ce,Badhwar, Vanderbilt) 

Coffee Break 

PI Reports (Kimes,Smith,Strahler) 

Reception, MIT Faculty Club 


Wednesday, January 30, 1985 


09:00-10:30 

10:30-1 0:45 

10:45-12:15 

Lunch 

13:00-15:00 

15:00-15:30 

15:30-16:00 

16:00-17:30 


PI Reports (Norman ,Paris,U!aby) 

Coffee Break 

PI Reports (Fung,Kong,Moore) 

Photo opportunity 

PI Reports (Lang,Conel,Irons,Kahle) 

Coffee Break 

Canopy model inversion (Goel) 

Presentation and discussion on proposed artificial canopy experiment 
(Chair: Gerstl) 


Thursday, January 31, 1985 


09:00-09:30 

09:30-1 ();00 
10:00-10:30 
10:30-11:00 
11 : 00 - 12:00 

Lunch 

13:00-14:30 

14:30-15:00 

15:00-17:00 


Invited Speaker Dave Staelin on Retrieval Techniques for Atmo- 
spheric Profiling 

Report on ISLSCP-FIFE (John Norman) 

NASA Headquarter Management Perspective (Bob Murphy) 

Coffee Break 

Discussion on project accomplishments and future research direc- 
tions (Chair: Smith and Ulaby) 

Research topic discussions and working papers 
Coffee Break , . 

Written summary for project accomplishments and future research 
recommendations 
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Msg: XGIF-1963-1511 


Posted: Thu Jan 31, 1985 10:28 AM EST 
From: REMURPHY 

To: DDEERING 

Subj: RE: Message to SRAEC Pis 

DON, Please read this note to the troops for me. 

I am sorry that I am not able to be with you at the 4th (can it really 
be the 4th?) meeting of the SRAEC principal investigators. My new job has had 
me on the road more than half of the time. I just couldn’t do this one, both in 
terms of the backlog at my desk and the impact on my family. 

I just returned from Geneva where ! secured the blessing of the World 
Meteorlogical Organization for the International Satellite Land Surface 
Climatology Project (better known to its progenitors as ISLSCP). This was a 
major step forward in the process of orienting NASA’s program to scientific 
goals. There is a description of ISLSCP activities in the current Space Science 
and Applications Notice which you should all have received. If you didn't see 
it, Don Deering, Bob Fraser and John Norman know everything you need to know 
about ESLSCP. There are also several reports by working groups which are 
available from Ron Witt (for retrospective studies) or Piers Sellers (for the 
field e?{periraents). Both are at Goddard. That program should provide some of 
the focus for the application of the your SRAEC efforts in the future. (You may 
recall that 1 said that cultural crops where not the focus anymore. Here is one 
of the alternatives.) 

We are attempting to generate some other foci for the program. Berrien 
Moore of the University of New Hampshire is chairing a subcommittee of the 
Earth System Sciences Committee which deals with the Land Processes. We had a 
very successful meeting at UNH in November, and Jim Smith is organizing a 
follow-on meeting nest month at Goddard to discuss measurement and modeling 
strategies. The Berriens committee meets again in March at Oregon State. Out of 
these meetings 1 expect some more fruitful directions. 

The work of the SRAEC group will be at the center of all that we will 
be doing. 1 hope to have a meeting of all Pi’s in the total program sometime 
later this year so that all can benefit from the type of interaction you are 
having now. 

A not quite pleasant topic to bring up is the IEEE special issue. A few 
of you (and you know who you are!) still owe me your papers. Soon! (please?) 

What we have so far looks great. 

To return to pleasant things, I will be looking forward to receiving 
your renewal proposals over the next few months. There is much left to do, and 
you guys can do it. If you have questions about the nature of desireable work 
you can talk with Don or me. I will conduct the peer review from NASA HQ, but 
will do so in close coordination with Don and I will look to him for 
programatic guidance. 

As for the future of the Land Processes Program, I think that it looks 
good. This year we are (I hope) at the nadir of our funding with SMS less than 
fiscal year 84. Total useful funding this year is about 12 MS for Geology, 
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Hydrology, Terrestrial Ecosystems, Remote sensing Science (thats you) and 
Applications. For next year, we have received high priority within NASA for an 
augmentation. Listen to President Reagan’s budget message next week and see if 
he mentions us by name (HAS). 

If we do receive that augmentation their can be some improvement in our 
average funding levels as well as the addition of new activitcs and some 
further upgrading of facilities. Things like the artificial canopy can proceed 
out of our existing funding base if someone redirects their curent work if we 
don’t get the augmentation, or it can be an add on if we get it. (This, of 
course, asumes that it passes muster in terms of the programatics and peer 
review. Your input through Don from this workshop will shape the programatics; 
peer review will be according to pur usual capricious standardsl) 

As you know, we had a successful shuttle flight this past year. While 
we got only a fraction of the SIR-B data that we hoped for, there is much good 
science ahead in the analysis of it. SIR-B will be reflown in early 1987 and 
SIRC will follow about a year later. We are hopeful of flying the Shuttle 
Imaging Spectrometer Experimet (SISEX) also in the 80’s. Our aircraft 
capabilities are gradulaly being upgraded and we know have semi- operational 
bidirectional measurement capabilities from the air as well as the Airborne 
imaging Spectrometer (AIS) , the Thermal Imagging Mapping Spectrometer (TIMS), 
a quad-pol SAR and a 4 beam L-Band pushbroom radiometer. 

I know I’ve left somethings out, and I am typing this while I am in 
transit to yet another meeting. So please excuse the rambling and incomplete 
nature of it. Keep up the good work. I hope to see you next year if not sooner. 
Perhaps at the IGARS meeting at UMASS. 


Bob 
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